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Summary 


An equipment is described used to establish automatically the statistical distribution of 
the time intervals during which the sound level of an acoustical phenomenon exceeds a 


prefixed level. 


Some results are reported that have been deduced from the examination of traffic noise 


and from the analysis of different voices. 


Sommaire 


On décrit un appareil permettant de déterminer automatiquement la distribution statisti- 
que des intervalles de temps pendant lesquels le niveau sonore d’un phénoméne acoustique 


est supérieur 4 un niveau fixé 4 l’avance. 


On donne quelques résultats déduits de l’examen des bruits de circulation et de l’ana- 


lyse de différents types de voix. 


Zusammenfassung 


Es wird ein Apparat fiir die automatische Aufnahme der statistischen Verteilung von 
Zeitabstanden beschrieben, waihrend welcher eine akustische Erscheinung einen héheren 


als den vorher festgelegten Schallpegel hat. 


Es werden dann einige Ergebnisse angefiihrt, die mit einem solchen Apparat bei der 
Aufnahme des Verkehrslirms sowie bei der Analyse von verschiedenen Stimmen erzielt 


wurden. 


1. Introduction 


The study of irregular phenomena as annoyances 
of radioelectrical origin, variations of the line vol- 
tage, noise level in a factory, etc. may be subjected 
to statistical analysis, which allows one to sum up 
with a few data a behaviour otherwise difficult to 
describe. This kind of analysis becomes easy if the 
data are obtained by employing one of the various 
methods used for their automatic investigation. 

Considering an acoustical phenomenon (such as 
traffic noise, speech, etc.) different kinds of statisti- 
cal analysis may be employed, each one concerning 
a given aspect of the phenomenon itself. For in- 
stance one can study the statistical distribution of 
the intensity or of the sound level [1], the distri- 
bution of the emitted frequencies, the possible cor- 
relation between frequency and intensity, and so 
on. In these investigations no account is taken of 
the parameter time, that is of the order in which the 
events succeed each other in our phenomenon. 

In this paper we consider the question under 
another point of view: the statistical distribution of 
the duration of given events which occur in the 
phenomenon studied; more exactly, we consider 
“events”, the time intervals during which the sound 
level exceeds a prefixed level. 


A problem showing some analogy with this one is 
found in telephone technique, and concerns the sta- 
tistical distribution of the duration of the telephonic 
communications, of occupation of the telephonic 
exchange, of waiting for the communication. 

By a statistical examination of this kind one can 
establish the average number of events occurring in 
a given time and their probable duration. 

In the case of acoustical phenomena, these sta- 
tistical data may be useful to define a given kind of 
noise, since its annoyance effect is different, de- 
pending whether the limiting level is exceeded sel- 
dom or frequently, for short or long time inter- 
vals, 

This analysis may be carried out graphically, 
though with a heavy outlay of time, when it is pos- 
sible to record graphically the intensity of the acous- 
tical phenomenon being studied. In Fig. 1a we re- 
port the analysis of the sound level of the studied 
noise as a function of time: it has been performed 
at the standard sound level of 80 dB. This level is 
exceeded during various time intervals (in Fig. la 
they are particularly evident); we are concerned 
with the statistical distribution of these intervals, 
that is we measure the duration of each one and 
plot the number of intervals as a function of their 
duration. 
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Fig. 1. Voltages in different stages of an equipment for 
the automatic measurement of the statistical 
distribution of time intervals as a function of 
time; (a) microphone output, (b) output of the 
stage clipping the voltage, (c) output of the 
squaring stage. 


If treated from an analytic point of view, the 
question appears rather complicated. A problem 
showing some analogy with ours has been examined 
by Rice [2], who studies the function of distribu- 
tion of the time intervals between two successive 
zeros in a white noise. This author * Lhe 
problem seems to be quite difficult and apparently 
nobody has as yet given a satisfactory solution. We 
will give some results which are related to the gene- 
ral problem and which give an idea of the form of 
the distribution for the region of small spacings be- 
tween the zeros”. 

The computation of the distribution of the time 
intervals may be made automatically, employing 
electronic or mechanical equipment [3]. 

Every device includes an electronic or mechanical 
stage to clip the microphone output (see Fig. 1 a) 
and reduce it to the form of Fig. 1b. The impulses 
of various forms produced by it are squared in the 
following stage and transformed into impulses of 
constant amplitude and different width (Fig. 1c). 
The various equipments differ for the method by 
which the width of every impulse is computed. 


Says: 


2. Description of apparatus 


The principle on which the equipment used is 
based is the following: when the sound level exceeds 
the given limiting value, it becomes possible for the 
impulses produced by an auxiliary generator, all 
equal and equally spaced in time, to pass through. 
The duration of every event is then evaluated from 
the number of impulses passing in the correspond- 
ing time interval. 
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These impulses are in some way the unit by which 
the intervals are measured: consequently they are 
only approximately evaluated. The scanning time 
has to be established as a function of the pheno- 
menon to be examined: in fact, if the scanning time 
is too long the evaluation of the time intervals is 
too approximate, while employing a very short 
scanning time the evaluation is more exact, but it 
is necessary to have a large number of registers so 
as not to count only the events of a given duration, 
and exclude the longer ones. 

However, a similar choice is necessary in most 
statistical measurements. 

For the automati¢é computation a mechanical 
counter is used, which appears to be more con- 
venient than an electronic equipment, since the 
scanning times must suit acoustical phenomena 


(Fig. 2). 


Fig. 2. Mechanical counter. 


Now we give some details of the equipment which 
has been designed and built. 

The microphone output is connected to a level 
recorder, including an additional contactor; on it 
a moving contact goes solidly with the index draw- 
ing the behaviour of the phenomenon on the paper 
(Fig. 3). The contactor is divided in two parts at 
the level established for the statistical computation: 
when this is exceeded, the moving contact closes the 
circuit of relay P, which actuates the counter: the 
relay is kept excited as long as the index does not 
fall lower than the limiting level. 

In the counter, relay T (fed by a low-frequency 
generator) scans the time according to the wanted 
rhythm (variable from a few seconds to a few 
hundredths of a second). Relay P, driven by the 
moving contact of the level recorder, actuates a 
group of master relays, which give the impulses 
coming from T a path to a chain of relays. When 
P is excited, every impulse actuates successively 
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Fig. 3. Schema of the equipment for the automatic meas- 
urement of the statistical distribution of the time 


intervals. 


one of the relays of the chain. The same relays 
driven by P open a path also to the registers, each 
one connected to one relay of the chain. When this 
is started, the corresponding register is ready to 
start. When relay P opens it raises the blockade of 
the commanding relays, the ones of the chain return 
to their original conditions and only the last register 
reached by the impulses does start (Fig. 12). 

The chain of relays has ten meshes and ten cor- 
responding registers. When a time interval to be 
computed begins, P is excited: the first pulse com- 
ing from the scanning relay actuates the first mesh 
of the chain and prepares the first register to start 
(but it does not actually start). If P continues to 
be excited, the second pulse actuates directly the 
second mesh of the chain: register one is discon- 
nected, the second register is prepared to start. 

So it goes on until the time interval finishes, the 
pen of the level recorder disconnects relay P and 
the last register used does start, indicating the du- 
ration of the corresponding time interval. 

At the end of the analysis the first register indi- 
cates how many events have occurred with duration 
of one period, the second how many of two periods 
and so on. The last register indicates the events 
lasting then periods or more. 


3. Representation of the data 


The data read off the registers can be represented 
in various diagrams, the abscissa of which is the 
duration of the events expressed in seconds. (The 
graduation of this scale is not continuous as the 
chosen method allows of computing events of nine 
different durations, each one being a multiple of the 
scanning time.) 

A first representation of the results reports the 
number of events read on the registers; yet it ought 
to report the number of events in the time unit, so 
that the results of different reliefs can be compared: 
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the following diagrams show the number of events 
per minute. 

Another representation gives the sum of the num- 
bers read on the registers: they mean the number 
of events per minute whose duration is equal or 
shorter than the corresponding abscissa. 

This kind of representation seems more con- 
venient, as the sum contributes to eliminate some 
errors produced in the evaluation of the time inter- 
vals. As the evaluation is not exact every event may 
be computed by one register or by the following 
one. 

Still another representation gives the percent of 
time during which the events have duration indicat- 
ed by the corresponding abscissa. 

An analytical relationship exists among the cur- 
ves of these three diagrams: if y,=/(t) represents 
the second curve, the first is given by y, = df(t) /dt, 
the third by y;=t: df(t) /dt. 


4. Applications 


The equipment which has been described has been 
used for a number of investigations to show its pos- 
sibilities. 

First of all a traffic noise has been studied, in an 
animated point of the town. The time during which 
the survey has been made is about half an hour, the 
scanning time ranges from a few tenths of a second 
to one second: we consider the time intervals dur- 
ing which the sound level exceeds 80 dB. 


Duration ——— 


Fig. 4. Statistical analysis of a traffic noise: number 
of events per minute as a function of their du- 
ration. 


Fig. 4 reports the number of events per minute, 
whose duration is indicated by the corresponding 
abscissa. 

The diagram of Fig. 5 has been deduced employ- 
ing a good magnetic record of the phenomenon: 
then its analysis can be repeated many times, with 
different scanning times. Fig. 5 shows the number of 
events per minute whose duration is equal or shor- 
ter than the one indicated by the abscissa. The data 
obtained for different scanning times (represented 
by different symbols) agree in a curve with asintho 


oo2 
0 2 4 6 Si asa. aK) 
Duration ——— 
Fig.5. Statistical analysis of the traffic noise for dif- 


ferent scanning times: number of events per 
minute whose duration is equal to or shorter 
than the one indicated by the abscissa 
(@@e@els, xxx 0.58, 000 0.33 s). 


tic behaviour: this means that in the studied pheno- 
menon there are no events whose duration is longer 
than a given value, indicated by the abscissa. 


Percent of time —~ 
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Duration ——= 


Fig. 6. Percent of time during which in the traffic noise 
occur events of a given duration, indicated by 
the abscissa. 


Fig. 6 represents the percent of time during which 
the events have the duration indicated by the abs- 
cissa: the curve shows a maximum, corresponding 
to the more probable duration of the events, about 
two seconds. 

Fig. 7 gives, in the same form of Fig. 5, the 
results of an investigation on the traffic noise: for 
the computation a scanning time constant of 0.5s 
has been used; the computation has been repeated 
for three different sound levels. 

Other results of statistical computations have 
been deduced by the investigation of different voi- 
ces, reading prose. The level for which the statistical 
distribution is sought varies from 60 to 70 dB, the 
scanning time from 0.5 to 0.066; every investiga- 
tion corresponds to 2 or 3 minutes of speech. 

In a similar way Fig. 8 gives the results of an 
analysis at constant sound level (65 dB) and con- 
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Fig. 7. Statistical analysis of the traffic noise: number 
of events per minute whose duration is equal 
to or shorter than the one indicated by the abs- 
cissa for different sound levels. 
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Fig. 8. Statistical analysis of three different voices for 
both sound level and scanning time constant: 
number of events per minute whose duration is 
equal to or shorter than the one indicated by 
the abscissa. 


stant scanning time for three different voices, all 
reading the same prose: while male voice B exceeds 
65 dB only for time intervals shorter than 1s, 
female voice A and male voice C have a quite dif- 
ferent behaviour and include also events whose 
duration is about 2 and 3s. 

On a prose read by a male voice the statistical 
analysis has been repeated for different levels; the 
scanning time has been kept constant at 0.2s. At 
high level (70 dB) the events are very short, with a 
maximum duration of 0.5s; the diagrams at lower 
levels (65 and 60 dB) do not show any saturation, 
even for time intervals of 2s (Fig. 9). 

Still using magnetic records we have carried on 
the statistical analysis of the time intervals at con- 
stant sound level and scanning time for a prose read 
by a female voice: before the computation the 
microphone output passes through a set of one- 
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Fig. 9. Statistical analysis of a male voice for diffe- 

| rent sound levels: number of events per minute 
whose duration is equal to or shorter than the 
one indicated by the abscissa. 


third-octave band-pass filters: a set of diagrams is 
deduced by this measurement, each one correspond- 
ing to a different band of speech frequency. The 
diagrams 


various deduced are represented in 


Fig. 10. 


0 04 08 . 12 AS Ti 28) 
Duration ——— 


Fig. 10. Statistical analysis of a filtered female voice 
passed through one-third octave band-pass fil- 
ters: number of events per minute whose du- 
ration is equal to or shorter than the one indi- 
cated by the abscissa. 


This diagram can convey an idea about the fre- 
quency distribution of the voice examined: in 
effect the frequency bands in which the components 
are weaker have but short peaks exceeding the limit- 

ing levels. The curves drawn at 160, 800, 1000 and 

1250 c/s show this behaviour. In other frequency 
bands, mostly about 500 and 600 c/s, a number of 
events occur, with various duration, up to more 
than 2s, 


a 
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Fig. 11. Statistical analysis of a female voice: number 
of events per minute whose duration is less 
than 0.5s as a function of the frequency. 


From the set of curves of Fig. 10, of very diffe- 
rent behaviour, one can deduce the diagram of 
Fig. 11: it shows the number of events per minute 
whose duration is less than 0.5s as a function of 
the frequency of the filter. 

A large number of other investigations and ap- 
plications may be made with the equipment describ- 
ed: which indeed is adapted to a great variety of 
phenomena and particular cases. 


Acknowledgement 


The author wishes to acknowledge the assistance 
received from ing. A. Caccta, of the S.T.I.P.E.L., 
and from ing. G. Miterto, of the $.T.E.T., building 
the equipment. 

Research executed with the financiary help of the 
High Authority of the C.E.C.A. (Communité Euro- 


péenne Charbon et Acier). 
(Received December 15th, 1956.) 


Appendix 


Fig. 12 reports the schema of the mechanical 
counter: it comprehends a number of relays (indi- 
cated by capital letters, the corresponding contacts 
of which are indicated with small letters) and ten 
registers Z. 

Before counting any time interval the different 
contacts have the position indicated by the schema: 
only relay D is excited, reaching the ground by con- 
tacts x and p. 

As the level recorder closes by its pen the circuit 
of relay P, p changes its position, D is disconnected 
from the ground and, as the excitation leaves it, 
contact d, is closed. This contact connects to ground 
(by contact p) relay I, whose contact i, connects to 
ground relay V, which closes contact v, . 

In these conditions the registers and the relays 
of the chain begin counting the pulses which are 
considered as time units in this measurement, scan- 
ned by T and an auxiliary generator. t periodically 
connects to ground, by contact d, and p, the chain 
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Fig. 12. Schema of the mechanical counter. 


of relays A. A, is excited by the first pulse: in 
starting it closes contact a, , after this also B, starts 
and changes the position of contact b,. A, closes 
also the contact a, of register Z,, which still does 
not start because its circuit is interrupted from 
ground by contact i, . 
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The second pulse, by the new position of b,, 
reaches directly A, , which inserts B, in the chain of 
operations, repeating what has happened in the first 
mesh. In the circuit of register Z,, contact ay dis- 
connects it, while another contact ay inserts register 
Zp. : 

The operations continue in other meshes until 
relay P is excited: if the time interval to be com- 
puted comes to the end by then, p changes its posi- 
tion, I is disconnected from ground and its contact 
i, connects to ground (by contacts vy, a: and ag) 
register Z, which this time does start. 

At the same time the circuit of relay V is inter- 
rupted: only Dis again connected to ground, and 
everything is in condition to start again for the com- 
putation of a new time interval. 
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THE ABSORPTION OF SOUND 
IN ARGON, NITROGEN AND OXYGEN 


by W. Temesst and H. D. Parsroox 


Acoustics Laboratory, University of Liverpool 


Summary 


Some improved values for sound absorption in argon, nitrogen and oxygen in the fre- 
quency/pressure range 0.45 to 12.6 Mc/s atm are presented and compared with theoretical 


values. 


Sommaire 


On donne et compare avec les valeurs théoriques quelques résultats améliorés de mesure 
de l’absorption du son dans l’argon, |’azote et l’oxygéne dans la gamme fréquence/pression 


allant de 0,45 a 12,6 MHz/atm. 


Zusammenfassung 


Es werden einige verbesserte Werte fiir die Schallabsorption in Argon, Stickstoff und 
Sauerstoff im Bereich Frequenz/Druck yon 0,45 bis 12,6 MHz/atm mitgeteilt und mit 


theoretischen Werten verglichen. 


1. Introduction 


The Stokes-Kirchhoff equation [1] for the absorp- 
tion of sound in a gas, and its modification by the 
“relaxation” theory [2] to correlate with the “ano- 
malous” absorption found in certain gases, have 


figured large in the relevant literature. The theory 
is, however, not entirely satisfactory and the work 
now reported was concerned with the validity of the 
Stokes-Kirchhoff equation, and the anomalous ab- 
sorption reported in some monatomic and diatomic 
gases. ; 
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The Stokes-Kirchhoff equation is generally ex- 
pressed in the form 


2 x? f2 y — 
a2 (424? 2) 


i 
oc 


“Gee ae 


is the amplitude absorption coefficient/em, 
f the frequency (c/s), 

o the ambient density (gm/cc), 

e the sound velocity (cm/s), 
q 
y 


where @ 


the viscosity (poise), 
the ratio of the specific heats (C,/C,), 
the specific heats at constant pressure and 
constant volume (cal/gm °C) 
and x the thermal conductivity (cal/ems °C). 
The dissipative losses are assumed to be entirely 
due to viscous and heat conduction processes; losses 
from other processes [3] are expected theoretically 
to be small. The above form of the equation is ob- 
tained by assuming that the sound pressure is al- 


Caan 


ways in phase with the condensation, or, alterna- 
tively, that a volume change can occur without loss 
of energy. Such an assumption leads to an expres- 
sion relating the second or dilatational viscosity 
and the shear viscosity 7/ of the fluid: 27 +3” =0. 
It has been suggested [4] that a better agreement 
between experimental results and theory would exist 
if the dilatational viscosity «’’ could be set equal 
to zero. (The two viscosity coefficients, 7 and ww”, 
are related to the bulk viscosity coefficient 7 by 
2y+3 mu =3y7 [17].) 

From the Stokes-Kirchhoff equation it can be 
shown that the sound absorption in ‘an ideal gas will 
be proportional to /?/p and independent of tempe- 
rature. The effect of a variation of temperature in 
a non-ideal gas is complex since nearly all the para- 


meters in the equation are to some extent tempera-— 


ture-dependent; approximate calculations, using pu- 
blished values for the parameters, show that the 
absorption rises with increasing temperature for 
nitrogen, oxygen and argon. 

Relaxation theories [2] have been developed to 
explain for certain gases the observed absorptions 
in excess of those predicted by the Stokes-Kirchhoff 
equation. This theory is essentially phenomenologi- 
cal and it has proved difficult to assess its signifi- 
cance. In monatomic gases the transfer of energy 
in the form of acoustic waves takes place by col- 
lision processes in which the translational energy 
changes of atoms are shared; at the energy levels 
involved no other form of energy would seem pos- 
sible. However, in a diatomic gas, several new de- 
grees of freedom appear and the energy of the mole- 
cules may be distributed between. the translational, 
vibrational and rotational modes. The interchange 
of energy between these modes can, at certain fre- 
quencies, exercise serious effect on sound absorp- 
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tion, and it is to relaxation phenomena, that is, to 
the time-dependent interchange of energy between 
the translational and other possibly dissipative mo- 
des, that the anomalously large absorptions observed 
in some polyatomic gases have been attributed. 

An examination of the published results up to 
1951 [18] for absorption in gases shows consider- 
able scatter in results at ultrasonic frequencies, and 
it seems unlikely, even in the best results, that an 
accuracy greater than + 10% can be claimed. This 
would seem to be inadequate for correlation of ex- 
perimental results with theories especially remember- 
ing that the changes in value sought can be less than 
the above error. 

The majority of absorption measurements in gases 
at the higher frequencies has been made using the 
ultrasonic interferometer. Pumper [5] gives a criti- 
cal account of the work before 1938; his main con- 
clusion is that there are considerable errors inherent 
in the experimental methods used and that these 
errors account for the discrepancies between results 
obtained by different workers. Pumper also describes 
experiments which demonstrate that the interfero- 
meter can be made to yield differing results by 
slight modifications of technique, modifications which 
would, in theory, not affect the operation. A con- 
siderable number of interferometric measurements 
have been made since 1938, but there is little evi- 
dence that Pumprr’s criticisms do not apply equally 
to these. Krasnoosuxin [16] has developed a theory 
of the interferometer taking into account the non- 
uniformity of the acoustic field. This theory predicts 
effects similar to those observed by Pumper. KELLER 
[6] in 1940 used an optical (schlieren) method; 
the scatter of the results obtained would appear to 
indicate that the method is not of a very high pre- 
cision. Parker, Apams and Sravseru [7] (1953) 
have described a “direct”? measurement technique 
in which the ultrasonic absorption is determined by 
measuring the sound intensity in a gas as a function 
of the distance from the sound source. Rochelle salt 
crystals were used as sound source and receiver at 
frequencies around 66 kc/s. Greenspan [8] has de- 
scribed a two crystal “interferometer” which is in 
effect a direct transmission method suitable for con- 
ditions of high absorption. So far as is known, this 
method has been used only in monatomic gases. 
These two latter approaches seem to yield results of 
a much higher consistency than those obtained pre- 
viously. 

In the work reported here, an improved method 
of measurement has been developed in the frequency 
range 400 kc/s to 3 Me/s, and it is thought that, 
from the results obtained in nitrogen, oxygen, argon 
and dry air, more adequate discussion of theories 
may ensue. 
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2. Apparatus 


The apparatus consists essentially of two trans- 
duceres, a source and a receiver, mounted facing 


each other in a thermostatically controlled gas-tight 


enclosure. The receiving transducer is mounted on 
a precision thread operated through a sealed bear- 
ing, permitting the path length between the two 
transducers to be varied from 0.2 to 10cm. The 
electrical input to the source transducer is pulsed 
(a) to prevent standing waves between the two 
transducers, (b) to permit separation of wanted sig- 
nals from those travelling other than through the 
gas, and (c) to avoid overheating the transducers 
and the gas. Quartz and barium titanate discs ope- 
rating in assumed thickness mode vibration (piston- 
like fashion) have been used as transducers. 

A typical transducer for both source and receiver 
is a quartz disc 25 mm diameter and 2.5 mm thick. 
This would have an expected fundamental thickness 
mode of frequency 1 Mc/s. The reverse side of the 
transducer is cemented to a block of polythene which 
acts mainly as a resistive loading reducing the un- 
damped transducer “Q” from greater than 10* to 
approximately 15 [9]. The transducer is some 75 
wavelengths in diameter and when vibrating in a 
near piston-like manner produces an almost plane 
wave sound field in the gas. Reference will be made 
later to the elimination of errors due to the devia- 
tions from a plane wave field. 

A modulated R.F. amplifier producing 500 rect- 
angular pulses per second of adjustable duration 
(6 to 60 ts) is driven by a simple R. F. oscillator. 
With the damped transducers a high stability oscil- 
lator was not found necessary. A tuned air-cored 
transformer with variable turns ratio is used to 
match the amplifier to a source transducer. 

The receiver-transducer output is amplified in a 
stable wideband amplifier (80 dB maximum gain), 
filtered to reduce hum and noise, and displayed on 
a D.C. oscilloscope. A Y shift control with a moni- 
toring “substandard” meter connected to the input 
of the D.C. amplifier permits measurement of the 
voltage required to move the waveform by an 
amount equal to the peak to peak height. The me- 
thod is essentially a null method, since it only as- 
sumes that the D. C. displacement voltage is directly 
proportional to the amplitude of the R. F. pulse — an 
assumption justified experimentally. A further as- 
sumption, which would appear to be justified, is that 
there is a linear relation between the received sound 
signal and the transducer output. 

Gases are admitted to the working enclosure (vo- 
lume about 1.3 litres) through a cold trap contain- 
ing either liquid air or solid carbon dioxide in ace- 
tone. The gas pressure is read from an open-ended 
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mercury manometer. A two stage rotary vacuum 
pump with an ultimate pressure of not more than 
0.1 mm was used to pump out the system and, since 
the lowest pressure used was about 5 cm, this re- 
presents a maximum impurity content of about 0.2%. 
Repeated “flushings” of the system before use with 
the gas to be examined considerably reduce this 
impurity figure. Care was taken in the construction 
of the system to eliminate materials likely to show 
serious outgassing. 

Supplies of argon obtained from “British Oxygen 
Gases Ltd.” were of 99.97% stated purity. Oxygen 
and nitrogen from the same source were of com- 
mercial purity, presumably about 99%. 

To provide some indication of the possible effects 
of impurities, oxygen and nitrogen were prepared 
chemically and the properties of the two different 
supplies compared in each case. Oxygen was pre- 
pared by the action of distilled water on “A. R.” 
quality sodium peroxide and dried by passing 
it through phosphorous pentoxide and a trap im- 
mersed in solid carbon dioxide in acetone. The 
reaction involved in the preparation is a compara- 
tively simple one in the sense that there are no 
secondary products and it would appear that gas of 
a very high purity should result. Nitrogen was pre- 
pared by the warming of a solution of ammonium 
chloride and sodium nitrate and passed through a 
solution of potassium dichromate acidified with ace- 
tic acid. This latter solution is claimed [10] to re- 
move any of the higher oxides of nitrogen which 
may be formed! Phosphorus pentoxide and solid 
carbon dioxide in acetone traps were added to re- 
move water. The reaction involved in the prepara- 
tion of nitrogen is a complex one and produces 
traces of ammonia and of various oxides of nitro- 
gen. The only impurity which might not be removed 
by the purification system is a trace of nitrous 


oxide. Further work is in hand to examine more 


closely any impurity effects. 

In operating the gas system the following proce- 
dure was employed. The enclosure was pumped out 
for about half an hour and the gas then admitted 
slowly through the cold-trap, up to a pressure of 
about one atmosphere. This filling served to “flush” 
the system. The enclosure was then re-pumped for 
about five minutes to remove the “flushing” gas and 
filled to the maximum required pressure, which in 
any case was limited to 80 cm mercury above atmo- 
spheric pressure by the construction of the enclosure. 
After an absorption measurement at the highest 
pressure, some gas was pumped out and the next 
measurement taken and similarly down to the lowest 
workable pressure. This lower limiting pressure is 
the pressure at which the output signal becomes 
comparable to the amplifier noise level. 
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3. Experimental procedure 


Initially interest centred on the method as a 
means of measurement of absorption in gases with 
particular reference to the errors involved. For this 
purpose argon was chosen because (a) it was easily 
obtainable in quantity and high purity, (b) previous 
absorption figures for this gas seemed more consis- 
tent than for some others and (c) the oc value for 
the gas is convenient for experimental work. Both 

quartz and barium titanate transducers were used. 

With 25 mm diameter quartz transducers at ap- 
proximately 932 kc/s the received pulse amplitude 
was measured at various distances from the source. 
Fig. 1, showing a plot of the natural log. of the 
amplitude against the distance of the receiving trans- 
ducer from the source, is a typical result for argon. 


Ln (amplitude) ——— 


“0 16 32 48 mm 64 
Transducer separation (arbitrary zero) ——>— 


Fig. 1. Received signal as function source — receiver se- 
paration (argon, 932 kc/s, 0.97 atm). 


The Amplitude A of the sound field at x from the 

source can be represented by 
A=A, exp(—ax) exp(—fx), 

where Ay represents the amplitude at the source, 
a the amplitude absorption coefficient in nepers/em 
and f the apparent absorption due to the departure 
of the beam of sound from a plane wavefront (i. e. 
to beam divergence). A series of measurements 
made in argon is shown in Table I and the graph 
of a+ plotted against 1/p produces the straight 
line passing through the origin shown in Fig. 2. 


Fig.2. Absorption in argon as a function of pressure 
(quartz transducers, 932 ke/s). 
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This is the result to be expected if the transmission 
is essentially plane wave (i. e. if 6=0) and the true 
absorption is inversely proportional to the pressure. 
The slope of this graph gives ap and _ hence 
a p/f? = 1.89 x 10-18, where p is measured in atmo- 
spheres (one atmosphere is defined as equivalent to 
76 cm of mercury at 0° C), and f in ¢/s. 

Having established that, for the particular trans- 
mitter used, the beam divergence effects are negli- 
gible, each single measurement can be used to work 
out ap/f? and the series of values so obtained used 
to derive the standard deviation of the mean. In this 
case the mean value of a p/f? is 1.89 x 107", with 
a standard deviation of 0.01 x 10° 1°. 

A further series of measurements were made 
using barium titanate transducers of 32 mm dia- 
meter and approximately 6 mm thickness, resonant 
at 405 ke/s. The results are shown in Table II and 
plotted as apparent absorption:against 1/p in Fig. 3. 
This graph intercepts the absorption axis at 0.103 
and this is treated as the apparent absorption due 
to the divergence of the ultrasonic beam, It has 
been shown [14] that, with a transducer of the di- 
mensions used, the expected loss of amplitude due 
to divergence of the beam is considerably less than 
that apparent here. The result obtained is probably 
due to the barium titanate vibrating in a complex 
mode, and hence at different amplitudes at different 
points on its surface. The slope of Fig. 3 yields 
a p/f?=1.91 x 10718, which agrees with the result 
at 932 ke/s assuming the Stokes-Kirchhoff equation 
applies. 


nepers/cm 


0.3 
| 


lo —~— 


Fig.3. Absorption in argon as a function of pressure 
(barium titanate transducers, 405 kc/s). 


4. Experimental accuracy 


The mercury manometer used for the determina- 
tion of pressures could be read to 0.1 cm or better 
and since pressure readings are derived from the 
difference of two observations, an overall accuracy 
of not less than 0.2 cm could be obtained. Thus all 
pressures higher than 20cm of mercury have a 
measurement accuracy of better than 1%. 
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The temperature of the waterbath was maintained 
to +0.2° C by means of a “Sunvic’’ thermostat and 
a circulating pump; this system was operative at 
temperatures above room temperature, 20°C, At 
lower temperatures a bath of melting ice provided 
ie Comte Bs 

Frequencies were measured on a frequency meter 
with crystal calibration, the crystal being compared 
with the B. B. C. 200 ke/s transmission. The overall 
accuracy is better than 0.1%. 

A wide-band R.C. coupled aperiodic amplifier 
was used for signal amplification, an amplifier of 
this type with heavy negative feedback having ex- 
tremely low non-linear distortion at the signal levels 
used (about 2 V peak to peak output). The actual 
figures for this distortion were not measured but can 
be estimated at 0.2 to 0.4%. 

It was essential with the method used that the 
gain of the amplifier and the amplitude of the trans- 
mitted pulse should remain absolutely stable during 
a run (approximately 10 minutes). To achieve this, 
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all important power supplies were electronically 
stabilised, with the exception of the heater supplies 
to the values which were derived from large capacity 
lead-acid accumulators. The stability was such that 
when the equipment was warmed up the output 
signal did not change detectably over a period of 
10 minutes. 

It seems unlikely that the total systematic error 
will be appreciably greater than 1.5%. 


5. Results 


In the following results a is the amplitude ab- 
sorption per_cm, p is the pressure in atmospheres 
and / is the frequency in c/s. 


a) Argon 

Tables I, II and III give the absorption values for 
argon at 932 ke/s, 937.5 ke/s, 405.0 ke/s, 1.172 
Mc/s and 3 Mc/s at 25° C. Values are corrected as 
necessary for any divergent beam effect. 


Table I 


Absorption values in argon 


Argon: 25°C; 932 and 937.5 ke/s 


| | 
Frequency (ke/s) | 932, 932 932 932 932 932 | 932 | 932 
Pressure (atm) 1.28 1.28 1.22 1.22 1.06 1.06 0.987 | 0.987 
a + B (nepers/em) 0.130 0.132 0.134 0.135 0.153 0.149 0.163 0.167 
(« + B) pif? (x 10) 1.92 1.94 1.89 1.90 1.87 1.82 1.85 1.90 
Frequency (ke/s) | 932 932 | 932 932 | 937.5 937.5 |. 937.5 937.5 
Pressure (atm) | 0.987 0.980 | 0.980 0.965 | 0.551 0.551 | 0.500 0.500 
a + £ (nepers/cm) 0.168 0.169 0.168 0.169 | ‘0.808 °| ' 0.808_) 0:336 0.339 
(« + B) p/f® (x 1013) 1.89 1.91 1.89 190 1.92 1.99) 31) 4.68 1.85 
Mean value of («+ 6) p/f? = 1.887 x 107%. 
Standard deviation of mean = 0.01 « 10-®. 
Table I 
Absorption values in argon 
Argon: 25°C; 405 ke/s 
Pressure (atm) | 1.01 0.970 0.724 0.562 0.376 0.260 0.232 0.176 0.131 
a + f (nepers/em) 0.134 0.136 0.143 0.162 0.198 0.252 0.244 0.284 0.340 
Table II 


Absorption values in argon 


Argon: 25°C; 1.172 Me/s 


Pressure (atm) | 2.08 1.98 192°) 181 1.77 1.20 | 1.05 0.781 0.664 | 0.505 
« (nepers/cm) 0.120 | 0.122 | 0.1389 | 0.133 | 0.150 | 0.231 | 0.247 | 0.325 | 0.391 | 0481 
ap/f? (x 1038) 1.82 1.76 1.94 1.76 1.93 2.02 1.88 1.85 | 1.89 1.84 
Mean value of «p/f? = 1.869 x 10-8, 
Standard deviation of mean = 0.03 x 1073. 
Argon: 25°C; 3.000 Me/s 
Pressure (atm) | 2.00 1.61 1.36 1.14 0.978 0.830 0.539 | 0.416 0.326 | 0.238 
a (nepers/em) 0.84 1.04 1.24 1.46 1.73 SOT ames to 4.05 5.17 |e 
x p/f2( x 10") 187 (01.862 45 1.87 1 1.8B eae Vad 187 | 187 | 87 9) ena 
Mean value of « p/f? =) 1-863' x. 10-3; 
Standard deviation of mean = 0.01 x 107%, 
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b) Oxygen 

Absorption in oxygen was measured at 405.0 ke/s, 
932 ke/s, 1.172 Me/s and at 3.000 Mc/s. In all cases 
the absorption was found to be proportional to f2/p. 
The results are set out in Table IV. 

To compare the commercial gas with that pre- 
pared chemically, absorption values were measured 
at a series of pressures in the two samples of gas. 
The results (Table V) appear to indicate that there 

is no difference between the two samples. 


c) Nitrogen 

Table VI sets out the absorption values in nitro- 
gen at frequencies of 704.2 kc/s, 932 ke/s, 1.172 
Mc/s and 3.000 Mc/s. Over the range of pressures 
and frequencies involved nitrogen displays a varia- 
tion of absorption with frequency and pressure si- 
milar to that found in oxygen and argon. Nitrogen 
prepared chemically was examined at 943 kc/s and 
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gave ap/f?=1.68 x 10~! which is in close agree 
ment with the value of 1.72 x 10713 obtained at a 
similar frequency using the commercial gas. 


60 °C 


i ay (ee EE, 
0 20 40 


Temperature — 


Fig.4. Temperature dependence of absorption in oxy- 
gen and nitrogen. 


Table IV 


Absorption values in oxygen 


Oxygen: 25°C; 405 ke/s 


Pressure (atm) 0.804 | 0.772 | 0.65 | 0.512 | 0.387 | 0.310 | 0.260 | 0.181 | 0.134 | 0.093 
a (nepers/em) 0.043 | 0.041 0.053 | 0.064 | 0.085 ODTE W201 29s O75 0.245165 02365, 
a p/f? (x 101%) | 2.26 | 1.93 | 2.10 | 2.00 | 2.00 2.10 | 2.06 1.94 2.00 |. 2.08 
| | ! i 
Mean value of «p/f? = 2.05 x 10-8 
Standard deviation of mean = 0.03 « 10-4 
Oxygen: 21°C; 931 ke/s 
Pressure (atm) | 207 | 1.845 | 1.64 | 1.24 | 1.01 | 0.807 | 0.624 | 0.482 | 0.308 | 0.228 
 @ (nepers/em) 0.084 | 0.092 0.109 0.140 | 0.172 | 0.217 | 0.272 0.360 0.556 0.740 
a p/f2(< 108) 2.00 1.96 eh 2. 0Vere| E90. [2.028 7.05 v) 1, 99-° ers, “los 
Mean value of «p/f? == 1209 105° 
Standard deviation of mean = 0.01 x 10-°¥ 
Oxygen: 25°C; 1.172 Mc/s 
Pressure (atm) 1.684 | 1.186 1.04 / 0.835 ) 0.795 | 0.695 | 0.595 | 0:504 | 0.408 | 0.351 
% (nepers/em) 0.162 0.214 | 0.262 | 0.315 | 0.329 | 0.407 | 0.452 | 0.53 | 0.64 0.74 
op/f? (x 10!) 1.99 2.01 1.98 1.91 Ol Sy |e 196) 2) 1.95 1.90 1.89 
Mean value of « p/f* = "1:96 x 10-2 
Standard deviation of mean = 0.02 « 107. 
Oxygen: 25°C; 3.000 Me/s 
Pressure (atm) | 2.01 | 1.65 1.37 1.504 | 0.844 | 0.592 | 0.469 | 0.426 | 0.326 
a (nepers/cm) | 0.92 | Lil 1.36 1.76 2.19 3.25 | 3.97 4.40 5.81 
ap/f2( x 1038) | 2.05 | 2.04 | 2.03 | 2.06 | 204 | 213 | 209 | 2.07 | 2.10 
Mean value of «p/f* = ORD LOT, 
Standard deviation of mean = 0.01 x 107%. 
Table V 
Absorption in “cylinder” and “chemically prepared” oxygen 
Oxygen: 25°C; 943 ke/s 
“Cylinder” “Chemically prepared” 
Pressure (atm) | 1.11 | 0.974 | 0.768 | 0.491 | 0.263 | 1.11 | 0.975 | 0.848 | 0.490 | 0.266 
a (nepers/cm) | 0.159 0.183 0.236 0.371 | 0.676 0.16 0.183 0.209 0.366 0.637 
a p/2( x 10") 198 2.01 | 2.04 | 202 197 199 | 201 1.97 1.99 | 1.88 
Mean value of «p/f? for “Cylinder”? oxygen = 2.00 x10. 


Mean value of «p/f* for “Chemically prepared” oxygen = 1.97 x 107, 


an 
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Table VI 
Absorption values in nitrogen 
Nitrogen: 25°C; 704.2 ke/s 


Pressure (atm) | 1.90 | 1.66 | 1387 | 115 | 0.952 | 0.789 | 0.71 ; ; 
o. (nepers/em) | 0.047 | 0.052 | 0.065 0.074 0.089 | 0.108 | 0.12 0.196 | 0.262 
ap/f2( x 1018) J 1a2 | 174 | 7s | Lm | 1.70 | Lm | 176 1.74 | 1.77 


Mean value of a p/f* == 17455 L038 
Standard deviation of mean = 0.01 x 107. 


Nitrogen: 23°C; 932 ke/s 


Pressure (atm) 73 | 1.51 | 1.25 0.545 | 0.478 
x (nepers/cm) 0.073 | 0. 0.087 | 0.097 | 0.12 0.274 | 0.312 
op/f2( x 1038) U7.) LTO 2) 178s) 1862) LI 1.7L Saal 
Mean value of «p/f? == O02 pO 
Standard deviation of mean = 0.01 x 107 


Nitrogen: 25°C; 1172 ke/s 


Pressure (atm) 20 Fa 
% (nepers/em) | 0.114 0.124 
ap/f?( x 1038) | 1.74 


Mean value of «p/f? 
Standard deviation of mean 


Nitrogen: 25°C; 3.000 Mc/s 


Pressure (atm) 2.00 | 157 | 1.34 | 1.08 | 0.844 | 0.677 | 0.527 | 0.395 | 0.332 

x (nepers/em) 0.78 | 0.99 | 117 | 144 | 1.89 | 2.39 | 2.99 | 414 | 4.70 | 5.61 

ap/f2( x 1018) | 273 | 7B | 174) 176 | 77 | 179° | 1.74) | 1-8) ee 
Mean value of « p/f? == G LO: 


Standard deviation of mean = 0.01 x 10-8. 


Table VII 


Variation of absorption with temperature in oxygen and nitrogen 


Oxygen | Temperature (°C) 0 | 20 | 25 50 60 

920—950 ke/s |  ap/f2(x 1033) 1.94 | 2.00 | 2.105 

3.000 Me/s | ap/f®( x 1018) 1.95 | 2.06 2.135 
| | 

Nitrogen | Temperature (°C) | 0 20 25 50 60 

920—950 ke/s |  ap/f2(x 1018) | 1.665 1.69 | | 1.74 

3.000 Mc/s ap/f2( x 1013) 1.69 | 1.755 | 1.81 

Table VIII 


Absorption values in dry air 
Air: 25°C; 935 ke/s 


Pressure (atm) | 0.893 

« (nepers/cm) | 0.175 

a p/f? (x 1018) | 1.79 

Pressure (atm) 0.382 

« (nepers/cm) 0.416 

a p/f? (x 1018) 1.82 / 
Mean value of «p/f? = 1.856 x 10-3, 
Standard deviation of mean = 0.01 x 10-4, 

d) Dry Air and oxygen at 3.000 Me/s and at 920—950 ke/s. 


The absorption in dry air at 935kc/s is set out In this latter case the barium titanate transducer 
in Table VII. (The results presented above comprise used had a large temperature coefficient of fre- 
about 25% of the total number taken. The remainder quency which made working at a constant frequency 
are in good agreement and serve merely to demon- impossible. The a p/f? values are set out in Table 
strate the repeatability of the measurements.) VIII and plotted in Fig. 4. The results may be ex- 
pressed approximately in the form that in oxygen 
the absorption is proportional to the 0.48 power of 

Some preliminary measurements were made of the absolute temperature and in nitrogen to the 0.31 
the temperature coefficient of absorption in nitrogen _ power. 


Temperature effects 


3 


: 
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6. Discussion 


For argon the Stokes-Kirchhoff absorption is given 

by Marxuam, Bryer and Linpsay [18] as 

a p/f? =1.88 x 10~** (or 1.85 x 10713) 

and the mean of the measured values corresponds 
closely to this (1.89 x 10-4). Using currently avail- 
able thermodynamics constants one obtains a value 
of 1.84 x 1071 at 25° C. Since similar figures were 
obtained at 66 kc/s by Parker et al. [7] and at 
11 Me/s by Greenspan [8], it would seem that the 
Stokes-Kirchhoff equation accounts completely for 
the measured results at frequencies between 66 kc/s 
and 11 Mc/s. The results of Kerirr [6] at 4.25 
Mc/s, although not very precise, are also in agree- 
ment. 

Table VII shows that the absorption in air is in- 
versely proportional to pressure and corresponds to 
an ap/f? value of 1.856 x 10743. This result is of 
no serious theoretical significance except that it 
shows that the absorption in such a mixture is ap- 
proximately that to be expected from the absorption 
of its components. It also shows that impurities such 
as small quantities of carbon dioxide produce no 
large effects. 

In oxygen and nitrogen the mean measured values 
are 2.02 x 10 * and 1.74 x 107" respectively. The 
Stokes-Kirchhoff values are not precisely available. 
Table IX sets out thermodynamic constants publish- 
ed in some current works of reference and demon- 
strates that errors of the order of at least 3 — 4% are 
possible in the calculated value. In spite, however, 
of these variations the difference between the meas- 
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ured and calculated values seems to be significant 
and again this view is supported by the measure- 
ments at lower frequencies by Parker et al. [7]. 
Expressed in the form dxp/as, we find for nitro- 
gen 1.34 and for oxygen 1.27 and it is necessary 
to consider why these ratios depart from unity. Four 
explanations appear worthy of consideration: 


a) Experimental errors. The theory of the present 
method is relatively simple and does not offer any 
other clear interpretation of the observations. 


b) There is the possibility that a serious error is 
present in one of the thermodynamic constants in 
the Stokes-Kirchhoff equation. From the experimen- 
tal point of view the most difficult constant to 
measure accurately is the thermal conductivity. How- 
ever, to bring the calculated absorption up to the 
observed value requires a doubling of the thermal 
conductivity and an error of this magnitude seems 
extremely unlikely. An estimate of the thermal con- 
ductivity can be obtained from the other constants 
of the gas using the formula x= (9 y—5)nC, [21] 
and this yields values for x in fair agreement with 
those obtained experimentally. 


c) A third possibility is that most commonly put 
forward to explain the phenomena, viz. relaxation. 
Herzretp and Rice [2] suggested that a “slowness 
of energy exchange between the translational and 
the internal degrees of freedom of a gas molecule” 
could affect the absorption of sound waves. In both 
nitrogen and oxygen rotational energy levels are 
possible, and, subject to certain fairly drastic as- 


Table IX 


Numerical constants from some current works of reference used in the calculation of sound absorption, corrected 
where necessary to 25° C. 


Nitrogen | Oxygen 
ic pea Maximum ) 2 Minimum | Maximum | Minimum 

Velocity 352.7 x 10° 352 x 10? / 351 x 10? | 350 x 10? 
(cm/s) B | I B | I 
Density 0.1146 x 10? | 0.1308 « 10-? 
(gm/ec) | H | H | 
Viscosity 1775 x 1077 1762 x 10-7 2060 x 10-7 . 2050 x 10-7 
(poise) I H | B . Hi 
Thermal | 60.6 x 10-® 65.1 x 10* 62.7 « 10-* . 59.2 « 10-% 
conductivity I K KL H 
(cal/em s °C) | 
Specific Heat . 0.249 0.2477 / 0.218 
(cal/gm °C) I B Bi 
Ratio of the 1.404 | 1.401 
Specific Heats B B 

B: Beranex, L. L., Acoustic measurements. Wiley & Sons, New York 1949. 

Te International critical tables. McGraw-Hill Book Co., New York 1933. 

H: Handbook of chemistry and physics. Chemical Rubber Publishing Co., Ohio 1955. 

K: Kennarp, E. H., Kinetic theory of gases. McGraw-Hill Book Co., New York 1938, 

KL: Kaye, G. W.C. and Lasy, T. H., Physical and chemical constants. Longmans, London 1956. 


These figures produce the following values for a p/f*  10**: 


Max. 1.322, 
Max. 1.634, 


Nitrogen: 
' Oxygen: 


Min, 1.275; 
Min. 1.590. 
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sumptions, it is possible to calculate relaxation times 
for the transfer of energy between translational and 
rotational states. Parker et al. [7], whose data are 
substantially in agreement with that presented here, 
work out that for both oxygen and nitrogen the re- 
laxation time of the process involved is approxima- 
tely 5x 10~'° seconds at one atmosphere pressure. 
It is assumed in these calculations that only one 
relaxation time is involved (an assumption without 
an apparent foundation other than that it greatly 
simplifies calculation) and that the effects of Stoke- 
Kirchhoff absorption and relaxation are superposed 
linearly, a procedure described by Trurspety [15] 
as “absolutely wrong” in his paper on hydrodyna- 
mical theory. TRuESDELL points out that the effects 
might add linearly if they are of different orders 
of magnitude. It is noted by Parker et al. [7] that 
there is a serious doubt as to whether or not the re- 
laxation equation used is in any case applicable to 
a process involving rotational energies. 

d) The fourth possible cause of the deviations 
from Stokes-Kirchhoff absorption is the invalidity 
of the equation itself. The equation is derived by 
adding the separate effects of viscosity and heat 
conduction, but the same formula can be obtained 
as a first order approximation to the series solution 
of the Kirchhoff-Langevin [15] propagation equa- 
tion, an equation exact so far as the effects of vis- 
cosity and thermal conductivity are concerned, The 
approximate solution seems valid over the range of 
frequencies used in this series of measurements (see 
TruEsDELL [15], p. 698). Both derivations leave in 
doubt the question of the viscosity to be used since 
the shear and bulk viscosities both appear in the 
result and there is very little experimental evidence 
as to whether or not the bulk viscosity differs from 
zero. Kinetic theory predicts that for a monato- 
mic gas it should be zero, but at the moment there 
appear to be no experimental values for bulk visco- 
sity available in nitrogen and oxygen. If the bulk 
viscosity is not zero the Stokes-Kirchhoff equation 
takes the form: 


2 Se ie 


oc 


[ Rn itl) “| 


C 
where R is greater than 4/3 and from our results 


we can compute R= 1.79 for oxygen and 1.96 for 
nitrogen. 
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It is not at the moment clear whether the relaxa- 
tion theory (possibly in an improved form) or the 
bulk viscosity hypothesis or perhaps some other 
mechanism offers a correct interpretation of the 
observed facts. 
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BEITRAG ZUR DARSTELLUNG DES SPALTFELDES 
BEIM MAGNETTON 


UND ANWENDUNG AUF DEN WIEDERGABEVORGANG 


von G. ScHWANTKE 


Mitteilung aus der akustischen Abteilung des Heinrich-Hertz-Instituts 
fiir Schwingungsforschung, Berlin-Charlottenburg 


Zusammenfassung 


Mit Hilfe der Greenschen Methode zur Lésung des ersten Randwertproblems der Poten- 
tialgleichung wird eine explizite und yerhiltnismaRig einfache Darstellung des Potential- 
und Feldyerlaufs gewonnen. Der Potentialverlauf iiber dem Spalt wird dabei als linear an- 
genommen. Das gleiche Ergebnis wird noch auf einem anderen Wege unter der Verwendung 
der Fourier-Transformation erzielt, wobei sich einfache Integraldarstellungen des Feldver- 
laufs ergeben. Diese lassen sich auf kompliziertere Randbedingungen yerallgemeinern. Fiir 
die Behandlung des Wiedergabevorgangs nach der Reziprozitatsmethode erweisen sich diese 
Darstellungen als besonders geeignet. Es zeigt sich, da} die Integranden bereits die Wieder- 
gabespaltfunktion im wesentlichen reprasentieren. Auf Grund dieser Eigenschaft ergibt sich 
eine eindeutige Beziehung zwischen dem Potentialverlauf itiber demrSpalt-und den Null- 
stellen der Spaltfunktion. 


Summary 


With the help of Green’s method for the solution of the first boundary problem of the 
potential equation an explicit and simple demonstration of potential and field distribution 
is obtained. The potential distribution over the gap is taken as linear. In another method 
using Fourier transforms a simple integral demonstration of the field distribution is ob- 
tained which can be generalised to more complicated boundary conditions and can be 
adapted to reproduction phenomena using the reciprocity technique. The integrals give the 
function for the reproduction in the gap, a property which leads to a single-valued relation 
between the potential field over the gap and the null-points of the gap function. 


Sommaire 


En appliquant la méthode de Green pour la résolution du premier probléme de la valeur 
limite de l’équation du potentiel, on a obtenu une représentation explicite et relativement 
simple de la variation du potentiel et du champ. On a supposé 4a cet effet que la variation 
du potentiel est linéaire dans |’entrefer. On est arrivé au méme résultat par une autre voie, 
en utilisant la transformation de Fourier, ce qui donne des représentations intégrales 
simples de la variation du champ. On peut les généraliser aux cas de conditions aux limites 
compliquées. Ce mode de représentations convient particuliérement bien pour l’étude du 
phénomeéne de reproduction par la méthode de réciprocité. On trouve que les intégrands 
représentent déja essentiellement la fonction entrefer de reproduction. Il résulte de cette 
propriété une relation univoque entre la variation du potentiel dans l’entrefer et les points 
zéro de la fonction d’entrefer. 


1. Einleitung 


Der Feldverlauf tiber dem Spalt beherrscht den 
Mechanismus von Aufsprech- und Wiedergabevor- 
gang. Eingehende Arbeiten iiber diese Probleme 
haben daher Feldberechnungen an den Anfang ge- 
stellt (ScummsBaver [5], Westmisze [6], Greiner 
[9]). Die Struktur des iiblichen Kopfspaltes lat 
eine ebene Behandlung als gute Naherung zu, und 
so wurde als Berechnungsgrundlage durchweg die 
Methode der konformen Abbildung benutzt (Punt- 
mann [2], Wesrmusze [6], Greiner [7], [8], [9]). 
Da die Spaltbegrenzung als unendlicher Polygonzug 
darstellbar ist, ]aBt sich die erforderliche Abbildungs- 


funktion nach dem Schwarz—Christoffelschen Verfah- 
ren gewinnen. 

Allerdings ist dabei die Giiltigkeit der Potential- 
gleichung im ganzen Inneren des Gebietes zu for- 
dern, d.h. die Rechnung gilt nur bei Abwesenheit 
eines Bandes der Permeabilitat u4p41. Diese Be- 
schrankung versucht Gremer [7] zu umgehen, in- 
dem er lings einer Linie y=h — der oberen Band- 
grenze — H,=0 fordert, was leicht durch Spiege- 
lung erreichbar ist. Allerdings kann die Unstetigkeit 
der y-Komponente der Feldstarke langs — 2)... + 2% 
— der unteren Bandgrenze — nicht beriicksichtigt 
werden, was die Ergebnisse doch als nicht ganz 
sicher erscheinen 1aBt. 
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Abb. 1. Idealisiertes Spaltmodell. 


AuBer der Beschrankung auf uwgp=1 bzw. Ab- 
wesenheit des Bandes haftet der Methode noch ein 
weiterer Nachteil an: es gelingt nicht, Potential- oder 
Feldstarkenverlauf als einen geschlossenen analyti- 
schen Ausdruck oder durch eine im ganzen Gebiet 
giltige Entwicklung darzustellen. Man ist vielmehr 
auf punktweise Konstruktion mit Hilfe graphischer 
Verfahren angewiesen. Das hat folgenden Grund: 

Der Potentialverlauf g(w) sei bekannt in einer 
komplexen w-Ebene, in der z. B. ein unendlich aus- 
gedehnter Plattenkondensator liege; gesucht sei er 
in einer z-Ebene, die von der zu untersuchenden 
Spaltform berandet sei. Da die Potentialgleichung 
bei einer konformen Abbildung invariant bleibt, ist 
dazu nur die Abbildung w=w/(z) notwendig, so daB 
man bilden kann 


p(w) =p[w(z)]=y(z). 
Die Schwarz—Christoffelsche Methode liefert aber 


gerade die Funktion z=z(w), deren Umkehrfunk- 
tion sich zumeist nicht auf elementare Weise gewin- 
nen laBt. 

Fiir die Analyse von Aufsprech- und Wiedergabe- 
vorgang ist es jedoch wiinschenswert, eine explizite 
Darstellung der Feldstarke H (x,y) zu gewinnen, die 
bei relativer Einfachheit das wesentliche des Feld- 
verlaufs reprasentiert. 


2. Darstellung nach der Greenschen Methode 
2.1, Allgemeines 


Fiir die Losung der Potentialgleichung Ag = 0 in 
einem Gebiet J’ mit der Berandung s und den Rand- 


werten Ps gilt 
1 OG 
: [ Ps mds (1) 
x n 


(Normale in das Innere des Gebietes gerichtet.) 
G ist die Greensche Funktion des Gebietes J” und 
la8t sich nach folgender Vorschrift bestimmen: 
f(z,z’) leiste eine Abbildung von J” auf den Ein- 
heitskreis derart, da8 z in den Nullpunkt kommt. 
Dann ist 


G= —In|f(z,2’)|. (2) 
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Zur Gewinnung der Abbildungsfunktion f(z, z’) ist 
das Schwarz—Christoffel-Verfahren ebenfalls nicht zu 
gebrauchen, da es wieder gerade die Umkehrfunk- 
tion liefert. Damit ware also zunachst nichts ge- 
wonnen. 

Nun ist aber zu bedenken, dafi fiir die Zwecke 
des Magnettons nur das Feld oberhalb des Spaltes, 
nicht aber im Spalt selbst interessiert. Es liegt daher 
nahe, das Gebiet so zu vereinfachen, da zwar einer- 
seits der eigentlich interessierende Raumteil des 
Spaltfeldes erfaBt wird, andererseits aber G eine 
recht einfache Struktur erhalt. Beide Forderungen 
werden erfillt, wenn wir im Fall des unbelasteten 
Kopfes, bzw. 43=T1. als neues Gebiet die obere 
Halbebene wahlen, wobei zunachst noch offen ist, 
wie die Potentialvorgabe langs der neuen Berandung, 
d.h. der x-Achse, auszusehen hat. Eine solche Be- 
schneidung des Gebietes ist eben deshalb méglich, 
weil Gl. (1) Felder mit beliebigen Randwerten zu 
berechnen gestattet, im Gegensatz zur Methode der 
konformen Abbildung, die nur fiir stiickweise kon- 
stante Randpotentiale geeignet ist. 

Die Abbildung der oberen Halbebene auf den 
Kinheitskreis wird geliefert durch 


, 

Z—2 

f(z,z) =, 
Z-—2 


(Der * bezeichnet die konjugiert-komplexe Grofe.) 


Damit ist 


z =x +iy’ Integrationspunkt, z=2+iy Aufpunkt. 
Mit Gl. (1) folgt schlieBlich 


++o0 
He yin Oe 
ean=" | oe See (3) 


2.2. Ermittlung des speziellen Randpotentialverlaus 


Langs der Spiegelflachen als idealen magnetischen 
Leitern ist anzusetzen 
—%: s(x) = —,/2, 

+2%...+0: (zx) =+9,/2. 

Fir das Ubergangsgebiet —x ... +2») muB die L6- 
sung nach der Methode der konformen Abbildung 
herangezogen werden. Der graphisch konstruierte 
Verlauf wird dann zweckmafig interpoliert. Das 
Verfahren soll kurz skizziert werden. 

Aus Symmetriegriinden geniigt es, die linke Halfte 
der Spaltform von Abb. 1 zu betrachten. 

Die Abbildung der ¢-Ebene auf die z-Ebene derart, 
da8 sich Punkte gleicher Indizes entsprechen, liefert 


lous, Cones 
z=] in. 
E ie VE+1 } o 


wie verschiedentlich ausfihrlich dargelegt worden ist 
(Westmize [6], Greiner [7]). Dabei ist die halbe 


— WH... 


> 
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y 
f v= Im (2w+In2=2) — 
2- Ebene | w+ 
2v 
= are tan (—-—_—____]} —-2v—ax 5b 
Yj (, Sy 5) wee) 
slilmeiae? 22 x . Yom ° 
und nach Gl. (5a) ein zugehdériges Potential. Den so 
gefundenen Potentialverlauf zeigt Abb. 4. 
\Y 
2, yy 
n Wi éNaherung (6b) 


Abb. 2. Abbildung des Halbspaltes auf die obere Halb- 
ebene. 


Spaltbreite zu a gewihlt. ZweckmaSig wird eine 
weitere Abbildung 

VE=Vitin=w=utiv 
angeschlossen. Wir suchen dann die Bildpunkte der 
x-Achse lings z....0, d.h. der oberen Spaltkante, 


in der w-Ebene, fiir die wir jetzt den Potentialver- 
lauf ermitteln wollen: 


Im (z) = Re (20410) =0, 
+1 


2u+In w-1| _4 
w+) 
oder r,/re=e"; 1ry=|w+1|, ro=|w—1]. 


Da fiir Apollonische Kreise r,/rz=const gilt, geben 
deren Schnittpunkte mit zugehdrigen Parallelen zur 

v-Achse im Abstand wu das gesuchte Bild (Abb. 3). 
Andererseits ist das Feld mit der zugehdrigen 
Randbelegung in der ¢-Ebene leicht bestimmbar. Es 

gilt (siehe auch Greiner [7]) 
P= es arc tan i : a) = = 
2 

=%9 are tan(“ —” Zahedi 
2x — Quv 4 


Somit bestimmt jedes Wertepaar u, v der Bildkurve 
in Abb. 3 einen Ort lings der oberen Spaltbegren- 
zung: 


(5a) 


W2 05 
Abb. 3. Bildkurve der oberen Spaltkante. 


Naherung (6a) 


Abb. 4. Potentialverlauf tiber dem Spalt. 


2.3. Lésung von Gl.(3) mit Hilfe linearer Néherun- 
gen fur das Randpotential 


Die Interpolation des Verlaufs in Abb. 4 geschieht 
nun zweckmafig mit Hilfe solcher Funktionen, fiir 
die das Integral (3) elementar lésbar ist, z. B. also 
mit Polynomen in x. Um nicht zu unhandliche Aus- 
drticke zu erhalten, beschranken wir uns auf die ein- 
fachste Naherung: den linearen Verlauf. Es ist dabei 
zu bedenken, das bereits der Verlauf in Abb. 4 durch 
die praktisch nicht realisierbare ideal scharfe Spalt- 
kante die tatsachlichen Verhaltnisse nicht ganz rich- 
tig wiedergibt. Eine Abrundung der Kante wiirde die 
Ecken im Potentialverlauf an den Stellen x= +2 
ausbiigeln und dadurch der linearen Naherung ent- 
gegenkommen. Auferdem wirkt sich der Fehler der 
Naherung sowieso nur in der Nahe der oberen Spalt- 
begrenzung aus. 

SchlieBlich aber kann man neben der Naherung 
des linearen Potentialverlaufs tiber dem Spalt mit 
stetigem Anschluf an das Potential der Spiegel flache 


Gps ee NT ae ee (Ba) 
2° Fo 
eine lineare Interpolation der Form 
] 
pane —%<r< + 2% 5 eal | (6b) 
2 ¥ X% 


verwenden, die sich dem wahren Verlauf besser an- 

pat und die scharfen Spaltkanten als Unstetigkeiten 

des Potentials an den Stellen + 2%) beriicksichtigt. 
Im Fall (6a) ergibt die Integration von Gl. (3) 


p(x, y) -£o(1 - ) arc tan 70 4 
2 y 
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und daraus erhalt man die Feldstarke 


ibe oy end Ue ( arc tan 70 — 


Or” Qa, Yy y 
(8 a) 
ius op Se Tt: (9a) 
OY) > Qian! hp 
mit 7,2= (r+)? +y?, ro? = (x—%)?+y?. 


Im Fall (6b) ergibt sich entsprechend 


p(x, y) = Po_ (r+ =| arc tan ano a 
2ny X y 
+(7 _ = arc tan eed + Y In | (7b) 
\ Xo y a. ee 
und daraus 
H,=— 9 _| are tan OS are tance 
2Y Xo y Y 
1 1 
xy (y no, a ey 
H, « Con eee +ag(n— 1) (24 Zo) 
207 Xp rs Ta ry 
(9b) 


Es ergeben sich — insbesondere im Fall (6a) — 
bemerkenswert einfache Ausdricke fur die Feld- 
starke, die aber den Feldverlauf im wesentlichen 
richtig wiedergeben (siehe Abb. 5 und 6). Man tiber- 
zeugt sich auch leicht, da Potentialgleichung und 
Randwerte erfullt werden. 


1/9 


Abb. 5. Feldverlauf in der Naherung (6a) nach Gl. 
(8a) und (9a). 


Fir groBe Entfernung vom Spalt in x-Richtung 
mit y <z ist der Feldverlauf 


Hy~ —10 9 


7 Hy = © : 
2) eR! Po 


Be 7 


unabhangig von der Spaltbreite der gleiche wie bei 
beliebig kleinem Spalt. 
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Abb. 6. Feldverlauf in der Naherung (6b) nach Gl. 
(8b) und (9b). 


Die Methode ist auch anwendbar, um eine Nahe- 
rungslosung fiir den Fall des aufliegenden Bandes 
hoherer Permeabilitat zu gewinnen, wobei angenom- 
men wird 

1l< MB < Kopf + 


Man geht vom Gebiet des Parallelstreifens aus, wo- 
bei man H,=0 langs der oberen Bandgrenze y=h 
durch Spiegelung erreicht analog zu Greiner [7]. 
Uber den Potentialverlauf lings der x-Achse zwi- 
schen —2) und +2, la®t sich genaues in diesem 
Fall nicht ohne weiteres aussagen. Die Tendenz geht 
jedenfalls bei hoherer Bandpermeabilitat vom linea- 
ren Verlauf weg in Richtung einer starkeren Aus- 
bauchung, da die Feldlinien mehr in das Gebiet der 
hoheren Permeabilitat gezogen werden und damit 
starker aus dem Spalt herausquellen. Naherungs- 
weise kann man den Verlauf durch 


Po sinh (2/9) 
2 sinh 1 


Ps(X) = 


ty os 


Abb. 7. Naherung des Feldverlaufs fiir ug > 1. 
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_darstellen. Damit ist mit der zugehérigen Green- 
schen Funktion eine elementare Integration méglich. 
Die Formeln werden allerdings schon recht kompli- 

_ziert, so das wir hier nur die graphische Darstellung 

-angeben wollen (Abb.7). Fiir grofe Entfernung 
vom Spalt in x-Richtung erhalt man einen Feldabfall 
H ~e~!*'! fiir festes y . 


3. Felddarstellung mit Hilfe der Fourier- 
Transformation 


Wir beschranken uns auch im weiteren auf das 
Gebiet der oberen Halbebene. Diese kann als Grenz- 
fall eines Rechteckgebietes aufgefaft werden. Daher 
ist es sinnvoll, die gesuchte Losung der Potential- 
gleichung 4~=0 nach den ,,Kigenlésungen“ des 
Rechteckgebietes zu entwickeln. 


2 
Cee oF (2) Ciy). 


Oa? Oy? 
Das fiihrt in bekannter Weise auf Losungen der 
Form 


pret ety 


4 ist zunachst beliebig, doch beschranken wir uns 
auf reelle Werte. Lineare Superposition fihrt dann 
auf den allgemeinen Ansatz 


+00 
pay) = f e140) ell + Ay(A) e] dd. (10) 
Auf dem Rande gilt : 
(2,0) =9.(2) = [ eS) ai da) 
mit der Umkehrung S4 
S(a) ao foo.) ee CIT h) 


Die Willkiir in der Wahl der Randwerte @y ist hier 
jedoch dadurch eingeschrinkt, dal die Integrale kon- 
vergieren miissen. Dafiir geniigt es, wenn ys und S 
absolut integrabel sind. 

Ein wesentlicher Vorteil gegeniiber dem Verfahren 
in Abschnitt 2 besteht darin, daB die Giiligkeit der 
Potentialgleichung nicht mehr im ganzen Gebiet der 
oberen Halbebene gefordert werden mu, sondern 
die Giiltigkeit in Streifen parallel zur x-Achse ge- 

-niigt. Man kann damit also Schichtungen beliebiger 
-Permeabilitit behandeln. Dabei sind in jeder Schicht 
Darstellungen der Form (10) anzusetzen. Langs der 
Grenzen y=const treten je zwei Ubergangsbedin- 
gungen auf — namlich die Stetigkeit des Potentials 
und der Normalkomponente der Induktion —, die 

die Koeffizienten A,, A, und B,, By verkniipfen. 
Dazu kommt noch die Potentialvorgabe lings y =0 
und eine Endlichkeitsbedingung fiir y— © . Damit 
ist das Problem eindeutig bestimmt. 
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3.1. Einfachstes Beispiel: Feld ohne Band 
Beschranken wir uns zunachst noch auf den be- 
reits im Abschnitt 2 behandelten Fall des einheit- 
lichen Mediums. Dann veriangt die Endlichkeits- 
bedingung fiir y—> co eine Darstellung der Form 


+00 
(x,y) = [esa e lly d2, 
Tit Beet (12a) 
S(A) = wa [e*o.c@) dz. 
2% 


Man kann daraus eine rein reelle Darstellung gewin- 
nen, indem man ausnutzt, daB g.(x) und damit S(A) 
ungerade Funktionen sind: 


co 


p(x, y) = [sing x)e ¥Yoa(Ad) di, 


2 (12b) 


9 co 
o(A) = = f sin x) s(x) dx 
1 6 
0 
mit o(A) =2iS(A). Daraus gewinnt man dann leicht 
die Feldstarke 


peed bet Pwo) die eds) 


Hy 
0 


Bei der Wahl der Randwerte zur Bestimmung von 
o(A) ist die erwahnte Bedingung der Konvergenz 
zu beachten. Wir setzen daher zunachst eine endliche 
Spiegelkante an, wie es ja den tatsachlichen Verhalt- 
nissen auch besser entspricht. Als naheliegend neh- 
men wir hinter der Kantenlange / ein exponentielles 
Abklingen des Potentials an, das sich stetig differen- 
zierbar anschlieBt (Abb. 8): 


OS Pe ps(a) = 10 * 

2 2X 
% SxS p(t) =", (14a) 
BS His iea 3 p(x) = fa(e—2) +1). 44); 


fiir «<0 entsprechend symmetrisch. 

Natirlich kann man an Stelle von Gl. (14 a) tiber 
dem Spalt im Sinne der Niherung (6b) auch an- 
setzen: 


O<r<m: Y= sd (14b) 


2 ay 


A, 


Abb. 8. Randpotentialvorgabe zur Gewinnung der In- 
tegraldarstellungen. 
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Im iibrigen gilt die gleiche Vorgabe wie Gl. (14a). Mit 
Gl.(14a) erhalt man fiir die zweite Gleichung (12b): 


o(A) =0~(4) + 0a,i(A) , 


res Go sind x (15) 
tA Azaty 

tae » 240% sin Al + a? (2? — a?) cosAl 

Re a (a2 + 22)? 


mit o6(A) +0 fir 24> 0. 
Gl. (13) bildet mit Gl. (15) schon einen recht 


komplizierten Integranden. Nun interessiert aber 
fiir das Spaltfeld der Einflu8 der auBeren Spiegel- 
kanten meist sehr wenig. Die Ausdriicke lassen sich 
vereinfachen, wenn man den Ubergang a— 0 fir 
den unbegrenzten Spiegel unter das Integral ziehen 
kann. 

Das ist in der Tat moglich. Bildet man namlich 


co 


lim p(x, y) = lim lim yee Bie e Vio( Ay dh = 
a0 


a>) «+0 


e-¥ lim o(A) dd= 


a>0 


Tab gea( A) aA 


sin (A x) 


l| 
rans bea 


so existiert die rechte Seite, obwohl o,,(/) fiir 2=0 
einen Pol besitzt. Man kann daher schreiben 


Dos (tof) = [ sina) e¥ ax0(A) di. 
0 


Fur die Feldstarke erhalt man in gleicher Weise fir 
den Grenzfall des unbegrenzten Spiegels 


-++oo 
—HAy _ 1 ALY Ahh A 
oes 2 a 7 ox (2) di. (16a) 
bzw. 
= Havel) eos a 
= we (Aa) ema Moe (A)rds (16b) 
lake lPecrin 


mit Gl. (15) oder (18). 
Die Darstellung (16b) ist identisch mit Gl. (8 a) 
bzw. (9 a), wie ja auch zu fordern ist. Die Kombina- 

tion 
H, +1 H, Sayin Po | eae: anes di. 


1 Ly 
0 


(17) 


stellt naémlich gerade eine Laplace-Transformation 


co 


Ao(A) 


Bigs) e | 
+H, | |cosh (2 h) + (4g —1)-sinh(2 a)- sinh (A d) sin 
0 
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dar von der Form 


foe) 


Rak: sin ¢ 
t 


0 
Eine Aufspaltung in Real- und Imaginarteil liefert 
genau die Behauptung. Diese Ubereinstimmung dient 
als Bestatigung dafiir, daB die Grenziiberginge keine 
Fehler enthalten. 
Vollig analoge Darstellungen gewinnt man natiir- 
lich auch fiir die etwas abgednderte Randvorgabe 


(14b). Man hat in Gl. (16) lediglich 


——dt=arccots. 


PC Neons peu ss +(y—1) cosda| (18) 


YL ALo 
zu setzen. 

Die Darstellungen, insbesondere Gl. (16) unter 
Benutzung von Gl. (15), sind von erfreulicher Ein- 
fachheit. Sie konnen von grofem Vorteil sein, wenn 
der Feldverlauf weiteren Rechnungen, insbesondere 
weiteren Integrationen zugrunde gelegt wird. Kin 
Beispiel wird in Abschnitt 4 gegeben. 


3.2, Feld innerhalb eines Bandes mit 43 >1 im Ab- 
stand a von der Spiegelflache 


In diesem Fall versagt die Methode der konfor- 
men Abbildung, da an den Bandgrenzen statt der 
Giltigkeit der Potentialgleichung eine Unstetigkeit 
der Normalkomponente der Feldstarke zu fordern 
ist. Wir haben also je eine Losung der Form (10) 
anzusetzen fur die Gebiete I, II, II] (siehe Abb. 9). 
Ohne vom Prinzip her dazu gezwungen zu sein, ver- 
einfachen wir das Problem dadurch, da wir an 
Stelle des Ansatzes fiir das Gebiet II fiir die obere 
Bandgrenze 

Ipn/Oy =0 fur y= h 


verlangen. 


I , 
: 


Qs tae 


N 


Abb. 9. Zur Bestimmung des Feldes in einem Band im 
Abstand a. 


Bei einer wesentlich von Eins verschiedenen Band- 
permeabilitat, also etwa ug>5, fallt der Fehler 
sicherlich nicht sehr ins Gewicht, zumal fiir a>0. 

Fur den unteren Rand gilt die gleiche Vorgabe 
wie friher. Die Einzelheiten der Rechnung seien 
iibergangen. Man erhalt schlieBlich 


£05 (3 2)-008P 1 (h— pia. (19) 
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Die Auswertung der Integrale gelingt nun nicht mehr 
in einfacher Weise. Mit Hilfe der Residuenmethode 
mogliche Reihenentwicklungen konvergieren verhalt- 
nismafig schlecht, so daB man praktisch auf nume- 
rische Methoden angewiesen ist. Das soll hier nicht 
im einzelnen ausgefiihrt werden. Wie im Abschnitt 4 
gezeigt wird, besitzt Gl. (19) unmittelbare Bedeu- 
tung fiir den Wiedergabevorgang. 


4. Anwendung auf den Wiedergabevorgang 


Das Problem des Wiedergabevorgangs wurde erst- 
mals in korrekter und befriedigender Weise durch 
Westmisze [6] gelost. Die im Abschnitt 3 beschrie- 
bene Methode erweist sich nun fiir dieses Problem 
als besonders geeignet. Sie fiihrt auf einfache Nahe- 
rungen fiir den Frequenzgang und gestattet weitere 
Verallgemeinerungen und einen vertieften Einblick 
in die physikalischen Zusammenhange. 

Nach einer einfachen Reziprozitatstiberlegung von 
Wesrmisze gilt fir den magnetischen Flu8 © durch 
den Horkopf 


o- | P-grad) dr (20) 
Band- 
ausdehnung 
bzw. 
D=D,+D,= 
h +oco 
aay) ay f to] SE Paley) 4 0 Pile) 
a eae: ‘ 


Es bedeuten 


_ 


P die (vektorielle) Bandmagnetisierung, 

g das Spaltfeldpotential fiir das reziproke (d. h. 
das in den Abschnitten 2 und 3 berechnete) Feld, 
bezogen auf einen Wicklungsstrom 1, 


b die Breite des Bandes. 


Das Verfahren von Scuwarz [3] beriicksichtigt 
bei der Integration der Bandmagnetisierung die 
Gewichtsfunktion in Form der Feldverteilung nicht. 
Dagegen laBt sich der etwas modifizierte Ansatz von 
Scummpaver [4] fiir den betrachteten Spezialfall 
auf Gl. (20) zuriickfiihren. 

In unserem Fall nun ist die Feldstarke als Fourier- 
Integraldarstellung von der Form 


+oo 
—Oy/dx_ ere? dd 


= : (21) 
—Op/ey J Fyly4) 
gegeben. Auferdem setzt man in Gl. (20) 
P,(x,y) =Poz cos(kx—t) , (22) 


P, (x, y) = Poy cos(kx—@ t+) 


ein, wobei die Magnetisierungsamplitude Py ebenso 
wie die Phasenlage 6 zunichst als unabhangig von y 
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angesehen wird — eine Voraussetzung, von der man 
sich bei spaterer genauerer. Kenntnis des Aufsprech- 
vorgangs leicht befreien kann. 

Die Kombination von Gl. (20) und GI. (21) fiihrt 
zunachst auf ein dreifaches Integral, zweifach iiber 
das unendliche Intervall erstreckt. Nun besteht aber 
— F(A) als stetig vorausgesetzt — auf Grund der 
Fourier-Darstellung einer Funktion die Identitat 


+oo +00 
Gays i, F(a) a | eile de. 
2a 


Dadurch kann man das dreifache Integral auf ein 
einfaches reduzieren. Bedenkt man noch, da bei 
unserer Symmetrie der Randvorgabe gilt 


Fr(k) =Fy(—k), Fy(k)=—F,(—&), 


so folgt 
h 
®, oh Bon cos Oe F,(y, k) dy. (23) 
Dy — Po, sin(wt—0) J F,(y,k) 
a 


Die Gleichungen (23) zeigen, dafi — abgesehen von 
der unwesentlichen und leicht durchfithrbaren In- 
tegration tiber die Dickenausdehnung des Bandes — 
der Integrand der Fourier-Darstellung des reziproken 
Spaltfeldes die Spaltfunktion des Wiedergabevor- 
gangs bereits reprasentiert. Damit sind aber die Dar- 
stellungen von Abschnitt 3 von neuem Interesse, un- 
abhangig davon, ob der Feldverlauf explizit ausge- 
wertet werden kann oder nicht. 


4.1. Gewinnung einfacher Naherungen fiir die Spalt- 
funktion 


Fir den Fall wg =1 entnimmt man Gl. (12) 


Fe(y,k) 3 
= | eM ho(k). (24) 
F,(y, k) =, 
Das ergibt mit Gl. (23) 
—kd 
De eth dar a — bo (k) 
kd 
x [Porcoswt+Po,sin(wt—d)]; (25) 


(a Abstand vom Spiegel, b Breite des Bandes, d 
Dicke des Bandes). 

Da die in o(k) eingehende Konstante @ gemal 
Gl. (20) auf den Wicklungsstrom Eins bezogen sein 
soll, hat man zu setzen 
Rsp 
Rise? 

(w Windungszahl, Rs, magnetischer Widerstand des 
Spaltes, R,., magnetischer Widerstand des Gesamt- 


Po=w 


kreises). 
Im einfachsten Fall a=0 (aufliegendes Band) 
und kd <1 (diinnes Band) ist der Wiedergabe- 
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frequenzgang allein durch den Ausdruck k o(k) ge- 
geben, wobei o gema Gl. (12) die Fourier-Trans- 
formierte des Randpotentials fiir das reziproke Feld 
darstellt. Fiir die beiden Naherungen des Potential- 
verlaufs tiber dem Spalt Gl. (6a) und Gl. (6b) 
(siehe auch Abb. 4) liefert Gl. (15) bzw. Gl. (18) 
fiir den Grenzfall des unbegrenzten Spiegels zwei 
genaherte Spaltfunktionen 


k 600(k) = %9 SinVedg) (26a) 
Tea IC 
k 600(k) = Syeaess + (y—1) cox ka | (26 b) 
Tanya enone y 
y= 151 


Gl. (26a) entspricht genau der elementaren Spalt- 
funktion, wie sie fiir erste Abschatzungen gebrauch- 
lich ist (siehe auch Gucxensure [1]). Sie kann nach 
Wesrmisze unter der Annahme reiner Langsmagneti- 
sierung und fir ein homogenes Feld auf einfache 
Weise abgeleitet werden. Sie gilt aber auch fiir den 
Fall reiner Quermagnetisierung. Gl. (26b) bedeutet 
eine Verbesserung der Naherung, was sich am deut- 
lichsten an den Nullstellen zeigt. Diese sind unabhan- 
gig von der zuletzt gemachten Annahme kd <1, die 
fiir technisch realisierte Falle nicht durchweg gilt. 

Fur Gl. (26a) liegt die erste Nullstelle bei kxy =a, 
bei Gl.(26b) dagegen bereits bei k2)~0,9 2, was 
dem Experiment besser entspricht. Will man dieser 
Tatsache Rechnung tragen, trotzdem aber den ein- 
fachen Verlauf (26a) beibehalten, so hat man eine 
um etwa 10% groBere Spaltbreite x) einzusetzen. Das 
fiihrte zur Unterscheidung der geometrischen von 
einer “effektiven“ Spaltbreite (siehe Abschnitt 4.3). 

Den Verlauf der Spaltfunktionen Gl. (26a) und 
Gl. (26b) zeigt Abb. 10. 

Fiir den Spiegel endlicher Lange ist nach Gl. (15) 
6x0.(k) um 6a,1(k) zu erganzen. Das fiihrt zu einem 
Tiefenabfall und einer Welligkeit bei den tiefen Fre- 
quenzen, die auch vom Experiment her wohlbekannt 
sind. 


4.2. Eine Verallgemeinerung der Wiedergabespalt- 
funktion 


Auf Grund der Feldberechnung GI.(19) kann man 
nun auch eine Bandpermeabilitat > 1 bei Abtastung 
im Abstand a erfassen. Man erhalt aus Gl. (23) mit 
Gl. (19) 


®, 
=nbd- 


2, 


ko(k) 
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4x 
u—e 
kig-u, y=") 


sin u 


[ J Sey ticosu] 


Abb. 10. Vergleich der aus den Naherungen (6a) und 
(6b) folgenden Spaltfunktionen; kd <1; 
a=0', 


Nahe von Eins gilt Gl. (26). Prinzipiell gestattet die 
angegebene Methode ‘auch die Konstruktion einer 
Wiedergabespaltfunktion, die beide Falle enthalt, 
der Gewinn rechtfertigt jedoch den Aufwand kaum. 

Der Vergleich von Gl. (26) und Gl. (27) zeigt, 
da eine hohere Bandpermeabilitat fiir die Wieder- 
gabe ungiinstig ist und einen weiteren Hohenabfall 
bewirkt. In den gebrauchlichen Bereichen mit wp <20, 
ka<1,kd<1 wird aber der Einflu8 im allgemei- 
nen zu vernachlassigen sein. 

Ein indirekter Einflu8 kann durch die Abanderung 
des Randpotentials und damit von o(k) bewirkt 
werden (siehe Abschnitt 4.3). 


4.3. Zusammenhang zwischen Randpotential und den 
Nullstellen der Spaltfunktion 


Wie schon bemerkt, sind die Nullstellen der Spalt- 
funktion eindeutig durch die von o(k) gegeben, da 
die anderen Faktoren nicht verschwinden. o(k) ist 
aber tiber eine Fourier-Transformation eindeutig 
durch das (fiir das reziproke Feld geltende) Rand- 
potential langs y =0 bestimmt. 

Der Vergleich von Gl.(26 a) und GI.(26b) zeigt, 
wie empfindlich die Nullstellen auf kleine Anderun- 
gen des Randpotentials reagieren. Setzen wir ftir 
den allgemeinen Fall eine Fourier-Reihe 


1 a x 
Y;= vat ae Se én sin nx — ): 


0 n=] 0 
~% SaS+ay (28 a) 
fiir das Randpotential iiber dem Spalt an, so folgt 
daraus nach Gl. (12) 


ko x(k) = (28 b) 
a Bo E k “0 _9Onk %y sin k x y Re énn(—I1)* 
1 k xo n=1 (k x)? — (nx)? h 
sinh kd Po, cos mt 
kd 


. (27) 


td Poy sin(wt—0d). 


Gl. (27) enthalt mit Gl. (19) noch die Voraussetzung, daB die Normalkomponente der reziproken Feldstarke — 
an der Bandoberseite verschwindet, ist daher erst fiir hdhere jp geeignet, z.B. up>5. Fiir wp in der 


swe 


a 
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An der m-ten Nullstelle des ersten Terms verschwin- 
den auch alle Summenglieder bis auf das m-te, wel- 
ches dort den Beitrag 
[O—xm Em |} 

liefert. Schreiben wir kx =u und die Nullstellen 

an.( | =0 mit u,=mx+4Anu, 

X% 
so folgt in linearer Naherung (¢,, < 1) 
Amt = (—1)™ m2 x? em, 
bzw. fiir die Nullstellen der Spaltfunktion selbst 
Uy =m ni l+(—1)"xé&_,m] . (29) 

Jede der Nullstellen liefert also einen Fourier-Ko- 
effizienten des Randpotentialverlaufs (fiir das rezi- 


- proke Feld) tiber dem Spalt und ihre Gesamtheit 


damit eindeutig den Verlauf selbst. Die experimen- 
telle Bestimmung der Nullstellen aus dem Wieder- 
gabefrequenzgang ist recht genau moglich. Band- 
eigenschaften wie Breite, Dicke, innere Dampfung 
usw. — abgesehen von der Permeabilitat — haben 
keinen EinfluB. Ebensowenig gehen die Parameter 
des Aufsprechvorgangs ein. Praktisch ist allerdings 
die Ordnung der erreichbaren Nullstellen begrenzt 
durch die Schwachung der Aufzeichnung infolge der 
exponentiellen Banddampfung. Das gilt insbesondere 
bei kleinen Spalten; hier ist jedoch auch die Wieder- 
gabespaltfunktion oberhalb der ersten Nullstellen 
nicht von Interesse. 

Umgekehrt kann man natiirlich bei bekanntem 
Potentialverlauf — etwa nach Abb. 4 — die Nahe- 
rung fiir die Wiedergabespaltfunktion verbessern, 
indem man Gl. (28b) in Gl. (26) bzw. Gl. (27) 
einfiihrt. 

Fiir die technische Anwendung ist die erste Null- 
stelle der Spaltfunktion am .wichtigsten. Die aus 
Gl. (29) folgende Verschiebung zu einem kleineren 
Wert als x driickt sich — sofern man die elementare 
Spaltfunktion (26a). zugrunde legt — als eine 
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scheinbare Spaltverbreiterung 4,2) aus. Fir diese 
folgt in linearer Naherung: 


A,u 


A, 
val) a ee 
7 


X 
Man kann auf diese Weise unmittelbar aus dem Po- 
tentialverlauf auf die effektive Spaltbreite schliefen. 
Die Fourier-Analyse des Verlaufs Abb. 4 fiihrt mit 
é;=0,04 zu einer Spaltverbreiterung um 12%, was 
annahernd in Ubereinstimmung mit der Erfahrung 
steht. Mit zunehmender Bandpermeabilitat, d.h. zu- 
nehmender Steigung der Magnetisierungskennlinie 
im jeweiligen Remanenzpunkt, ist eine wachsende 
Ausbauchung é, zu erwarten und damit eine weitere 
scheinbare Spaltverbreiterung. Auch in diesem Sinne 
wirkt sich also eine hohere Bandpermeabilitét auf 
den Wiedergabevorgang durch eine Benachteiligung 
der hohen Frequenzen unginstig aus. Praktisch spielt 
dieser Effekt allerdings kaum eine Rolle, da die 
obere Frequenzgrenze bei den heute moglichen klei- 
nen Spalten allein durch die Eigenschaften des Ban- 
des bestimmt wird. 


Herrn Prof. Dr.-Ing. L. Cremer danke ich fiir das 
fordernde Interesse an der Arbeit. Mein Dank gilt 
auch Herrn Dr. Guckensure fir wertvolle Hinweise. 


(Eingegangen am 9. November 1956.) 
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Summary 


The theory of Zwrxxer and Kosten for the propagation of sound waves in porous, flexible 
materials is extended by the introduction of known expressions forthe complex density 
(Cranpatt) and bulk modulus (Kosren) of air in cylindrical_pores. 

The complicated solution of the differential equations for the case of a flexible layer 
on a rigid wall is simplified and then discussed. Measurements corroborate the theoretical 
results. An increase of the sound absorption coefficient in the low frequency region is the 
most spectacular effect of the flexibility of the medium. 


Sommaire 


On généralise la théorie de Zwixxer et Kosten concernant la propagation des ondes 
sonores dans les matériaux homogénes, poreux et flexibles en introduisant les expressions 
connues de la densité (Cranpatt) et de la compressibilité complexes (Kosten) de lair 
dans des pores cylindriques. 

La solution compliquée des équations différentielles dans le cas d’une couche flexible 
appliquée a une paroi rigide est simplifiée puis discutée. Les résultats expérimentaux sont 
en bon accord avec la théorie. L’effet le plus spectaculaire de la flexibilité du milieu est 
une amélioration de l’absorption aux basses fréquences. 


Zusammenfassung 


Die Theorie von Zwixxer und Kosten fiir die Fortpflanzung von Schallwellen in homo- 
genen, pordésen, nachgiebigen Schallschluckstoffen wird erweitert durch Einfiihrung von 
bekannten Ausdriicken fiir die komplexe Dichte (Cranpaut) und Kompressibilitat (Kosten) 
fiir Luft in zylindrischen Poren. | 

Die recht komplizierte Loésung fiir den Fall einer nachgiebigen Schicht auf einer starren 
Wand wird vereinfacht und diskutiert. Messungen bestatigen die Theorie. Der wohl am 
meisten in Erscheinung tretende Effekt auf Grund der Nachgiebigkeit ist eine wesentliche 


Steigerung der Schallabsorption bei niedrigen Frequenzen. 


1. Introduction 


It is well known that in the case of plane waves 
travelling through a simple, homogeneous material 
the acoustic properties of the medium can be fully 
described by two complex quantities, e. g. the (com- 
plex) propagation constant y and the (complex) 
wave impedance WW. Another possible set of para- 
meters is the (complex) bulk modulus K and the 
(complex) density @. These two sets of parameters 
are interrelated by the equations y =j @/VK/o and 
W—-VKo. 

In the case of sound waves travelling through the 
air in a homogeneous, porous material with rigid 
skeleton the air is the medium to which the state- 
ments of the first paragraph apply. Its parameters 
in this case can be derived from basic physical pro- 
perties. We refer in this respect to the fundamental 
theories of Kirchhoff (propagation in cylinders at 


high frequencies [1]), Cranpati (expression for @ 
in cylinders [2]) and Kosren (expression for K in 
cylinders [3]), as well as to the phenomenological 
theories of Zwikker and Kosren [4]. 

In many cases the theories mentioned were ap- 
plied to materials of which the skeleton was so flex- 
ible that it seems dangerous to assume that the 
frame remains immobile when a sound wave pene- 
trates into the material. In fact there have been 
published absorption characteristics of low-density, 
flexible, porous materials that cannot be explained 
by the simple “rigid” theory. These materials should 
be dealt with as consisting of two components, two 
media, both taking part into the vibration and 
vibrating more or less independently. A general 
theory for these materials has been drawn up by 
Zwikker and Kosten [5]. Unfortunately the latter 
results are very complicated indeed so that they are 
hardly of any practical value without simplification. 


J ae ee ~*~ o>) Zz* a 
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It is the aim of the present article to derive some 
simplified expressions and to démonstrate their value 
by means of a few experimental results. 

It is astonishing that, although flexible sound ab- 
sorbing materials are used extensively, the literature 
on the theoretical aspects of wave propagation in 
such materials is scanty indeed. Apart from [5] we 
only know of one publication [6] bearing upon this 
subject; it follows the line of thoughts of [5] and 
contains simplified formulae analogous to those of 
the present article; in our opinion there is serious 
doubt as to the correctness of the differential equa- 
tions on which [6] is based so that we prefer to 
base our present work on [5] exclusively. The same 

‘symbols as used there by Zwixxer and Kosten will 

be used here (see list at the end of this paper). It is 
impracticable and unnecessary to repeat here in de- 
tail the theory mentioned. We shall only recall the 
line of thought, thus enabling the reader to under- 
stand at least the meaning of the results arrived at 
in this paper without reading [5]. For a full under- 
standing we refer the reader to the original publi- 
cation. 

If a sound wave penetrates into a porous flexible 
material the air inside as well as the solid frame 
(skeleton) will be set into vibration. The wave pro- 
pagation is governed, therefore, by two equations 
of motion 


3 


eo, at + 8(04— 09), (1) 
SPs 0, 282 + s(o V1) (2) 
and by two equations of continuity 
Ps Ky+ u ma (Ky — Po) oe 
Ov9 (3) 
+ (1—A) Kya 
Sh) (K, —P;) 2214 h Kyo. (4 
. ie ( h) (Ky Py) meth Kx (4) 


If a plane progressive simple harmonic wave is 
to satisfy the four differential equations, all depen- 
dent variables must be of the form exp(j @t—y2), 
where y is the propagation constant of the wave; 
consequently we may replace 3/9¢ by jw and 9/dx 
by —y. The equations (1)... (4) then simplify 
to four linear homogeneous algebraic equations in 
Pi; 1; Pe; V2 having a non-trivial solution only if 
the determinant of the coefficients equals zero, which 
leads to an equation in y having two different com- 
plex solutions for 7, say 7 and 7”. One can there- 
fore distinguish two different wave types at the 


same pulsatance @. Each wave type, if present in 


the air, is also to be found in the skeleton so that 
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at one single pulsatance w four waves can be dis- 
tinguished in one direction and in general also four 
reflected waves. 

When one of the two values for y is substituted 
in the four homogeneous equations in p;, V1, Po, V2 
we can compute all possible ratios between these 
quantities, e.g. p,/v; and ps/vy, i.e. the wave im- 
pedances of the skeleton and the air respectively 
for that y. There are four of them, W,’, V,", W2’ 
and W.’, since there are two media and two wave 
types. 

Also the ratio p;/p. can be computed if y is 
known, which means that the amplitudes of the 
waves in the skeleton and the air for one and the 
same wave type are not independent of each other. 
This comes down to saying that there are only two 
independent waves at the same in one direction 
(and two independent reflected waves). In a specific 
case their amplitudes should be calculated from 
the four boundary conditions of the layer. 

We must refer to the original publications, of 
which a survey is given in [5], for a more exten- 
sive presentation of what has just been given in 
outlines. 

vo, Wy, Wy", Wo, Wo” and other quantities 
can be expressed along the lines described as func- 
tions of the coefficients of the differential equations 
and w. 

Several of these coefficients, e. g. s and K,, are 
rather complicated functions of other constants of 
the material, such as the air resistance, the porosity 
and the structure factor, and of the pulsatance w. 
These relations will be derived in section 2. We 
shall then (section 3) give simplified solutions ap- 
plicable in the case of very porous and flexible 
layers on a rigid wall and describe finally some 
practical results that corroborate the theory (sec- 
tion 4). 


2. The coefficients in the differential equations 
for the wave motion 


In the equations of motion (1) and (2) s is the 
so-called coupling coefficient, a constant of the ma- 
terial, and, therefore, independent of e. g. v, . Equa- 
tion (2) must be satisfied always and everywhere, 
so also at a point in the wave where v, happens to 
vanish for all values of ¢. The equation to be satis- 
fied there is obviously [compare eq. (2) ] 

— Op,/Ox = (0, +5/j @) dv,/9t (5) 
if we confine ourselves to the harmonic steady state 
only. 

The expression in brackets in equation (5) can 
be considered as the generalised complex density 
of the air. Zwrkxer and Kosten [7], with the aim 
of expressing the coupling coefficient s in well- 
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known quantities, compared the equation of motion 
in this form with the more usual form, borrowed 
from the “rigid”? theory [4], containing the struc- 


ture factor k; and the specific air resistance ©, 


namely: 
— Op,/dx = (ks 0p +h? o/j w) Ovo/9t (6) 

thus arriving at 

S=j@0,(k,—1) +h? o. (7) 
There is, however, another approach to the problem 
that yields an expression for s that probably gives 
a somewhat better description of s as a function of 
frequency. This approach is based on the following 
form of the equation of motion instead of (6): 

— Op,/dx =k; ho* Ov./9t, (8) 


where 


Qe" =O |1— - ee (9) 
in which w= //8w ky @/ho. 

o* is Cranpatt’s expression for the complex 
density of a gas in a rigid cylindrical pore [2]. 
Combining (5) and (8) one obtains 


S=s/jo=k,ho* — @5° 


(10) 
0 VK2/x% 4 


Fig. 1. The complex quantity /K,/x P, as a function 
of u (in the first quadrant of the complex plane), 
and the complex quantity V/Oefr/ks 09 as a func- 
tion of w and ks Q)/0,; (as parameter) (fourth 
quadrant). With these curves and a chosen value 
of the complex velocity of sound c” = VKo/oert 
and of the wave impedance W.” = Ko oett. 
The dotted parts of the curves VOort/ks 09 are 
not based on numerical calculations. The curve 


labelled “rigid’’ applies to non-flexible materials 
(Cranpaui [2]). 


In Fig. 1 is shown a.o. Vo*/g) as a function 
of uw, through which 0* is implicitly a function of 
the frequency and the specific air resistance. One 
sees at a glance from Fig. 1 that @* can only de- 
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viate from Q considerably, say by a factor 10 or 
more, for small values of «, i. e. at low frequencies 
and/or for very dense materials (i. e. with high air 
resistance). 

Another coefficient in the basic equations (1) 
...(4) is Ky, the complex bulk modulus of the air 
in the porous material. We can borrow an expres- 
sion for this quantity from [3], viz.: 


K, =» P, | 2 ee 


Buy —j Jo(B uV —j) 


(11) 


where B= V7, Cy/Ay-=~.0.86 for diatomic gases (air) 
and_-is the frequency and resistance parameter 
introduced previously. The quantity /VK./x Py is 
also shown in Fig. 1 as a function of its variable u 
in the complex plane (for air in a rigid cylindrical 
pore). . 

Assuming that the specific stiffness K, of the 
frame is known, e. g. from measurements, as well 
as the static specific air resistance, we know all co- 
efficients of the differential equations (see also the 
list of symbols at the end of the paper), so that we 
can proceed to consider their solutions. 


3. Simplified solutions 


For many practical materials one will find that 
K,/Ks <1, 0,/02>20 and h~1. Such materials 
may be called very porous and flexible. 

As has been stated elsewhere in a slightly dif- 
ferent notation ([8] or [9]), without giving the 
lengthy and tiresome derivation, it can be shown 
by elementary calculation and very careful approxi- 
mation that the impedance of a layer of such a ma- 
terial with a porous surface and backed by a rigid 
wall can be written in the same form as for non- 
flexible materials, namely 


z=W,” cothy’ Ll. 


Since W,” and y” are the wave impedance and the 
propagation constant, respectively, of the air type 
wave in air, the latter equation might give the im- 
pression that the approximation has been so rigo- 
rous that the flexibility of the material is entirely 
neglected; the impedance seems to be governed by 
the properties of the air alone. This is not the case, 
however, as will appear from the following expres- 
sions, in which the quantity Q,¢p is also a function 
of the density of the skeleton, due to the fact that 
the flexibility of the material has been taken into 
account. 


7” =j@Voett/Ke, 
ed ea h Vectt Ko > 


(12) 
(13) 


where Qefs is the effective density of the air in the 
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material, defined as 
1/Octt = 1/ks 0* + 1/0, . 


The expression for @,¢ as implicitly given in [8] 
differs from (14) in so far as the right hand mem- 
ber of 


1/Octt = (02 + S)~1 + (S/0;) (02 +S) 


differs from 

U/kg 0" +1/o,; now g:+S5=k,o%, but S/ (e+) 
is not always equal to unity to a sufficient approxi- 
mation as would be necessary for the two approxi- 
mations to be identical. For k,o* > 0, however 


(14) 


S/(@2+5)—>1 and for k,o* ~@ the values of 
l/o, or (S/0;) (02 +5S)~1 are negligible compared 


with 1/k,0", so one may conclude that (14) is 
sufficiently correct for the present purpose. 

For the frequency tending to zero (9* > ~) we 
obtain 
(12’) 
(13’) 


Wo Vals. 
Wy” =hV0,Ky, 
which are reasonably good approximations for the 
correct expressions for this case 


7 =j OV (0 + @2)/(Ky + Ky) , 
W.” =hV (0; + 02) (Ky + Ko) 
(c.f. for example [4], p. 68). 

The latter equations express the fact that at fre- 
quency zero the coupling between the air and the 
flexible skeleton is so tight that they move together. 
The total density 0;+ 0, and the total stiffness 
K, +Ky govern the propagation. For the frequency 
tending to infinity (0*—> 09) we find that >; —> k, 0” 


~ and so we obtain 


(12”) 
(13”) 


y =joVk,o*/Ks, 
Wo =hVk, 0* Ky. 


These equations mean that the flexible skeleton 
is so inert that it no longer vibrates appreciably; 
air and skeleton are almost completely decoupled. 
The same can be said in the case that 0; > ~; this 
means, however, that the skeleton is immovable. 
We may conclude that the acoustic properties of 
homogeneous porous materials with rigid skeleton 
can be derived from the theory given here by put- 
ting 0,= ©, as was to be expected. 

o* and K, are for cylindrical pores known func- 
tions of « (see equations (9) and (11) and Fig. 1); 
other variables that enter into the problem are h, 
k, and @,, of which h might be put equal to unity 
in sufficient approximation for many materials. 

The constant k, deserves special attention because 
it is rather difficult to define its exact physical 
meaning. It might seem to be introduced to match 
theory with experiment. There must however be 
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something like a structure factor as introduced by 
ZwikkeR and Kosten (see e.g. [10]), because ma- 
terials like e. g. fibreglass certainly cannot be de- 
scribed solely as piles of straws. As long as the 
theory presented here is able to describe the be- 
haviour of very porous and flexible materials suf- 
ficiently it must be possible to measure ‘,. It can 
be shown now that 


k, = (0.375/sin y — 0.360)+ 04/05, 


where @ is given by twice the angle between the 
real axis and the tangent through the origin to the 


locus representing //Q¢f;/Q) in the complex plane, as 
computed from measurements of y” and W,” (re- 


member that V0et;/0)9 = V7" W»/j © Q9 holds). This 
procedure is followed for the measurements given 
in section 4, 

In order to give an insight into the influence of 
the various variables on y” and W,” we have shown 


on Fig. 1 the quantity V Gett/ kes Q9 (in the complex 
plane) as a function of the parameters “ (for fre- 
quency and resistance) and 0,/k;0Q) (for density 
and structure). Moreover one will find in this figure 
the function /K,/x Py for air, i.e. for B=0.86. 
From this it is easy to compute the values for y” 
and W,” as soon as the basic properties (constants) 
of the material under test are known; schematically 


y” =j » Vis (Vectt/ ks Oo) Co (VKo/# Po) 
and W.” =h Vig 09 co(V Oett/ks Qo) (VKo/% Po) - 


Inspection of Fig. 1 shows that the influence of the 
introduction of flexibility is a shift of Voce to the 
real axis and to a smaller absolute value. This shift 
is the more pronounced the lighter the skeleton is, 
and is rather important at low frequencies (small 
values of «). VK»o/xP, does not deviate a great 
deal from unity irrespective of the numerical value 
of 4, which means that what has just been said for 
Vecte applies also to the wave impedance W,”. 
A shift of W,” to smaller values of a more “real” 
character is benificial to a high absorption coeffi- 
cient, so that the conclusion of the theory is that 


at low frequencies an increase in absorption coeffi- 
cient must be expected from the flexibility of the 
material. It seemed worthwhile to verify experimen- 
tally whether this improvement is really observed. 
To this end experiments were carried out with the 
object of measuring y and the impedance of layers 
of flexible materials; in this way W is implicitly 
measured. 


4. Measurements 


In the case of rigid materials the propagation 
constant y and the wave impedance W are usually 
measured with the aid of an acoustic interferometer. 
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The application of this method to the case of flex- 
ible materials meets with the great difficulty that 
on the one hand the flexible sample should fit the 
interferometer tube in an airtight way, whereas on 


the other hand the tube should not hamper the. 


vibrating skeleton of this flexible sample. 

These requirements come to saying that there 
should be no contact between the sample and the 
tube, but the free space between sample and tube 
must at the same time be negligible. 

A vertically placed cylindrical tube filled with 
disks cut from the flexible blanket under test would 
seem at first sight to be the solution to this problem, 
if the diameter of each disk could be one or two 
millimeters smaller than the diameter of the tube. 
Only for short lengths of the test sample however 
this solution is indeed correct, as the density of the 
material would vary considerably over the sample 
due to compression of the lower layers by the upper 
ones if samples were used, of greater length than, 
say, 0.5 m. 

Of course, another solution can be found in the 
Then however 
contact between disks and tube walls cannot be 
avoided. 


normal horizontal interferometer. 


In view of these and other difficulties we applied 
three different methods for the measurement of the 
quantities of interest using horizontal as well as 
vertical wave tubes; a sketch may be found in 
Fig. 2. 

With method A, y” is found directly from meas- 
urements of the sound pressure at various places in a 
very long sample, assuming that p = pp exp(— 7’ 2). 
In method B, y” is computed from the measured 
complex ratio of the sound pressures pp and p; at 
the front and back surface respectively of a sample 
of thickness / backed by an infinite impedance, as- 
suming that py=p,coshy1. Finally, method C 
yields the impedance z which is known to equal 
W,” cothy’ lL. 

The three equations on which these methods are 
based can only be approximations, useful when it 
can be shown that the frame type wave (y’) is 
negligible in the given situations. This has been 
shown indeed and an account wil be published 
later. 

We shall now describe in more detail how the 
measurements were carried out. 

In method A a microphone of the probe tube 
type [11] is moved parallel to the axis of the 
horizontal interferometer tube through a triangular 
canal at the lower side of the disks. The cross- 
section of the probe canal was approximately 1 cm?, 
the diameter of the copper interferometer tube about 
10 cm. The sample consisted of many circular disks 
with a triangular notch (for the canal); the total 
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Fig. 2. The three methods used to measure: 

A) the propagation constant y from 
Px/Po = exp (—7 2) 

B) the propagation constant y from 
Do/pi=cosh yl, 

C) the acoustic impedance z=W cothy/. 

In all cases the diameters of the disks cut from 

a flexible blanket were slightly less than the dia- 

meter of the resp. tubes. 


length was about 1.5 m so as to prevent appreciable 
reflections from the end. Nevertheless accurate 
measurements below 80 c/s were not possible due 
to the small damping. 

In method B the pressures pp and p; were meas- 
ured by means of one microphone attached to a 
short probe tube that could be inserted alternately 
in two small holes drilled in the interferometer tube 
wall at the front and back surfaces of the sample. 
Also with this method measurements at very low 
frequencies (<60c/s) were inaccurate, this time 
because of the errors consequent on the leakage of 
sound through the chink between sample and tube. 


For the methods A and B it was possible to de- 


termine the phase differences rather accurately (to 
+ 2°) by direct reading on the screen of an oscillo- 
scope with circular time base. On this time base, 
derived from the sinusoidal loudspeaker signal, was 
imposed a small pulse, derived from the peak value 
of the microphone signal, thus serving as a “poin- 
ter”. 


Method C involves still more computation: the — 


acoustic impedance (z= W,” coth yl) of a sample 
of known thickness /, air resistance o and structure 
factor k, (derived from previous measurement, e. g. 


with methods A and B) backed by a rigid wall is — 


computed and compared with the values measured 


eS ee ee ee 


Pe ee ne ee 


ae ee) ee 


in the acoustic interferometer. The cross-sectional 
dimensions of the sample were 1 mm less than those 


‘of the interferometer (225 mm x 225 mm). No trans- 


verse waves were observed for frequencies up to 
2600 ¢/s. 

It will be clear that many sources of errors were 
present in the methods mentioned. One can there- 
fore expect that the results of measurements will 
not be very accurate, the third method probably 
being the best. 


Fig. 3. y Colj @ (=c,/c’) for fiberglas as a function of 
the frequency, indicated in c/s in the figure along 
the theoretical curves and near the measured 


points, 
Sample I: 0,;=20kg/m*, 6 =104 kg/m’ s, 
ks=1.2 (estimated) ; 
eee measurements with method A. 
Sample II: 0;=20 kg/m’, o=3-10kg/m3 s, 


ks=2 (estimated) ; 
mm measurements with method A. 


Fig. 4. y" cy/j @ (=cy/c”) for fiberglas as a function 
of the frequency in c/s. 0;=20 kg/m’, 
o=104kg/m*s, k,=2.0 (estimated) ; 

theoretical, 

e measurements with method B. 

——— theoretical curve neglecting flexibility. 


Fig. 5. y” ¢9/j @ (=c/c” ) for tinslagwool as a function 
of the frequency in c/s; 0, =85 kg/m’, 
o=5:105 kg/m? s, ks =2.5 (estimated) ; 

theoretical, 


« measurements with method B, 
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In the Figures 3, 4 and 5 is shown the dimen- 
sionless quantity ¢)/c” for four different materials 
(c’ being the complex sound velocity = j w/y’’). 
Theoretical curves are drawn for comparison. It is 
seen that especially in the lower frequency region 
the differences between “rigid” and “flexible” are 
great, This is of considerable importance to the 
design of new absorbing materials such as e. g. those 
made of foam plastic or sponge rubber. 

In Fig. 6 is shown the dimensionless specific im- 
pedance z/0) cy of a layer of fiberglas backed by a 


flexible 


Fig. 6. z/Q9 cy) for a layer of 10 cm fiberglas as a func- 
tion of the frequency; 0, =25 kg/m’, 
o=3:10!kg/m's, ks=2 (estimated) ; 

theoretical, 

e« measurements with method C, 

—— — theoretical curve neglecting flexibility. 
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Fig.7. The absorption coefficient at normal incidence 
as a function of the frequency as derived from 


Fig. 6. 
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rigid wall together with the theoretical expectation 
for a flexible as well as for a rigid skeleton. The 
measured points lie rather close to the “flexible” 
curve and in general the agreement in our opinion 


is excellent. This is the more important as the me-' 


thod used (C) is supposed to be the most accurate 
one. In Fig. 7, derived from Fig. 6, the rather good 
properties of porous and flexible absorbing layers 


can easily be recognized. 
(Received December 29th, 1956.) 


List of symbols 


More indicates frame type wave or air type wave 
respectively, 

lor 2 as suffix indicates frame or air property res- 
pectively, 

B sq. root of the Prandtl number; B equals 
0.86 for air, 

Cy specific heat of air at constant pressure, 

c complex velocity of sound, 

Co velocity of sound in air, 

y=jqa/c propagation constant, 

y viscosity of air, 

h porosity, i.e. ratio of the volume of the ac- 


cessible holes to the total volume of the po- 
rous flexible material, 


j=V-l, 

J Bessel functions, 

K,, Ks complex specific stiffness of the frame and 
of the air respectively, 

ks structure factor, 

#=C,/Cy ratio of specific heats of air, 

An thermal conductivity of air, 

Lu (8 w ks Qo/h 0)”, 

Re static pressure of air, 

Pix Pe force acting on frame or air per unit area of 


cross-section of the sample; py equals A times 
the real sound pressure, 

density of frame or air; 02:2=ho,, 

00 density of free air, 
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effective density of the air in the material; 
see eq. (14), 

o density of air in cylindrical pores, according 
to CranpALL; see eq. (9), 


S=s/jo, 

s coupling coefficient as defined in eqs. (1) 
and (2), 

Oo specific air resistance of the material, 

t time, 

U45V5 mean velocity of the skeleton and of the air: 
V2=1/h times greater than the velocity of 
volume displacement, 

W __ wave or characteristic impedance. 
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Bookreview 


Surveys in Mechanics; The G.I. Taytor 70th 
Anniversary Volume. Edited by G. K. Barcurtor 
and R.M. Davies. Cambridge University Press, 
Cambridge 1956, 2.10 s. Od. 


The name of Sir Grorrrey Taytor means the same 
to workers in fluid mechanics and acoustic today that 
that of Lord Rayteicu meant to an earlier generation. 
Both produced a continued series of works in these 
fields which will remain as sources of inspiration for 
many years to come, and both had the facility of de- 
signing simple experiments capable of giving precise 
results to illustrate theory. This volume is issued in 
celebration of G.I.’s seventieth birthday. The tribute 
to Rayteicu consisted of a six-volume edition of his 
scientific papers. An equally voluminous re-issue of 


G.I.’s papers might have been made, but these are 
fairly easily got at, seeing that most of them are in the 
Proceedings of the Royal Society and the editors had 
the better idea of inviting prominent workers to survey 
the present position in the fields of mechanics in which 
he had worked. 

This forms a stimulating compendium, too long to 
list in detail, but the two papers which touch acoustics 
most nearly are that on stress waves in solids by R. M. 
Davies and that on viscosity effects in sound waves of 
finite amplitude by M. J. Licurumx. A biographical note 
describing how Sir Grorrrey first became interested in 
the fields he has made his own should be read by every 
young scientist starting research. 


E. G. Ricnarpson 


ACOUSTIC (0,1) WAVE TRANSMISSION 
THROUGH THIN RECTANGULAR PLATES 


by E. G. Ercutrr * and R. F. Lampert 


Department of Electrical Engineering, University of Minnesota, Minneapolis 14, 
Minnesota, USA. 


Summary 


Transmission loss measurements of thin, rectangular aluminium plates employing the 
(0,1) acoustic mode of a waveguide are described. Theoretical transmission loss formulas 
are derived for the limiting cases of small and large ratios of frequency to plate thickness. 
In the case of large ratios, a “mass law” formulation is obtained which exhibits an angular 
dependence characteristic of the (0,1) mode. Complete sound transmission predicted for 
angles near grazing is verified by the experimental results for all plates studied. Some 
discrepancies between theory and experiment for small ratios are discussed and an attempt 
is made to correct the calculations empirically. The discrepancies between theory and 
experiment are here attributed to internal plate dissipation and imperfect boundary clamp- 
ing. 


Sommaire 


On décrit des mesures de la perte par transmission, dans des plaques rectangulaires 
minces en aluminium, du mode acoustique (0,1) d’un guide d’onde. On établit des formu- 
les théoriques de cette perte dans les cas extrémes d’une faible et d’une grande valeur du 
rapport entre la fréquence et l’épaisseur de la plaque. On obtient pour les grandes valeurs 
de ce rapport une «loi de masse» présentant une influence de ]’angle qui est caractéristi- 
que du mode (0,1). Pour toutes les plaques étudiées, les résultats expérimentaux ont con- 
firmé la transmission acoustique totale prévue pour des incidences presque rasantes. On 
examine quelques divergences constatées entre la théorie et l’expérience pour des faibles 
valeurs du rapport ci-dessus, et on essaie de corriger empiriquement les calculs. On attri- 
bue ces divergences a la dissipation 4 l’intérieur de la plaque et a un serrage imparfait 
aux bords. 


Zusammenfassung 


Es werden Dammungsmessungen beschrieben, die an diinnen rechteckigen Aluminium- 
Platten unter Benutzung der (0,1) Mode des akustischen Wellenleiters ausgefiihrt wurden. 
Fiir die Grenzfalle kleiner und grofer Werte des Verhialtnisses Frequenz zu Plattendicke 
werden Formeln abgeleitet. Bei groBen Werten dieses Verhiltnisses erhalt man ein ,,Masse- 
Gesetz“, das eine fiir die (0,1) Mode charakteristische Winkelabhangigkeit aufweist. Der 
Totaldurchgang, der bei nahezu streifendem Einfall theoretisch zu erwarten ist, wurde bei 
allen untersuchten Platten experimentell bestatigt. Einige Unterschiede zwischen Theorie 
und Experiment, die bei kleinen Werten des Verhiltnisses auftreten, werden diskutiert und 
die Formeln versuchsweise empirisch korrigiert. Die Unterschiede werden auf innere Ver- 
luste und unyollkommene Randeinspannung der Platte zuriickgefiihrt. 


1. Introduction 


Sound isolation by panels is a highly important 
aspect of architectural acoustics. Problems involv- 
ing sound transmission through panels and _ parti- 
tions are generally quite complicated and the cur- 
rent state of knowledge leaves much room for ad- 
justment between theory and experiment. Among 
the significant factors which are difficult to evaluate 
even experimentally are the finite size of the plate, 


* boundary and internal dissipation, plate thickness, 


=> 


* This work has been taken in part from a thesis 
submitted in partial fulfillment of the requirements for 
the degree of Master of Science at the University of 
Minnesota, 1956. 

Present address: Technische Hochschule, Institut fiir 
Schwachstromtechnik, Wien, Austria. 


and the realization of a particular type of incident 
wave. In the following we will concern ourselves 
with thin plates. By thin plates we mean plates 
whose thickness A is much smaller than any mecha- 
nical wavelength in the material, so that opposite 
surfaces move in phase. 

Within an acoustic waveguide sound waves at 
normal and oblique incidence can be realized very 
closely. For normal impedance measurements [1] 
the dominant (0,0) mode or plane wave is usually 
employed. For oblique incidence measurements [2], 
[3], [4] the higher order (0,1) mode may be em- 
ployed because it can be conceived of as being 
composed of two plane waves making an angle 
® =arc sin(f,/f) with the principal axis of the 
guide. Hence, higher order mode techniques should 
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also yield transmission values comparable with 
those obtained by oblique incident plane waves 
striking an infinite plane surface at an angle of in- 
cidence © with respect to the normal. Such an 


acoustic system is much simpler than those found’ 


in building structures but it was hoped that it could 
be treated somewhat more rigorously since the angle 
of incidence is under control. The problems present- 
ed by plate boundary impedance and dissipation 
are still inherent, however. 

Complete solutions of the elasticity equations for 
vibrations of thin circular and supported, rectangu- 
lar plates have been known since the investigations 
[5] of Kircunorr and Rayierex. However, when 
treating finite, rectangular plates having general 
boundary conditions one must resort to approximate 
methods of analysis in order to calculate the nor- 
mal mode frequencies. The method of Rayieren and 
Rrrz [6] allows the calculation of expressions for 
the frequency once the form of the corresponding 
mode has been guessed. By assuming plate modes 
as products of modes of beams subjected to the 
same boundary conditions at their ends as the plate, 
Warsurton [7], obtained expressions for the fre- 
quencies of all modes having their nodal lines par- 
allel to the edges. 

From a knowledge of the plate mode configura- 
tion and frequencies it is possible to calculate values 
for the transmission loss (TL) of sound through 
the plate. We consider in this paper only the (0,1) 


wave of the guide which is of interest to us because , 


of the expected dependence of TL upon angle of 
incidence. Measurements on aluminium plates of 
various thickness having both clamped and rubber 
mounted edges are presented. The results are then 
compared with calculated values based upon ap- 
proximate theoretical expressions developed in the 
following section. 


2. Transmission loss theory 


The formal expression for the vibrational pattern 
of thin plates is readily obtained in the form 


1 ~ ~ Wmn 
ay B - - 1 
oh m=0 a fh Orn — 0 


where @ is the angular frequency of excitation, 
the B,,, are the coefficients of wy», in the plate nor- 
mal mode expansion of the pressure acting on the 
plate, and the plate normal mode functions yy, are 
assumed to satisfy the homogeneous equation 


i Eh’ 


12(1 — a2) ie Orne bh} Wmn=0. 


(2) 

The normal mode angular frequencies W,,, de- 
pend in general upon the dimensions of the plate 
and the boundary conditions at the edges of the 
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plate. Calculations of these frequencies from the 
work of Warsurton [7] for a rectangular clamped 
plate of dimensions a = 34.15 mm and b =72.10mm, 
Fig. 1, yield values for the lowest modes as tabu- 
lated below: 


foo = 1.02 x 10? ch, 
for = 1.27 X10? ch, 
fog = 2.47 X 103 ch, 
foy =5.39 x 103 ch. 


(3) 


Fig. 1. Rectangular waveguide with plate sample. 


The international MKS system of units will be 
used for all calculations. For a supported plate of 
the same size the corresponding frequencies turn 
out to be 

foo S002 Cle, 
Lorde Clee. 
fog = 1.88 21105 ch? 
fey 4.05 10" his 


which lie somewhat below those of the clamped 
plate. 

Within a rectangular waveguide having rigid side 
walls and cross dimensions a and 8, Fig. 1, the com- 
plex amplitude of the acoustic excess pressure can 
be written as the sum of all possible propagated 
and attenuated modes. For the incident (left) side 
one writes (5) 


P(x,y,z)= >) > %an\Pnn® > Pie 


m=0 n=0 


(4) 


where Yn, =cos(mx/a) x*cos(nx/b) y are the cha- 
racteristic functions of the guide and the propaga- 
tion constant Ymn=@/¢y(1—fé, mn/f?)*” is real for 
f>fe,mn With the cut-off frequencies / 


fe, mn = 9/2[ (m/a)? + (n/b)?}*?., (6) 
For this study m=0, n=1 which yields for the 0,1 
mode cut-off frequency, f, =c9/2b. 


The complex amplitudes of the pressure and par- 


ticle velocity components are related by 


Kae FY nn Ping ? (7) 
where 
M2 —j eee tan (™ 2) 
agg . > ee 

ng Qo — nw tan (“2 v) 
b 
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_where the negative sign holds for waves travelling 
in the positive z direction. 


Suppose a thin, rectangular plate of thickness h 
is inserted across the guide as shown in Fig. 1 
If there are no waves reflected from the far ends 
of the guide, the acoustic field on the incident 
side consists of the incident 0,1 wave P§), and 
the reflected waves P), while on the transmis- 
sion side there exist only the transmitted wave P(\) 


and those higher order modes produced by scatter- 


ing of the incident wave upon the plate located at 
z=. The pressure acting upon the plate is the 
sum of all these components, i. e. 

P®) 


P(0) = Fi Xor + » S: tn (P ee a (9) 


=0 n=0 


Since the normal velocity is continuous at the 
plate it follows that 


* yl Ns = PO) — Pe) 


mn mn* 


(10) 


By combining eqs. (9) and (10) one calculates 


2» y Amn P a 


m=0 n=0 


P(0) = (11) 
for the total pressure at the plate. 

The transmission loss (TL) of the 0,1 wave will 
be expressed in decibels as 


| Pa) 


Pe (12) 


Ibe sl 20 logio 


Approximate expressions for TL valid for the 
limiting cases of large and small frequency ratios 
f/fo, will be derived *. If the excitation frequency is 


_ near or below the plate frequency fp; , one can con- 


; 


" 


] 


; 


| 
t 


7 


clude from physical reasoning that the motion of 
the plate will consist essentially of only the 0,1 term 
of the series solution, eq. (1). One therefore writes 


La (13) 


For simplicity we choose the normal mode pat- 
tern for the clamped plate in the form 


Wo = Sin” (#2) sin? (2) cos (22) : 


3 * Note added in proof: This also corresponds ap- 
proximately to large and small ratios of plate span to 


(14) 


_ plate thickness, b/h. Since the actual frequency range 


covered in the experiments is quite narrow the experi- 
mental results are plotted as a continuous function of fre- 


quency employing plate thickness as a parameter. Note, 


in the situation covered by Fig. 3b neither the “mass 
law” nor the low frequency TL expressions strictly 


apply, the former being used for rough comparison 


only. 
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The series coefficient By; in the expansion of 
P(0) turns out to be 


Bo = | [ P() yo dd / [ fyz,44 = 


16 


a les Pe POP O) 0 (15) 


and the normal velocity of the plate V=j wy be- 
comes 
16 


_ jor 1) — pes) (3) ene 
i oh(«5,— o*) a (Poi Po + Pos ang at Pe, 


(16) 


where the coefficents P®, have been eliminated by 
means of eq. (10). 

The choice of the mode pattern expressed by 
eq. (14) is preferred to a closer approximation be- 
cause only a few of the integrals of Wo1 %mn taken 
over the plate contribute to eq. (15) (see appen- 
dix). It is also interesting to note that the value of 
the normal mode frequency fy; obtained from the 
Rayrercu—Ritz method employing eq. (14) lies very 
close to the value calculated using the results of 
Warsurton. In fact, using eq. (14) one calculates 


fox = 1.30 x 108 ch (17) 


a value to be compared with the results of eq. (3). 

The normal velocity V, eq. (16), can also be ex- 
panded into a series of characteristic functions Y,,_ 
where each coefficient is proportional to the respec- 


tive pressure coefficient according to eq. (8). For 
the transmitted waves one obtains 
P&) = =@' 0 Vonsl Ymn = 
= eof [V tmn dAlynn | [2,44 . (18) 


The system of equations relating the transmitted 
and incident pressure coefficients may be expressed 


as 
L—e/yo, — &/Yo €/Yox —€/Yo, 
é/Yo3 1—e/o3 = —8/Yos é/Yos 
Ze/yo, —2e/yoy 1—2e/yoy 2e/yo 
—2 8/723 2/72 2efyog 1—2é/Yog 
Por —€/Yo1 
3 PY f €/Yos po 
Po 2 e/Y24 re ee 
PS — 2/795 
where ¢= — j? w? 09/3 0h (a4? — w?). 


From eq. (19) one readily obtains the desired 
transmission ratio 
Pt) l 1 1 2 2 
: a =e Iz if Thy: ae a ey 
FS Be \ya. yo Ya Ves 
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from which the TL, eq. (12), becomes 
| 2 

TL=20logiy 1 Be) ou 1)+ 
| 2 ar 


( 1h. ee ) 
“7 on ican sheeeg tet 
Y038 5 24 23] | 

In the frequency range lying between f, and the 
cut-off of the (0,3) or the (2,1) guide mode (which- 
ever is lower) the third term on the r.h.s. of eq. 
(21) is small compared with the second term and 
will be negative and imaginary, thus slightly in- 
creasing TL below fp; and slightly decreasing it 
above fy, . By dropping the third term from eq. (21) 


one obtains 
3 eh yor he Ms 7 ‘| 
205 or li 


a form which should be valid for frequencies near 
or below fo, . 

As was to be expected this result predicts perfect 
transmission of the sound wave at the resonance 
frequency of the plate. Similar minima should occur 
at the resonance frequencies fos, fo, , and fy; which 
also contribute to the transmission of the (0,1) 
sound wave. These higher plate modes were neglect- 
ed in the above calculation but there is some ex- 
perimental evidence of their contribution. 

Next we consider the limiting case of high fre- 
quencies where {> fim, m,n = N, where N is a 
sufficiently large number. In this case there are a 
large number of modes available to sustain a velo- 
city distribution proportional to the driving pres- 
sure distribution over the plate. Thus, the reaction 
of the plate becomes more and more that of localiz- 
ed mass elements, i.e. a surface of local reaction 
[8], [9]. In this asymptotic mass model the be- 
haviour is described by the relation 


P= oO ohVs 


(21) 


ties 101g.) i (22) 


(23) 


where P(x,y) and V (x,y) are the total pressure 
and velocity which may be developed in a series of 
characteristic functions like those of the waveguide. 
A straightforward analysis predicts the pressures 
ratio for any mode of the guide, 

Prin 


SB) 7 


mn 


-Q h¥m ne 


(24) 
2 Op 


which yields for the TL of the (0,1) sound wave 


TL=10 logy, f + (es ta) 


2 
2 09 ad 


This formula is the “mass law” relation for the 
TL of the 0,1 sound wave. It is interesting to note 
that the low-frequency approximation, eq. (22), ex- 
ceeds the mass law formula by the amount of 9.55 
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decibels at high frequencies. This difference is 
probably due to the omission of plate modes other 
than the (0,1) mode which also contribute to the 
TL at high frequencies. Both formulas exhibit the 
same behaviour near mode cut-off and predict zero 
TL at cut-off. This behaviour can readily be ex- 
plained in physical terms as an effect of the match- 
ing of the guide impedance in the z direction, 
Zo1 = 09 cL1 —f.2/f?] 1”, to the transverse impedance 
of the plate. It is clear that the mass law formula- 
tion is independent of conditions at the boundaries 


‘of the plate. The TL calculated from eqs. (22) and 


(25) will be plotted against frequency for various 
panel thicknesses and compared with measurements 
taken for several clamped and rubber mounted 
plates. 


3. Experimental apparatus and measuring 
technique 


3.1. Experimental apparatus 


The acoustic waveguide used in the experimental 
studies was a rectangular, copper-plated, brass tube 
of commercial 10cm waveguide stock. Its inner 
cross-dimensions were 34.2 mm by 72.1 mm, and 
its wall thickness was 2mm. The complete experi- 
mental tube consisted of three sections jointed to- 
gether by waveguide flanges, Fig. 2. The source sec- 
tion was terminated in a fiberglas wedge about 
1.0m in length. Two high frequency (University 
Model 4401 tweeters) loudspeakers were joined to 
two slits cut into the opposite side walls of the tube. 
The two speaker units were driven 180° out of time 
phase and thus generated a strong (0,1) acoustic 
wave [10] in the tube. 


\ [tf = standing wave tube 


Termination tube 
Plate mounting device 


Fig. 2. Experimental apparatus for measurements em- 
ploying (0,1) wave. 


The length of the middle (standing wave detec- 
tor) section of the waveguide was 1.8m. A slot of 
width 1.6 mm was provided at the centre of the top 
wall for the pick-up probe. The other three sides 
were embedded in asphalt in order to damp the 
wall vibrations. The transverse vibrating reed pick- 
up device [10] was mounted in a carriage which 


4 


= Ate 
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could be moved along the top of the detector sec- 
tion. Its balsa wood reed probe was about 0.5 mm 
thick, 10 mm wide, and protruded about 5 mm into 


the acoustic waveguide. By virtue of its construction 


the pick-up was sensitive to the transverse particle 
velocity associated with the (0,1) sound wave only, 


_ and hence did not respond to effects of even order 


modes, especially to the dominant (0,0) mode and 
to the (0,2) mode which may be present in the fre- 


~ quency range of interest. Over the frequency range 


extending from the cut-off frequency of the (0,1) 
wave at /, =2.4kc/s to the cut-off frequency of the 


(0,3) wave at 7.2 kc/s no other acoustic disturban- 


ces can be excited as progressive waves. Hence, 
when acoustic measurements are taken sufficiently 
far away from the source and from the test plate 
which terminates the detector section, no inter- 
ference from other higher order modes is to be ex- 
pected. Thus, the detected signal should follow the 
pattern of the (0,1) wave only. 

The third or terminating section of the waveguide 
was also slotted to permit measurement of the trans- 
mitted (0,1) wave by the reed pick-up device de- 
scribed above. This section was also terminated in 
a fiberglas wedge in order to prevent reflections 
from its far end. 

A special mounting device which held the test 
panels was inserted between the detector and term- 
inating sections. It consisted of two heavy, rectangu- 
lar steel frames of which the contact surfaces had 
been carefully machined in order to provide as uni- 
form clamping as possible along all edges of the 
plate. The flanges of this mounting were in turn 
isolated from the waveguide sections by one-half 
inch rubber pads having a hardness of about Shore 
50. The purpose of these rubber mounts was to pre- 
vent the transmission of mechanical vibrations be- 
tween tube sections thereby eliminating flanking of 
the test plate. Flanges, bonded rubber isolators, and 
mounting frames were machined carefully to follow 
the inner dimensions of the waveguide. 

The test plates were cut from aluminium foil and 
sheet material. All plates were cut to a size of 51 
by 89 mm. Metal joints were greased with petroleum 
jelly in order to make them air tight. The test plates 
were clamped either directly between the frames or 
mounted softly between rubber washers. The rubber 
washers were 3.4 mm in thickness and had a hard- 
ness of about Shore 50. 

The electrical driving system consisted of two 
speaker units driven by a GR Model 713A beat- 
frequency oscillator which was calibrated against a 
secondary standard oscillator (WE Model 6010 B) 
and found to be sufficiently stable over the time 


x necessary to complete an experimental run. A small 


strip of sheet iron cemented to the balsa wood 


E.G. EICHLER and R, F. LAMBERT: WAVE TRANSMISSION 


383 


reed of the pick-up device served as the armature 
of a telephone receiver magnet. The two coils of the 
detector were connected in series and the generated 
signal was fed to a battery-operated preamplifier 
(Ballantine Model 220) and from these to a selec- 
tive filter (SKL Model 302), the output of which 
was observed on a cathode-ray oscilloscope (Du- 
mont Model 304 H). Acoustic, electric and magnetic 
shielding had to be employed in order to suppress 
interference from various mechanical and electrical 
sources thereby eliminating erroneous signals not 
directly related to the 0,1 acoustic wave in the 
waveguide. 


3.2. Experimental method 


The relation given in eq. (7) states that the cor- 
responding excess pressure and transverse velocity 
coefficients are directly proportional to each other. 
Hence, the TL as expressed by eq. (12) may be 
rewritten in the form 

Voty 


At a distance from the primary and secondary 
sides of the plate, disturbing effects of other higher 
order modes will have been sufficiently damped by 
inherent high attenuation [11] since the frequency 
lies well below their respective cut-off frequencies. 
Transmission measurements made in this region 
may then be used to determine the TL of the (0,1) 
have. The transverse velocity amplitude will exhibit a 
maximum value in the standing wave pattern on the 
incident side a distance dy from the plate which is 
related to the incident travelling wave transverse velo- 
city 


Vy (dw) | =2| Viry ‘cosh(ady+J1")) e~/* (27) 


on the assumption that a? < #?. In the above 1’) 
is the real part of the pressure reflection coefficient 


(2) 
Pox 4 

a 
Poy 


a} ee e 20 otj r;) . 


(28) 


Let the transmitted normal velocity be measured 
at a distance D from the transmission side of the 
plate. We then have 


| Vi? (D)| =| Vory|e72. (29) 


Eq. (26) for the TL may now be expressed in terms 
of measurable quantities 


aD+I'y 


Vy (dy) a a 


TL=201 hc 
O£10 2V‘) (D) cosh (a dy + L'9) | 


(30) 


The attenuation constant 2 and the standing wave 
ratio S are sometimes needed as corrections to the 
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measured transmission ratio. The reflection coeffi- 
cient can be related to the standing wave ratio S by 
means of the approximate formula 

wy =ar cothS—ady, (31). 
where d,, is the distance from the plate to a cor- 
responding minimum point in the standing wave 
pattern. If this is substituted into eq. (31) one ob- 
tains 


() 5 
TL = 20 logy, | Vy (du) ee) 6 am) 
2Vy (D) 
S< ee — 2 = = ~~: 
1+exp(—2arcothS) 


x 


: (32) 


which now contains S as a correction to the measur- 
ed transmission ratio. For large standing wave ratios 
the third term in eq. (32) is close to unity and 
this correction may be neglected. 


3.3. Results and Discussion 


The results of the TL measurements according to 
eq. (32) are displayed graphically in Figs. 3 to 5. 
The material of the two thinnest plates was Reynolds 
aluminium foil while that of the thicker samples 
was Alcoa aluminium sheet 25—1/2H. All ma- 
terials could be considered to represent relatively 
pure aluminium and had probably been rolled in a 
cold state. 

Upon the charts theoretical curves for TL values 
according to eq. (25), or to eq. (22) for thicker pla- 
tes, are also presented for comparison purposes. The 
physical constants used in the evaluation are as fol- 
lows: Poisson’s ratio 6 =0.3, dilatational velocity of 
sound c= 5200 m/s, density @ = 2.7 x 10° kg/m? for 
aluminium, cy = 348 m/s and the density @) = 1.29 kg 
per m® for air. Temperature and humidity of air 
were measured occasionally and found to average 
about 28° C and 50%, respectively. 

There is general agreement between experiments 
and theory in the general behaviour of transmission 
loss for thin samples. For all samples there is agree- 
ment in the order of magnitude of the loss and the 
fact that TL approaches zero as the frequency ap- 
proaches the cut-off frequency f, from above. Fur- 
ther there is general agreement in that TL increases 
as the log of f?(1—f,?/f?) for values of f/fy; which 
are sufficiently large. There exist TL minima for 
thicker plates which can be attributed more or less 
satisfactorily to plate resonances. It is hoped that the 
direct TL measurements, according to eq. (32), 
were not in error by more than 2 dB. Data could 
be reproduced within 1 dB or less for TL values 
below 20 dB. This result was true only if the mount- 
ings were not altered in any manner. 

In Fig. 5 are presented curves of the lower nor- 
mal mode frequencies of a clamped, rectangular 
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plate as a function of plate thickness according to — 


eq. (3). At these points one should reasonably ex- 


Fig. 3. Transmission loss characteristic of (a) 0.05 mm, 
(b) 0.20mm, (c) 0.63 mm, (d) 1.01 mm, and 


(e) 130mm (—o— experimental, theo- 
retical, and — — — theoretical, corrected empi- 


rically). 


Fig. 4. Transmission loss characteristic of (a) 0.05 mm, 
(b) 0.20 mm, and (c) 1.30 mm rubber mounted 
aluminium plate (—o— experimental, ——— 
mass law). 


Sees 
ean 
Ae 


42mm 


0 02 04 Sa 08 ~~ 1.0 


_Fig.5. Plate normal mode frequency characteristic. 


pect the TL to exhibit low values. The experimental 


_ points plotted in Fig. 5 coincide with sharp TL mi- 


j 


-nima taken from Figs. 3 and 4. For the 0.20 mm 


plate, the thinnest which could exhibit a resonance 
above f,, the agreement with theory is quite good 


if one attributes the findings to the two steeper lines. 


For the 0.63 mm sample the measured points fall 
below the calculated ones. This result is to be ex- 
pected especially if the clamping is not perfect. For 


_ the two thicker plates the observed minima fall con- 


siderably below any possible calculated values. No 


_ rational explanation has been found for this result. 


It is interesting to observe from Fig. 3c that the 


_ agreement between the values predicted by eq. (22) 


E.G. EICHLER and R. F, LAMBERT: WAVE TRANSMISSION 


385 


and the experiments on the 0.63 mm plate can be 
made quite good if the resonance frequency fo, of 
the plate is taken to coincide with the observed TL 
minima occuring at /9,. The results of empirical 
corrections applied to calculations for the 1.01 mm, 
and the 1.30 mm plates are shown in Figs. 3d and 
3e. The fact that the calculated values are in gene- 
ral larger than the measured ones is attributed to 
an effect of neglecting higher order plate modes 
which also contribute to (0,1) wave transmission in 
the derivation of eq. (22). It may also be due to 
the neglect of dissipation in the plate system. 
Perfect transmission at plate resonance as pre- 
dicted by neglecting dissipation according to eq. 
(22) was not found. Substituting a complex fre- 
quency Wnn=Ompn(1+j%7) into eq. (22), permits 


‘one to evaluate the dissipation constant 7 from a 


knowledge of the TL at a minimum point. Consider, 
for example, the 1.01 mm plate, Fig. 3d, which ex- 
hibits a TL of 15.5 dB at fo; =4.8kce/s. One cal- 
culates from these data that 7 = 0.029, a value lying 
near the upper limit of rare, scattered data to be 
found in the literature [12]. Of course, the mount- 
ing undoubtedly contributes to the damping and 
the plate dissipation constant actually would have 
a somewhat lower value. 


4. Conclusions 


Transmission loss (TL) of higher order (0,1) 
acoustic waves through thin aluminium plates has 
been studied over the frequency range 2400 c/s, the 
cut-off frequency of the (0,1) mode through 7000 
c/s. These measurements deal with clamped and 
rubber mounted rectangular plates. TL data employ- 
ing direct transmission measurements were recorded 
and compared with calculated values over a range 
of panel thickness. The measurements followed a 
mass law approximation where agreement was 
reasonably expected, especially for thinner plates. 
At mode cut-off the TL approached zero resulting 
in almost perfect transmission. For thicker plates, 
inherent resonances caused wide fluctuation about 
the mass law curve. These fluctuations were less 
pronounced and occured at lower frequencies for 
the rubber-mounted samples. 

Comparisons of observed TL minima with plate 
resonance frequencies calculated using a RayLercH— 
Ritz approximation or the results of Warsurron 
show good agreement for the thinner plates. For 
the thicker plates, however, the observed frequen- 
cies fall well below calculated ones. This discrepancy 
may be due to imperfections in clamping and to 
dissipation within the plate or at its boundaries. 

Theoretical TL formulas for (0,1) acoustic wave 
transmission valid for large ratios of frequency to 
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plate thickness leading to a “mass law” relation, 
and for small ratios taking into account only the 
(0,1) mode of the plate, were derived. The agree- 
ment between theory and experiment for thicker 


plates leaves something to be desired, particularly’ 


because of the observed discrepancies between meas- 
ured and calculated plate resonance frequencies. 
However, if experimentally observed frequencies are 
used to correct empirically theoretical formulas, the 
comparison is brought into much closer agreement. 
This result brings some confidence to the accuracy 
of the calculations. An empirical dissipation cons- 
tant for the clamped plate material may be calculat- 
ed from TL data by employing a complex Hooke’s 
or Young’s modulus. The observed TL minima are 
here attributed to resonances between standing wa- 
ves, Some attention was given to “coincidence” ef- 
fects [9] which also could give rise to low TL, but 
it may be shown that these effects were not to be 
expected from the experimental arrangement used 
in this study. (Received October 29th, 1956.) 


Appendix 


The following integrals taken over the plate area 
A are needed for calculating the series coefficients 
in the derivation of eq. (15): 


| [ ¥oudA =3.ab/128, 
[ [144 =a8/2, If ¥o1 20144 =ab/16, 
| [ Hdd =a8/2, I [ Vor 20344 = —ab/16, 
[[ @idd=ad/4, I [ ¥o1%1d4 = —ab/16, 
| [ G44 =0b/4, I | Vor %2344 =a 8/16, 
and all other [ [ 91 %mnd4=0. 


Glossary of symbols 


a,b cross dimensions of rectangular plate or wave- 
guide 

a attenuation constant of 0,1 wave 

B phase constant of 0,1 wave 

Bun plate mode coefficients in the expansion of ex- 
cess pressure 

c velocity of dilatational sound waves in plate 
material 

Co velocity of sound in air 

dm,dy distance of transverse particle velocity minima 
and maxima respectively from the plate 
Young’s modulus 

n plate material dissipation coefficient 

f frequency of excitation 
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fe cut-off frequency of acoustic 0,1 mode of wave- 
guide 
finn normal mode frequencies of the rectangular 


plate or waveguide 
Ymn waveguide propagation constants 
Yo1 =f —j a propagation constant of 0,1 wave 


['=I)+jTI, reflection coefficient 


PD angle of incidence with respect to principal 
axis of waveguide 

h plate thickness 

Y%mn . waveguide characteristic functions 

m,n positive integers including zero 

Ne acoustic excess pressure amplitude especially 


upon plate 

1) pl2 3 : : 
P. es 1% ae Ie Ae) coefficients of the waveguide charac- 
teristic functions in the expansions of incident, 
reflected, and transmitted wave pressure, res- 


pectively. 
yp steady state normal displacement 
Wmn normal modes of clamped, rectangular plate 
a) density of plate material 
Qo density of air 
Ui, transmission loss in decibels especially of 0,1 
wave 
0 Poisson’s ratio 
transverse particle velocity standing wave ratio 
V normal velocity amplitude of the plate 
Vinn amplitude of particle velocity 
o angular frequency of excitation 
mn normal mode angular frequencies of the plate 
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ZUR MESSUNG DER DIFFUSITAT VON SCHALLFELDERN 


DURCH KORRELATION 


von P. DAmmic 
III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Theoretische Untersuchungen der Korrelationskoeffizienten fiir Schallfelder in sehr groBen 
Raumen und solchen, deren Abmessungen mit dem Abstand der bei der Messung verwen- 
deten Mikrophone voneinander vergleichbar werden, ergeben fiir inkoharent und unikoha- 
rent am Aufpunkt eintreffende Signale, da die Korrelierten die Diffusititsverhaltnisse im 
Raum nicht immer eindeutig beschreiben. Sie erweisen sich in vielen Fallen als notwendige 
Bedingungen fiir eine bestimmte Diffusitat, jedoch nicht als hinreichend. Die Empfindlich- 
keit des MeByerfahrens kann je nach den zufallig gegebenen Versuchsbedingungen schwan- 
ken, insbesondere ist sie im Bereich héherer Diffusititen gering. 

In Modellversuchen an einem im reflexionsfreien Raum aufgebauten Kugeloktanten aus 
Lautsprechern werden in einigen iibersichtlichen Spezialfallen diese Verhaltnisse untersucht. 
Theorie und Experiment zeigen gute Ubereinstimmung. 


Summary 


Theoretical studies of the correlation coefficients for sound fields in very large halls, 
and in halls of dimensions comparable with the separation of the microphones used for 
the measurement, show that such correlations do not uniquely describe the diffusive proper- 
ties of such rooms, and this whether the signals arriving at the pick-up points are coherent 
or incoherent. 

They are often characteristic of a given diffusivity but not completely so. The sensitivity 
of the technique of measurement becomes small especially when the diffusion is consider- 
able. These facts are illustrated for certain special cases in a model experiment with arrays 
of loud-speakers, forming an octant of a sphere. There is good agreement between theory 
and experiment. 


Sommaire 


Des études théoriques sur les coefficients de corrélation relatifs au champ sonore dans 
de trés grandes salles et dans des salles dont les dimensions sont comparables a la distance 
des microphones utilisés pour les mesures montrent que, pour des signaux arrivant de fagon 
incohérente et unicohérente au point de réception, les coefficients de corrélation ne dé- 
crivent pas toujours d’une facon univoque le caractére de la diffusion dans la salle. Ils 
constituent dans beaucoup de cas des conditions nécessaires mais non suffisantes d’une 
diffusion déterminée. La sensibilité du procédé de mesure peut varier suivant les conditions 
fortuites de l’essai; elle est faible, en particulier, dans la gamme des grandes diffusions. 
On a étudié ces effets dans quelques cas spéciaux intéressants au moyen d’essais sur ma- 
quette effectués sur un octant sphérique constitué de haut-parleurs et installé dans une 
salle exempte de réflexions. L’accord est satisfaisant entre les résultats théoriques et ex- 
périmentaux. 


Einleitung 


Es ist bekannt, daB die Nachhallzeit allein offen- 
sichtlich nicht geniigt, um die akustischen Figen- 
schaften eines Raumes vollstindig zu beschreiben. 
Daher versucht man seit langem, fiir die Raum- 
qualitaten weitere Kriterien zu finden. Diese Ent- 
wicklung ist keineswegs abgeschlossen [1]. Neben 
der Messung der Frequenzkurvenschwankung, Un- 


‘tersuchung der Frequenzabhingigkeit der Nachhall- 
 kurven, Analyse der Schallriickwiirfe von den Wan- 


~ den nach verschiedenen Gesichtspunkten hat sich die 


Bestimmung der Richtungsdiffusitét als besonders 
brauchbar erwiesen. Diese Messung erfolgt mit einem 
Richtempfanger; man bestimmt damit im Grunde 
die mittlere relative Streuung der aus den verschie- 
denen Raumrichtungen einfallenden Energie um den 
Energie-Mittelwert. Cook und seine Mitarbeiter [2], 
[3] sowie Gersuman [4] haben auf ein anderes 
Verfahren zur Diffusitaétsmessung hingewiesen. Sie 
korrelierten die Empfangssignale zweier Mikrophone 
im Raum und vermuteten, daf die in Abhangigkeit 
vom Mikrophonabstand gemessene Korrelierte wei- 
tergehende Aussagen tiber die Schallfeldqualitat ge- 
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stattete als Nachhallmessungen. Man hat in der Folge- 
zeit verschiedentlich versucht, dieses Verfahren an- 
zuwenden. Es unterscheidet sich von der Messung 
der Richtungsdiffusitat wesentlich dadurch, daf} man 
Energiewerte nicht an einem Punkt im Raum, dem 
Aufpunkt, mift, sondern vielmehr eine Energiebezie- 
hung zwischen zwei benachbarten Punkten herstellt. 
Das erste MeBverfahren ist zunachst richtungsselektiv, 
man bildet erst nachtraglich Mittelwerte aus den 
Messungen. Bei der Korrelationsmessung dagegen 
summiert jedes der beiden Mikrophone von vorn- 
herein uber den ganzen Raum, und es kommt eine 
Richtungsbevorzugung erst in zweiter Linie durch 
die Orientierung der Verbindungsachse der Mikro- 
phone, im folgenden kurz Empfangerachse genannt, 
im Raum zustande. Weil diese Mittelungen nicht 
stets vertauschbar sind, ist von vornherein nicht zu 
erwarten, da} die beiden Verfahren aquivalent wa- 
ren. Tatsachlich scheint gerade die notwendigerweise 
gerichtete Aufstellung der Empfanger im Raum bei 
der Korrelation eine Begrenzung der Anwendbarkeit 
des Verfahrens zu bewirken. 

Die vorliegende Arbeit versucht theoretisch und 
durch experimentelle Uberpriifung der Theorie in 
Modellversuchen auf diese Grenzen hinzuweisen. Die 
Verhaltnisse in Raumen sehr hoher Diffusitat sind 
an anderer Stelle studiert worden [3]. Aus experi- 
mentellen Griinden werden hier besonders Felder 
geringer Diffusitat untersucht. Es erwies sich sehr 
bald als zu schwierig, dabei in wirklichen Raumen 
zu arbeiten, weil dort der Aufbau des Schallfeldes 
aus den einzelnen Riickwiirfen von den Wanden zu 
unubersichtlich ist, als daf er zur Kontrolle einer 
Theorie hatte dienen konnen. Deshalb wurde im 
reflexionsfreien Raum des Instituts ein Raummodell 
aus 32 Lautsprechern aufgebaut, die auf einem Teil 
einer Kugelschale angeordnet waren. Sie ersetzten 
gewissermaBen die Spiegelschallquellen eines wirk- 
lichen Raumes. wahrend der dort noch vorhandene, 
koharent auf die Empfanger wirkende Direktschall 
von der Schallquelle entfiel [5]. Die Untersuchungen 
muBSten. sich natiirlich auf wenige tibersichtliche Son- 
derfalle von Schallfeldern beschranken. Wahrend in 
wirklichen Réiumen im allgemeinen weder die Am- 
plituden noch die Phasen der aus den verschiedenen 
Richtungen am Aufpunkt eintreffenden Wellenziige 
gleich sind, wurden hier, um experimentell definierte 
Verhaltnisse zu erhalten, die Schallquellen auf glei- 
che Amplituden eingeregelt. Verschieden groke Dif- 
fusitét wurde dann nur durch Anderung der Raum- 
winkelbereiche hergestellt, in denen die Lautsprecher 
betrieben wurden. Auch die getroffene Auswahl die- 
ser Winkelbereiche selbst stellt Sonderfalle dar. Be- 
ziiglich der Nullphasenbeziehungen der Schallquellen 
untereinander wurden die beiden Grenzfalle vélliger 
(statistischer) Phasenverschiedenheit und _ vélliger 
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Phasengleichheit am Aufpunkt untersucht. Die Ver- 
haltnisse in wirklichen Réumen lassen sich dann in 
die dadurch gegebenen Grenzen einordnen. 


1. Theorie 


Die Grundlagen der Theorie von den Korrelations- 
funktionen sind in der Literatur oft behandelt wor- 
den, z.B. [6], [7], [8], [9], [10]. Daher sollen 
an dieser Stelle lediglich einige Bezeichnungen und 
Tatsachen wiederholt werden, die zum Verstandnis 
erforderlich sind. 

Die Korrelierte zweier Zeitfunktionen w;(¢) und 
u;(t) ist das Faltungsintegral 


vif 
Be fuco ux (t£7) de. 


Dabei bedeutet t eine zeitliche Verschiebung der 
beiden Signale gegeneinander. Ist i=k, so spricht 
man yon Autokorrelation, fiir 7 == k von Kreuzkorre- 
lation. Praktisch mi®t man sogenannte Kurzzeitkor- 
relierte, d.h. die Integration wird nur innerhalb 
endlicher Zeitintervalle ausgefiihrt. Fir t=0 gibt 
das Integral eine mittlere Energie pro Zeiteinheit an, 
die den Signalen wu; und uz gemeinsam zuzuordnen 
ist. Die hier durchgefiihrten Messungen sind Kor- 
relationen von Schalldrucken an zwei bestimmten 
Stellen im Raum. Daher tritt neben der Zeitverschie- 
bung auch der Ortsvektor % als Parameter auf, bzw. 
es steht statt des Ortsvektors selbst der Abstand r 
der Mikrophone voneinander. Ferner ist das Integral 
hier an der Stelle t=O genommen. Normiert man 
dann noch wie tiblich mit den Effektivwerten von u; 
und u;, so erhalt man den ,,Korrelationskoeffizien- 
ten“. Fiir zwei Punkte P (Ortsvektor #,) und Q 
(Ortsvektor 2) soll festgesetzt werden: 


Ry (|e2—21|) = Rj (r) = 
iB 
Were Bs = 
T Ec © 1) Us (t, #9) dt 


a 0 


j : = T = Et Bat Lo at T = oe 
| 7 eac x4) dt Va fare Zo) dt 
0 0 


(ty Up) 


V (uy?) (u2”) 
Dann ist Rys(0) =1, und R schwankt zwischen +1 
und — 1. Bei R= +1 spricht man von Unikoharenz, 
bei R=O von Orthogonalitat der korrelierten Si- 
gnale. (u,u,) bedeutet das Zeitmittel des Produkts 
U; Uy. In den folgenden Abschnitten werden nun die 
theoretischen Korrelierten fiir die Schalldrucke an 
verschiedenen Stellen des untersuchten Raumes er- 
mittelt. AuBer der Beriicksichtigung der Kohirenz 
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oder Inkoharenz der am Aufpunkt eintreffenden 
Signale wird dabei unterschieden zwischen unendlich 
ausgedehnten oder wenigstens sehr groBen Raumen 
und solchen, deren Dimensionen mit dem Abstand 
der den Schalldruck empfangenden Mikrophone ver- 
gleichbar sind. 


1.1. Der unendlich ausgedehnte Raum 


Da die Raumdimensionen sehr grof} gegen den 
Empfangerabstand sein sollen, wird die Betrachtung 
fiir ebene Wellen durchgefiihrt. Diffuse Feldvertei- 
lung liegt dann vor, wenn in jedem Punkt des Rau- 
mes bei stationaérer Anregung in hinreichend langen 
Zeitabschnitten Wellen gleichmafig aus allen Raum- 
‘richtungen eintreffen, deren Amplituden gleiche Ef- 
fektivwerte haben und deren Nullphasen in jeder 
Finfallsrichtung ebenfalls gleichmafig statistisch alle 
moglichen Werte zwischen 0 und 2a annehmen (In- 
koharenz). Insbesondere wird dann auch die mittlere 
Energiedichte in jedem Raumpunkt gleich sein. 

Zur Berechnung des Korrelationskoeffizienten in 
der Umgebung eines Aufpunktes P wird eine Kugel- 
schale vom Radius 0) um diesen Punkt betrachtet. 
Alle Flachenelemente dQ sollen gleichmafig je einen 
Strahler mit der Schnelle v(t; 2) tragen, derart, dab 
dadurch im Aufpunkt P die Erregung 

A-cos(M@t+@o) dQ 

verursacht wird. Die Nullphase nehme alle méglichen 

Werte von 0 bis 22 statistisch an, jedoch gleich- 

mafig in hinreichend langer Zeit. Ein fester Null- 

phasenwinkel, der den Abstand 0 der Strahler vom 
_ Aufpunkt P beriicksichtigt, ist zur Vereinfachung der 

Schreibweise fortgelassen. Die gesamte Erregung im 

Kugelmittelpunkt P ist dann 


up = A [cos( t+ Vo) dQ, 
wobei das Integral tiber die ganze Kugelschale er- 
‘streckt wird. In einem benachbarten Punkt Q ergibt 
sich entsprechend ug = A f cos(W@t+@ag—k Ag) dQ. 
Dabei ist k= w/c die Wellenzahl der betrachteten 


Schwingung, 49 der sich gegeniiber P ergebende 
Gangunterschied. Es gilt also 


cosk Ag dQ 
‘ } (1) 


Abb. 1. Unendlich ausgedehnter Raum. Eintreffen ebe- 
ner Wellen an den Empfaingern P und Q. 
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den Wellenziigen von den einzelnen Strahlern zu- 
geordnet sind, denn es iiberlagern sich bekanntlich 
die Autokorrelierten inkoharenter Anteile eines Si- 
gnals additiv, wahrend die Kreuzkorrelierten ver- 
schwinden. Da der Kugelradius grof gegen den 


Empfangerabstand sein soll, 0) >r, wird mit 
Ag =r-cos Vg (Abb. 1) 
[ | cos(ker cos?) sin } dd dp 
Rpg = —— . 


[ [sind ad dp 
Dieser Ausdruck Jat sich ftir den Fall, daB die 


Variablen # und  voneinander unabhangig sind, 
also in festen Grenzen laufen, umformen in 


Rpg= (2) 
ms (cos 0; + cos V5) | sin Ei 


cos > (cos 0, —cos Ds) 


at (cos 0, — cos By) 


fiir Pe? SO. 
Diese Funktion ist vom Azimut ~ unabhangig. Fur 


den ganzen Raum gilt insbesondere mit 0 <# <2 


(Abb. 2, Kurve a) [3] 


Ree ee (3) 
kr 


Abb. 2. Theoretische Korrelierte fiir vollstindig diffuses 
Schallfeld (a), unikohérente Einstrahlung aus 
dem Gesamtraum (b) und aus einem Kugel- 
oktanten (c). 


Treten die Nullphasen der Elementarstrahler auf 
den Flachenelementen der Kugeloberflache nicht mit 
statistischer Verteilung auf, fordert man vielmehr 
gleichphasig schwingende Strahler, so erhalt man 
den Fall der Unikoharenz. Die Messung der Rich- 
tungsdiffusitat mit dem Richtmikrophon gestattet 
keine Aussage tiber derartige Unterschiede in der 
Nullphasenverteilung, weil man dort fiir jede Rich- 
tung einfach Effektivwerte, z. B. des Schalldrucks, auf- 
nimmt. Die Korrelationsmessung ftihrt jetzt jedoch 
zu anderen Ergebnissen. Fiir den Kugelmittelpunkt P 
darf man im betrachteten Spezialfall ~g =O setzen, 
die Erregung wird dort also jetzt up = A if cos @tdQ, 
im Nachbarpunkt Q entsprechend 


ug=A | cos(wt—k Ag) dQ. 


Pte 


‘ ; ’ 2 eT ee 
ee A et 
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Fiihrt man die zeitlichen Mittelungen aus, so ergibt sich 
coskAgdQ — f dQ ; 
Rpq= Siete foe - = (Rpg) inkon. - E (P,Q) « 


fag \([cosk 4ad2) +(f sink dodQ) 
Rpg entsteht also durch Multiplikation des Korrelationskoeffizienten fiir den inkoharenten. Fall mit einer 
ortsunabhangigen Funktion, die die regelmaBige Struktur des Feldes charakterisiert. Der Nenner ist seinem 
Wesen nach der Effektivwert einer Feldgréfe, z. B. des Schalldrucks, und daher positiv zu nehmen. Rpq 
behilt deshalb auch fiir den koharenten Fall stets das Vorzeichen von (Rpq) jnkon.- Zu einem tibersichtlichen 
Ausdruck gelangt man wieder, falls ® und @ in festen Grenzen laufen. Dann gilt wie oben Ao =r + cos Vo 


fiir Q) >7, und es wird 


sin E (cos 0; — cos 0.) 


Rpg =cos fy (eos 0; + cos ,) = = = (4) 


| sin iF (cos #, — cos 04) 


a Ina Ss" (cos ~ cosy) S(2n+1) x, 


oy, cos T (cos 0, +cos 0,| fiir 


“ 


(5) 


Set (2n+1) x<*" (cos; —cos Hy) S(2n+1) ax. 


Diese Funktion ist wieder vom Azimut g unabhan- 
gig. Vergleicht man Gl. (4) mit dem fiir inkoharente 
Erregung erhaltenen. Ausdruck (2), so sieht man 
wegen |sinz| <2 fiir x20, daB unter den ge- 
machten Voraussetzungen stets 


(Rpa) koh. | ae |(Rpa) inkoh, | (6) 
ausfallt. Fir Einstrahlung aus dem ganzen Raum 
wird mitO <0 ax 

QnnSkrS(2Qn41)a 
(23D) ap seh r 2 Ge) ae 
also sttickweise konstant (Abb. 2, Kurve b). 

Die Ausdriicke (2) und (5) hatten sich als unab- 
hangig vom Azimut ~ erwiesen. Es lassen sich daher 
die Korrelationskoeffizienten fiir Einstrahlung aus 
solchen Teilgebieten des Raumes sofort angeben, 
die dadurch entstehen, dafi der Polarwinkelbereich 
0 <%<nx wie bisher beibehalten wird, @ jedoch 
z. B. von 0 nur bis y,<2a lauft. Aus dem Gebiet 
0OSd) <n, 9, <¢~ SL 2K soll also keine Strahlung 
im Aufpunkt P eintreffen. Fir solche Gebiete stim- 
men die theoretischen Korrelierten mit denen fiir in- 
koharente bzw. kohiarente Einstrahlung aus dem 
Gesamtraum tberein, d. h. 


Rpg= Wyte (7) 


sinkr 
kr 
Diese Bedingungen erweisen sich daher fiir diffuse 
Felder wohl als notwendig, sie sind jedoch nicht 
hinreichend. 
Ein weiteres Beispiel ist der Kugeloktant 
0<0<2, 0S95*. 


Fiir inkoharente Einstrahlung folgt aus Gl. (2) dafiir 


(Rea) inkoh, = > (Rpq) Paihe = tee, 


d. h., es ergibt sich auch hier wieder der gleiche Aus- 
druck wie fiir Einstrahlung aus dem Gesamtraum. 
Fur koharente Einstrahlung liefert Gl. (5) dagegen 
(Abb. 2, Kurve c) (8) 


a PA) eo Q2n-2nSkrS(2n41)a 
Rpg 9. cos=e un 
— 2 (2n+1l)aSkrS2(n+]) a. 


1.2. Raum, dessen Dimensionen mit dem Empfanger- 
abstand vergleichbar werden 


Bisher konnte stets vorausgesetzt werden, da} die 
Normalen der bei den Empfaingern ankommenden 
Wellenziige im Aufpunkt P und im Nachbarpunkt Q 
unter dem gleichen gegen die Empfangerachse ge- 
messenen Polarwinkel @ eintrafen. Bei kleineren 
Raumen ist das nicht mehr ohne weiteres der Fall. 
Weil diese Verhaltnisse bei den Versuchen vorlagen, 
ist es erforderlich, zum Vergleich zwischen Theorie 
und Experiment den theoretischen Verlauf der Kor- 
relationskoeffizienten auch fiir nicht mehr parallele 
Wellennormalen zu kennen. Die Rechnungen des vor- 
hergehenden Abschnitts galten fiir den Fall, da die 
unendlich weit ausgedehnte Kugelschale um den Auf- 
punkt P auf allen Flachenelementen gleichmaBig mit 
Strahlern belegt war. Jetzt soll ein Netz von diskre- 
ten, punktformig gedachten Strahlern auf einer Ku- 
gelschale vom Radius @) um den Aufpunkt P ange- 
ordnet sein, so wie es bei den Experimenten tatsach- 
lich der Fall war. Die Polarachse der Kugel stimme 
zunachst mit der Empfangerachse iiberein, und die 
Strahler mégen auf den Breitenkreisen der Kugel 


FI 
@ 


aa 
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befestigt sein. Bei der Berechnung der Gangunter- 

schiede sind dann die Strahlergruppen gleicher ,,geo- 

graphischer Breite“ gleichberechtigt, man hat m ver- 

schiedene Gruppen zu je mn; (A=1,...,m) gleichen 
lh 

Strahlern, also insgesamt n — >’ nj Strahler. Im An- 


pam 
A=1 


satz der Amplituden ware eigentlich deren Abnahme 
mit 1/0 zu beriicksichtigen. Sie ist jedoch unerheb- 
lich gegeniiber der Auswirkung der Abstandsande- 
rungen auf die Phasen, die hier allein betrachtet 
wird. 

Handelt es sich zunachst um inkoharent schwin- 
gende Strahler, so mége jeder von ihnen im Auf- 
punkt P eine Erregung A-cos(@t+q@i,), im Nach- 
barpunkt Q entsprechend A-cos(m t+ Min —k Ai) 
verursachen. Die Nullphasen ~j, sollen wieder mit 
statistischer Verteilung, jedoch gleichmafig in hin- 
reichend langer Zeit auftreten. Dann ist die gesamte 
Erregung im Aufpunkt P 


m n mY 
Up= >. A -cos(@ t+ Pin) 
A=1 w=1 
und in Q 
m n 2 
ua= > a A-cos(Wt+Qin—k A). 
Il pt 
Man erhiilt also als Korrelationskoeffizienten R PQ 
des endlichen Raumes 


Vi 


m 
>. nicos k A; 
Rpg= *—" Y'ni Rea”. (9) 


n—_ 


Mie igel all destunendlichen Raumes ist auch hier 


Rpa das arithmetische Mittel der zu den einzelnen 
Strahlern gehoérenden Korrelierten. 

Fiir den anderen Grenzfall gleichphasig schwin- 
gender Strahler werden die Nullphasen 7;,,=0 ge- 
setzt. Dann ergibt sich analog 


D = n 

Req= Cpe > jm m 
1S’ S'ning cos k(Ai,— Ax) 
/ — — 


A=1 w=1 


mi Minaliace:-F (P,Q) . 


d.h. wie beim unendlich ausgedehnten Raum wird 
der Korrelationskoeffizient fiir koharente Einstrah- 
lung durch Multiplikation des inkoharenten mit 
einer Funktion gebildet, die vom Ort der Empfanger 
abhangig ist und die Schallfeldverteilung charakteri- 
siert. Die Gangunterschiede 4 als Funktion des 
Polarwinkels und des Empfangerabstands r erge- 
ben sich nach Abb. 3 jetzt zu 


(10) 
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A= ai(l RO 5 er cos 9) — 0) = 
QoS 
wr (cos) +," sin? ? ) (11) 
2 


Damit konnen die Gl. (9) und (10) je nach der 
im Einzelfall vorliegenden Strahleranordnung nume- 
risch ausgewertet werden. 


Abb. 3. Zur Errechnung des Gangunterschieds fiir einen 
Raum mit endlicher Ausdehnung. 


1.3, Begrenzungen der praktischen Verwendbarkeit 
des Mehverfahrens. 


Eindeutigkeit und Empfindlichkeit 


Bei der Berechnung der Korrelationskoeffizienten 
fiir den unendlich ausgedehnten Raum hatte sich 
schon gezeigt, daB z. B. die Bedingung (sin kr)/kr 
als Kriterium fiir vollige Diffusitat nicht hinreichend 
ist. Es lieBen sich Felder konstruieren, die die glei- 
che Funktion ergeben, ohne vollstandig diffus zu 
sein. Im Zusammenhang damit sollen die Verhalt- 
nisse bei inkoharenter Uberlagerung im unendlichen 
Raum noch etwas naher untersucht werden. 

Fiir die praktische Brauchbarkeit des Mefiverfah- 
rens muf} man verlangen, daf fiir eine vorgegebene 
Strahleranordnung auf der Oberflache der unendlich 
groBen Kugel um den Aufpunkt P, d.h. fiir vor- 
gegebene Diffusitat des Feldes in der Umgebung 
von P, das aus der Korrelationsmessung abgeleitete 
Ergebnis unabhangig ist von der zufalligen Orien- 
tierung der Empfangerachse im Raum. Das ist aber 
nicht immer der Fall. 

Bis auf einen Normierungsfaktor entspricht nach 
Gl. (1) die Korrelation dem Integral 


Jn= | cos(k+r9) dQ: = [cos (h reosV) dQ¢ 
G G 
| dQE=C, 
G 
iiber ein Gebiet G der Flache G,, tiber das Strahler 
gleichmaBig verteilt sein sollen. k=k-n ist der Aus- 
breitungsvektor einer in P einfallenden Welle, die 
Indizierung bei d{2¢ deutet an, daB er mit d{2 varia- 
bel ist. ry sei der Verbindungsvektor der Empfanger. 
Wahlt man nun bei unverandertem G einen anderen 
Vektor += 7,, so ist dies gleichbedeutend mit einer 
entsprechenden 


mit 


(12) 


Schwenkung der Empfangerachse 
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oder einer Anderung der Lage der Strahler im 
Raum bei unverdndertem Diffusitatsverhaltnissen. 


Der Vektor “) (Abb. 4) liege in der z-Achse eines 


en 


Abb. 4. 


Abb. 5. Kugelzone im Polar- 
winkelbereich 7p bis ?. 


raumfesten Koordinatendreibeins, die Richtung von 
k sei beztiglich der Koordinatenachsen durch die 
Komponenten 


(sin? cos gy, sin? sing, cos 7) 
und die von 7, entsprechend durch 
(sin O cos ®, sin O sin ®, cos O) 
gegeben. Dann gilt 
cos a=cos# cos @ + sin? sin © cos(y — ®) 


fiir den Winkel a zwischen & und ry und Gl. (12) 
geht mit |7)| = [71| =r tber in 


77 — [ cos(E-r4) dQ%= | cos {kr[cos 9 cos O+ 
G G 


+ sin} sin O cos(p— ®)]} dQ;z. 

Das ist fiir festes kr im allgemeinen == Jy. Fir den 
Fall, daB die Integration tiber die gesamte Kugel um 
den Aufpunkt P erstreckt wird, also im Fall des vol- 
lig diffusen Feldes, hangt J nicht mehr von der Rich- 
tung des Vektors 7 ab. Dann tritt also insbesondere 
der Fall ein, da8 die Korrelation von der Richtung 
der Empfangerachse unabhangig wird. 

Als weiteres Beispiel werde noch die Integration 
iiber koaxiale Kugelzonen gleicher Flache, jedoch va- 
riabler Lage auf der (Einheits-) Kugel betrachtet 
(Abb. 5). Die Flache einer solchen Zone ist 


F=2x(cos %)—cos #) , 


es werde bei festem F und 1%) als Parameter also 
stets cos } = cos 0, — F/2 x gewahlt. Fiir eine derar- 
tige Strahleranordnung ist nach Gl. (2) 


Rpa= c05| cos Dy= i) (s=*22)/ “re (13) 


Es ergibt sich also nicht eine Funktion von F allein, 
sondern auch von #. Je nach Lage der Kugelzone 
nimmt daher Rpg verschiedene Gestalt an. Ist um- 
gekehrt Rpg als Kurve etwa durch eine Korrelations- 
messung gegeben, so kénnte man zuniichst daran 
denken, F aus dem Faktor 
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( ._ krF krF 
sin ster 

41 4 
zu bestimmen. Das ist in den praktisch vorkommen- 
den Fallen, die durch die Werte von #) und F fest- 
gelegt sind, jedoch nicht moglich, weil diese Funk- 
tion nicht wirklich den Charakter einer Einhiillen- 
den der cos-Funktion annimmt und man daher nicht 
z. B. auf die Nullstellen extrapolieren kann. Nur fiir 


Do = werden die Verhaltnisse iibersichtlich. Dann 


ist namlich (14) 

Ate TOR Dep Te (i me a 

Rpg= sin —— cos —— =| sin ‘ 
kr PF 4m~ 42 2n 21 


Diese Form des Integrals unterscheidet sich von Gl. 
(3) nur durch eine Mafistabinderung der Abszisse. 
Es ware in diesem Fall daher ohne weiteres denkbar, 
die GroRe F, die als ein Ma fiir die Diffusitat an- 
gesehen werden kann, z.B. aus der Verschiebung 
der ersten Nullstelle zu bestimmen. Man kann jedoch 
nicht, etwa durch eine Drehung der Empfangerachse 
im Raum, beliebige, selbst koaxiale Anordnungen 
einer gleichgrofen Strahlerflache in diesem Fall trans- 


. . vi . 
formieren; wenn die Kugelzone %)= > nicht ent- 


halt, gilt statt Gl. (14) der allgemeinere Ausdruck 
Gl. (13). 

Es ist nun die Frage, ob man eine von der zu- 
falligen Orientierung und Lage unabhiangige Aus- 
sage tuber die Diffusitat erhalt, indem man Messun- 
gen in allen modglichen Raumrichtungen ausfihrt 
und dann den Verlauf der daraus entstehenden mitt- 
leren Korrelierten betrachtet. Es werde dazu wieder 
das raumfeste Koordinatensystem von Abb. 4 ver- 
wendet. Die Mittelung hat die Form 


Rene 1 D Reade = 
Ax 


[ cos[kr cos( ky r)] dQe 


1 é ‘ 

W~—--—— —— dQr. 

An faaz 
G 


Darin ist [ a2 Gy eine Konstante und 
é 


dQ; =sin 3 dd dp, 


dQ; =sinO dO d@. 
Ferner gilt wieder fiir beliebig gerichtetes r : 
cos(k, Tr) =cos # cos @ + sin # sin © cos(p— P). 
Fiihrt man bei festem & zuerst die Integration nach 
dQ; aus, so wird dieses Integral von der Richtung 


von k unabhangig und ist gleich 4x(sinkr)/kr, es 
bleibt also 


1 snkr,,+ sinkr 
Rpg = —— dts 
si af kr. >|. aanea 
C 


? 


unabhéngig von der Verteilung der einfallenden 
_ Wellen. Man gewinnt also auf diesem Wege keine 
_Aussage iiber die Diffusitatsverhaltnisse. Im iibrigen 
_konnte man, selbst wenn die Mit- 
_telbildung nicht zu einem unab- 
_ hangigen Wert fiihrte, von einem 
 solchen Mefiverfahren gegeniiber 
_ der ganz ahnlichen, aber technisch 
_ sehr viel einfacheren Messung der 
_ Richtungsdiffusitat kaum einen grundsiatzlichen Fort- 
_ schritt erwarten. 

- Aufer durch die Mehrdeutigkeit wird die prak- 
: piss Verwendbarkeit des MeBverfahrens dadurch 
 eingeschrankt, da die Empfindlichkeit nicht iiberall 
“igleich groB ist, daB sie insbesondere in den Berei- 
chen groSerer Diffusitaten, wie sie in der Praxis 
; besonders haufig vorkommen, stark abfallen kann. 
Das soll am Beispiel des Korrelationsintegrals (13) 
 erlautert werden. Als ein Ma8 fiir die Empfindlich- 
_ keit beziiglich Anderungen der Diffusitat, d.h. fiir 
_ die Verschiebung der Gl. (13) entsprechenden Kor- 
 relationskurven bei Anderungen der mit Strahlern 
_ besetzten Flachen, wird die Verschiebung der Null- 
stellen von GI.(13) in Abhangigkeit von J) und F 
betrachtet. Die erste Nullstelle des cos-Anteils liegt 
bei (kr)g= t2/2[cos 9) — (F/4x)]. Fir kleine 
Anderungen von #, und F gilt 


2(kr) ? 
IU 


cos 


Rpq(t) = 


. 


: 
y 


d(kr)y=+ (sin Be diy + ft): 
41 


_ Nimmt man an, da (kr) 9 vorgegeben sei und F zu- 
nachst einen festen Wert habe, so hat man die grofte 


Nullstellenverschiebung fiir Oo=5 zu erwarten. Da- 
fiir wird dann d(kr) y= + (22?/F?) dF. Die Ab- 
solutempfindlichkeit nimmt fiir wachsendes F also 
sogar quadratisch ab. Die gleichen Verhiltnisse lie- 
gen fiir die Nullstellen des sin-Anteils vor; dafiir 
- gilt d(kr) = F (422/F2) dF unabhingig von ¥y. 
Praktisch bedeutet das, da die Korrelationskurven 
- gerade im Bereich groSer Diffusitéten nur wenig 
_ von der theoretischen Kurve des vollstandig diffusen 
Feldes, (sinkr)/kr, abweichen. Darauf haben in 
etwas anderem Zusammenhang auch Cook und seine 
- Mitarbeiter hingewiesen [3]. 


a 4, Ergiinzungen 


Es bleibt noch kurz zu untersuchen, welche Veran- 
derung die bisher errechneten Korrelierten erfahren, 
wenn man bei der Korrelation Laufzeitverschiebungen 
#0 zulaBt: 
a - 
a ania) = ee) walt)? 
Viup?(t) > ug? (t) > 
Statt des nur vom Verbindungsvektor der Punkte P 

und Q abhingigen ,,Korrelationskoeffizienten“ erhilt 
an dann eine orts- und zeitabhangige Korrelations- 
ktion. Fiihrt man die Rechnung fiir den inkoharent 
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erregten unendlichen Raum durch, so ergibt sich statt 
Gl. (2) unter den dort gemachten Voraussetzungen 
(Kugelradius 9) > Empfangerabstand r ; # und — von- 
einander unabhiangig) 


k 


5 (cos 3, + cos Jy) —w : 


sin E (cos }; —cos 73) 


- (15) 
2 (cos 7, — cos My) 


Fiir den ganzen Raum gilt mit 0 < # < =x insbesondere 


Rpq(t) =cos w poo : 
kr 

Gl. (15) unterscheidet sich von Gl. (2) also durch eine 
Phasenverschiebung wt des cos-Faktors im Zahler, die 
sich so auswirkt, daB z. B. die Korrelationsfunktion fiir 
den ganzen Raum aus Gl. (3) einfach durch Multipli- 
kation mit der Autokorrelierten cos w 7 der erregenden 
Schwingung hervorgeht. Die Einfiihrung der Laufzeit- 
verschiebung liefert hier also keine die Schallfeldyertei- 
lung im Raum betreffende neue Information. Im allge- 
meineren Fall der Gl. (15) vermag sie die wesentlichen 
Mangel der Korrelationsmessung, insbesondere die dort 
auftretende Mehrdeutigkeit von Rpaq, z. B. wegen der 
Unabhangigkeit vom Azimut qm, die Abhangigkeit von 
der Orientierung der Empfangerachse im Raum und die 
wechselnde Empfindlichkeit in den verschiedenen Werte- 
bereichen nicht aufzuheben. 

Im Anschlu8 hieran la8t sich auch die Frage beant- 
worten, ob es zu Diffusitatsmessungen etwa ausreichend 
ist, statt der bisher stets betrachteten Kreuzkorrelatio- 
nen zwischen den Empfangssignalen zweier um den Ab- 
stand r voneinander entfernten Mikrophone Autokorre- 
lationsmessungen mit einem einzigen Mikrophon am 
Aufpunkt P durchzufiihren, wobei dann die Laufzeitver- 
schiebung t als Parameter auftritt. Dazu hat man, wenn 
wieder der Fall des inkoharent erregten unendlichen 
Raumes betrachtet wird, in Gl. (15) r—0O gehen zu 
lassen. Man erhalt dann 


Rpp(t) =cos@mt, (16) 


unabhiangig von den bei der Anregung beteiligten Win- 
kelbereichen. Derartige Messungen sind versuchsweise 
mit Rauschbaéndern verschiedenen Frequenzumfanges 
durchgefiihrt worden. Fiir solche Bander geht Gl. (16) 
durch entsprechende Summation tiber die beteiligten 
Frequenzen tiber in 

: ie 

sin Aw 


“ 


Rpp(t) =cos Wm T (17) 


T 
Aw 9 
wenn @m der Mittelfrequenz, 4@ der Breite des ver- 
wendeten Bandes entspricht. Auch darin sind Aussagen 
iiber die Schallfeldverteilung nicht mehr enthalten. 
Autokorrelationsmessungen lassen jedoch unter ge- 
wissen Bedingungen Schliisse zu auf die Laufzeitunter- 
schiede und Amplituden koharenter Einzelvorginge, aus 
denen sich die korrelierten Signale zusammensetzen [11]. 
Betrachtet werde z. B. ein von einer Schallquelle abge- 
strahltes Signal f(t), das einmal direkt und zum zwei- 
ten auf dem Umweg iiber einen Reflektor am Auf- 
punkt P eintrifft. Beide Signale interferieren dort und 
moégen eine Erregung up=f(t)+a/f(t—o) erzeugen, 
wobei o die durch den Umweg iiber den Reflektor her- 
vorgerufene Zeitverzégerung der beiden Signale gegen- 
einander und a <1 eine Konstante bedeutet. Dann 
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wird mit Spp(t) =<f(t)- f(t+7)? (18) 

opts) = (ht2) SeP(s) +2 [Sep (r-+0) + Spr(r—0)} 
(1+a7) Spp(0) +2 a Spp(o) 

Gilt nun yon hinreichend groBem |r| an Spp(t)~0,, 

wie es z. B. fiir weiBes Rauschen oder auch fiir Rausch- 


binder geniigend grofer Bandbreite erfiillt ist, so folgt 
aus Gl. (18) fiir eine kleine Laufzeitverschiebung Az 


Rpp(0+ At) =Spp (0+ At)/Spp (0) , 
Rpp(o + Ar) =a Spp(0+ At) / (1+?) Spp(0) = 


+4 __Rpp (0+ Ar) . 
+a 


; (19) 


Im Korrelogramm tritt dann also zweimal die Auto- 
korrelierte Spp(z) des Signals f(t) auf, und zwar ein- 
mal am Nullpunkt der Zeitachse selbst, zum anderen 
an der durch die Echolaufzeit 0 festgelegten Stelle. Die 
Amplituden verhalten sich dabei wie 
Rpp(0) _ 1+a? (20) 
Rpp (a) lo} 


und daraus folgt Rpp(0)/Rpp(o) 2 2. 


2. Experimenteller Teil 
2.1. Apparatur 


Zur experimentellen Uberpriifung der theoretischen 
Ergebnisse wurden die Verhaltnisse in einem Schall- 
feld dadurch angenahert nachgeahmt, da im re- 
flexionsfreien Raum des Instituts ein Kugeloktant 
aus Lautsprechern aufgebaut wurde. 32 Lautsprecher 
(LP 100/19/70) waren auf den Breitenkreisen eines 
entsprechenden Teils einer Kugeloberflache von 
2,60 m Radius aufgehangt, derart, daB die mittlere 
Belegungsdichte auf der ganzen Flache etwa kon- 
stant war. Der Winkelabstand der Breitenkreise von- 
einander betrug annahernd 15° (Abb. 6). Die MeB- 


frequenz war 2,5 kHz; um die Ausbildung stehender 


Abb. 6. Lautsprecher und Mikrophone im reflexionsfreien Raum. 
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Wellen an der Mikrophonaufhangung zu vermeiden, 
wurde mit einem Rauschband von 200 Hz Band- 
breite um diese Mittelfrequenz gearbeitet. Die Laut- 
sprecher konnten gruppenweise parallel geschaltet 
werden, je nach dem Raumwinkelbereich, fiir den 
die Korrelation gemessen werden sollte. Mit dieser 
Betriebsweise, zum Teil mit RC-Phasenschiebern in 
den Zuleitungen zur Korrektur der Nullphasen, 
wurde der Fall der koharenten Uberlagerung der 
Signale realisiert. Um auch inkoharente Signale zu 
erzeugen, muBten die Nullphasen fiir jeden einzel- 
nen der beteiligten Lautsprecher statistisch variiert 
werden, Das war dadurch zu erreichen, das die vor 
dem reflexionsfreien Raum gelegene groBe Versuchs- 
halle (Volumen etwa 750 m?, Nachhallzeit 2,2 s bei 
2,5 kHz) stationar mit dem 2,5-kHz-Rauschband er- 
regt und 14 Lautsprecher von ebenso vielen in der 
Halle unregelmaSig verteilten Mikrophonen itiber 
getrennte Verstarkerkanale betrieben wurden. Zur 
Anregung der Halle diente dabei ein Dodekaeder 
von Hochtonlautsprechern mit annahernd kugelfor- 
miger Richtcharakteristik (nach Késters und Harz). 
Korrelationsmessungen mit gefiltertem Rauschen lie- 
ferten Anhaltswerte fiir die erforderlichen Mindest- 
abstande der Mikrophone voneinander (Abb. 7). 
Weil dabei der koharent empfangene Direktschall 
storen wurde, war der Lautsprecher hinter einer 
Blende aus Sperrholz aufgestellt. Dadurch lag im 
verwendeten Frequenzbereich der vom Direktschall 
herrtthrende Anteil an der Energiedichte des Schall- 
feldes um etwa 30 dB unter der reinen Nachhall- 
energie. Dies ist auch aus Abb. 7 ersichtlich, und 
ebenso, daB fiir Mikrophonabstande r>40 cm die 
empfangenen Signale als voneinander unabhangig be- 
trachtet werden konnen. Fir 
den eigentlichen Versuch mit 
unabhangig erregten Laut- 
sprechern im reflexionsfreien 
Raum wurden die Mikrophon- 
abstande noch wesentlich gr6- 
Ber gewahlt. 

Zur Bildung der Korrela- 
tionsintegrale 

T 


1 
(uj uz) = pp | wea 


0 


wurden die Eingangssignale 
u; und u,, zweier hochwertiger 
Kondensatormikrophone 
(Telefunken) mit kugelf6r- 
miger Richtcharakteristik zu- 
nachst in einem Ringmodu- 
lator multipliziert [12] und 
dann durch ein Fluxmeter 


(AEG) integriert. Der Pha- 


ae ae 
Ss wy ly 
oat 

t is - 
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: “Abb. 7. Versuchshalle. Korrelogramm fiir Anregung mit 


Rauschen. Terzband bei 2,5 kHz. 
Messung mit Blende, 
a Dicsouns ohne Blende. 


_ senabgleich der beiden getrennten Verstarkerkandle 
des Korrelators geschah mit RC-Phasenschiebern, die 


in dem einen Kanal Phasendrehungen von 0 bis 2 
zulieBen. Laufzeitverschiebungen in den Korrela- 
tionsintegralen, z. B. bei der Realisierung von Gl. 
(17) oder (19), konnten durch ein Magnettonband- 
gerat mit zwei verschiebbaren Hork6pfen, hergestellt 
werden. Fiir kontinuierliche Registrierung der Kor- 
relierten durch einen Pegelschreiber mit linearem 


_ Potentiometer (Briiel & Kjzr) diente dann als Inte- 


grator ein RC-Glied. 


2.2. Diskussion der Mebergebnisse 


Durch die mit der beschriebenen Anordnung aus- 
gefiihrten Versuche werden die theoretischen Uber- 
legungen gut bestatigt, soweit eine unmittelbare 
Ubertragung méglich ist. Bei der Beurteilung der 
MefSgenauigkeit darf nicht vergessen werden, dab 
es sich bei diesen Experimenten um Modellversuche 
handelt, um grobe Anniherungen an die Verhilt- 
nisse im wirklichen Schallfeld, die selbst in dieser 
Form z. T. erheblichen technisehen Aufwand erfor- 
derten. So war es z. B. nicht méglich, die Korrelierte 
fiir den Fall der gleichmaSig mit inkoharenten Strah- 
lern gelegten Achtel-Kugel genau nachzubilden. In 
Abb. 8 sind zum Vergleich die Funktion 


R,=(sinkr)/kr (Kurve a) 


und die Kurve R, eingetragen, die durch Summation 
der Korrelierten fiir die einzelnen Lautsprecher nach 


Gl. (9) und Gl. (11) theoretisch entstehen miiBbte 


Abb. 8. Kugeloktant mit inkohirenten Strahlern. Kor- 
relierte fiir gleichmaSige Belegung mit Strah- 
lern (a) und fiir Belegung mit Strahlernetz (b). 
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(Kurve b). Beide Funktionen sind einander im Ver- 
lauf zwar ahnlich, weichen jedoch immerhin merk- 
lich voneinander ab. Die eingetragenen Mefipunkte 
entsprechen inkoharentem Betrieb der Lautsprecher, 
sie passen sich der Funktion R, gut an. In allen 
weiteren Diagrammen sind die theoretischen Kurven 
fiir inkoharente Anregung nach Gl. (9), ftir koha- 
rente nach Gl. (10) berechnet, jeweils unter Verwen- 
dung von Gl. (11). Als Abszissenmafstab ist der 
mit der Wellenzahl & multiplizierte Mikrophonab- 
stand r gewahlt. Die MeBergebnisse sind als Kreise 
eingetragen, soweit sie sich auf inkoharenten Betrieb 
der Lautsprecher beziehen, als Punkte, soweit sie 
von gleichphasiger Erregung herstammen. 

Fiir die in Abb. 9 dargestellten Korrelationskurven 
wurden die Lautsprecher auf je einem Breitenkreis 
des Kugeloktanten angeregt. Neben den theoretischen 
Kurven erscheinen in den Diagrammen sowohl die 
MeBpunkte fiir inkoharente als auch fir phasen- 
gleiche Anregung. Wie man leicht einsieht, stimmen 
namlich die theoretischen Kurven hier fiir die beiden 
Falle tberein, Die Mefpunkte nahern sich den. theo- 
retischen Werten gut an. Es ist aus diesen Kurven 
bereits ersichtlich, dafi sie im Bereich kleiner Win- 
kel ? fast die gleiche Form haben, wahrend fiir ? 
nahe 90° beim Ubergang von Kurve zu Kurve, z. B. 
in der Lage der Nullstellen, sehr gut sichtbare Unter- 
schiede auftreten, entsprechend der groReren Emp- 
findlichkeit des Mefiverfahrens in diesem Bereich. 
Daran wiirde wegen der erforderlichen Normierung 
auch nichts geandert, wenn man in den Bereichen 


Abb. 9. Korrelierte fiir kohaérente Lautsprecher auf Brei- 
tenkreisen. Polarwinkel #=0° (a), 15° (b), 
30° (c), 45° (d), 60° (e), 75° (f), 90° (g). 
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kleiner Winkel # die Dichte der Strahler auf den 
Breitenkreisen erhoht. Der Unterschied der zu er- 
wartenden. Werte fiir den unendlich ausgedehnten 
und den hier realisierten endlichen Raum wird be- 
sonders gut bei der Kurve fiir # = 90° sichtbar. Fiir 
den unendlichen Raum hitte die Korrelierte iiberall 
den Wert R= +1, die gemessenen Werte weichen 
jedoch in Ubereinstimmung mit der theoretischen 
Kurve fiir den endlichen Raum von etwa kr=8 an 
(entsprechend einem Mikrophonabstand von 18 cm 
bei 2,5 kHz) merklich von +1 ab. 

Einige weitere Beispiele fiir die wechselnde Emp- 
findlichkeit in den verschiedenen Winkelbereichen 
sind in Abb. 10 gegeben. Es sind Korrelierte fir 
Anregung der Lautsprecher in Gebieten entsprechend 
den oben behandelten Kugelzonen. Im oberen Teil 
der Abbildung sind die Korrelationskurven fiir die 
Zonen zwischen ?=0° bis 15°, 15° bis 30° und 
0° bis 30° fiir koharente Erregung dargestellt (Kur- 
ven a, b und c). Die mit Strahlern besetzten Flachen 
verhalten sich in diesen drei Fallen wie 1 : 3: 4, die 
theoretischen Kurven ebenso wie die gemessenen 
Werte unterscheiden sich jedoch nur unmerklich von- 
einander. Immerhin ist besonders an den theore- 
tischen Werten sichtbar, dali das Anwachsen der 
Flache nicht einmal durch die Tendenz der Kurven- 
anderung dargestellt wird; die der groBten Flache 


zi 


SY 


0 2 4 6 8 10 2 14 16 eels 


Abb. 10. Korrelierte fiir kohiarente Lautsprecher auf 
Kugelzonen. Polarwinkel #=0° bis 15° (a), 
0° bis 30° (b), 15° bis 30° (c), 15° bis 45° 
(d), 60° bis 75° (e), 75° bis 90° (f), 60° bis 
90° (g). 
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entsprechende Kurve (b) liegt zwischen den beiden 
anderen. Das entspricht der oben erlauterten Tat- 
sache, daB die Korrelierten aufer von der Grobe 
der mit Strahlern belegten Flachen auch noch yon 
deren Lage auf der Kugel abhangen. 

Im unteren Teil der Abb. 10 sind die Korrelierten 
fiir drei Zonen im Bereich gréBerer Winkel # dar- 
gestellt, namlich fiir #=60° bis 75°, 75° bis 90° 
und 60° bis 90° (Kurven e, f und g). Die entspre- 
chenden mit Strahlern besetzten Flachen verhalten 
sich hier wie 1: 1,1: 2,1, unterscheiden sich also 
nur in geringerem Mafe als die vorhergehenden Bei- 
spiele. Zwischen den zugehorigen Kurven bestehen 
jetzt aber sehr groBe\Unterschiede. Neben der gro- 
Beren Empfindlichkeit in diesem Winkelbereich wird 
hier auch die Abhangigkeit der Kurvenform und der 
MeBpunkte von der Lage der Zonen deutlich. Wieder 
liegt die der gréBten Flache entsprechende Kurve 
zwischen den beiden anderen (g). Dies ist mit an- 
deren Worten ein Beispiel dafiir, daB sich die Grobe 
der Strahlerflache, entsprechend einer gewissen Dif- 
fusitat, nicht eindeutig in einer Korrelationskurve 
abbildet: Zu den Zonen # = 60° bis 75° und 75° bis 
90°, deren Flachen. sich nur wenig unterscheiden, 
gehoren vollig verschiedene Korrelierte. Sie stimmen 
weder in den Nulldurchgangen noch in den Extrem- 
werten auch nur annéhernd iiberein. Noch deutlicher 
wird dies durch Vergleich der Kurven d und f aus. 
Abb. 10. Es sind das die theoretischen Kurven und 
die MeBergebnisse fiir die Winkelbereiche 3 =15° 
bis 45° und #=75° bis 90°. Diesen Winkelberei- 
chen entsprechen Flachen auf der Kugel, die sich 
wie 1 : 1 verhalten. Von einer Ubereinstimmung der 
Korrelierten ist jedoch gar keine Rede. Zur Gestalt 
der Kurven selbst ist zu bemerken, da sich bei f) 
und g) bereits der Ubergang zu Rpg = + cos(kr/2) 
gema} Gl. (8) und Abb. 2 bemerkbar macht. Die in 
den Diagrammen schraffierten Gebiete sind durch 
die Messung nur schwer exakt wiederzugeben, weil 
dort der Nenner von Gl. (10) sehr nahe an Null 
liegt und daher nur sehr schwer gemessen werden 
kann. 

Es wurde gezeigt, daB auch eine Drehung der 
Empfangerachse im Raum bei sonst unverdnder- 
ten Diffusitatsverhaltnissen vollig andere Korrela- 
tionskuryen liefern kann. Dies wird deutlich durch 
Vergleich der Kurven in Abb. 11 und Abb. 9. Kurve f 
der Abb. 9 stellt die Korrelierte fiir Anregung der 
Lautsprecher auf dem Breitenkreis mit = 75° dar. 
La8t man die Anordnung der Lautsprecher unver- 


Abb. 11. 
Korrelierte fiir kohirente 
Lautsprecher auf Breiten- 
kreis #=75°. Empfanger- 
achse normal zur Polar- 

achse der Kugel. 


Kugelzonen. Polarwinkel =0° bis 15° (a), 


oe 12. Korrelierte fiir koharente Lautsprecher auf 
0° bis 30° (b). Empfingerachse normal zur 


Polarachse der Kugel. 


; 
i | 
- indent und dreht nur die Empfangerachse aus der 
- Polarachse der Kugel in die dazu normale Lage, so 
_ ergibt sich der durch MeSwerte markierte Kurven- 
_ verlauf in Abb. 11. Das hat, wie man leicht einsieht, 
~ die Ursache darin, da im ersten Fall alle Laut- 
_ sprecher unter dem Polarwinkel 0=75° auf die 
4 Empfanger wirken. In der zweiten Stellung liegen 
_ jedoch die gegen die gedrehte Empfangerachse ge- 
_messenen Winkel im Bereich zwischen # =15° bis 
_ 0 —90°. Zwei weitere Beispiele ergeben sich durch 
_ Vergleich von Abb. 10 und Abb. 12. In Abb. 10 
waren, Korrelationskurven fiir Kugelzonen im Polar- 
- winkelbereich = 60° bis 90° und 75° bis 90° an- 
" gegeben (Kurven g und f), die einen charakteristi- 
_ schen Verlauf zeigen. Abb. 12 stellt durch MefSwerte 
- gwei dazu sehr abnliche Kurven dar, die dadurch 
_ entstanden sind, da8 nach Drehung der Empfinger- 
_achse in die zur Kugelpolarachse normale Lage die 
- Lautsprecher in den Zonen zwischen den urspriing- 
_ lichen Breitenkreisen 3=0° bis 30° (Kurve b) und 
§=0° bis 15° (Kurve a) erregt wurden, Dadurch 
_ergeben sich auch fiir die neue Lage der Empfanger- 
i Winkelbereiche von 3 —60° bis 90° bzw. 
# =75° bis 90°, woraus sich die groBe Ahnlichkeit 
des Kurvenverlaufs erklart. Da8 sie sich nicht vollig 
sntsprechen, liegt daran, dais die Lautsprechervertei- 
lung innerhalb der Bereiche nicht genau iiberein- 
‘stimmt und daher auch die Uberlagerung gema} 
‘Gl. (10) etwas anders ist. Bemerkenswert ist aber 
die Tatsache, da® trotz der Ahnlichkeit der Kurven 
die den Bereichen entsprechenden Fliichen véllig ver- 
schieden sind. Im ersten Beispiel (Abb. 10 g und 
12 b) verhalten sie sich wie 3,7 : 1, im zweiten Bei- 
spiel (Abb. 10f und 12a) sogar wie 7,6:1. Damit 
zeigt sich wieder sehr deutlich, da8 die Korrelations- 
urven kein eindeutiges Abbild der Feldzusammen- 
zung liefern und dal} insbesondere Riickschliisse 
aus den Korrelationskurven auf die sie erzeugende 
Feldverteilung kaum méglich sind, solange man nicht 
die Orientierung der Empfangerachse zum Feld kennt. 
Gerade dariiber aber ist bei praktischen Messungen 
nichts bekannt. ! 
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Zur Erlauterung der Erscheinungen waren bisher 
stets solche Korrelationskurven verwendet worden, 
die durch koharenten, genauer phasengleichen Be- 
trieb der Lautsprecher entstanden. Ebenso lieSe sich 
alles mit inkoharentem Schall zeigen, wenn also die 
Nullphasen der von den Lautsprechern abgestrahlten 
Signale statistischen Schwankungen unterworfen sind. 
Nach Gl. (9) und (10) unterscheiden sich die Kor- 
relierten fiir koharente und inkoharente Uberlage- 
rung der Signale von den einzelnen Schallquellen 
durch eine vom Abstand der Empfanger abhiingige 
Funktion. Sie bewirkt eine Anderung des Kurven- 
verlaufs, la8t aber die Vorzeichen und die Nullstellen 
ungeandert. Dies pragt sich aus in den in Abb. 13 
wiedergegebenen Kurven. Im oberen Teil sind die 


kre ——~ 


Abb. 13. Korrelierte fiir inkoharenten Betriek der Laut- 


sprecher. Lautsprecher auf Kugelzonen. Polar- 
winkel #=45° bis 75° (a), 60° bis 75° (b), 
60° bis 90° (c), 75° bis 90° (d). 


Korrelierten fiir Kugelzonen zwischen 9 = 45° bis 
75° und #=60° bis 75° wiedergegeben (Kurven a 
und b), wobei die Mikrophonachse mit der Kugel- 
polarachse iibereinstimmte. Die voll ausgezogene 
Kurve entspricht den theoretischen Werten fiir in- 
koharente Einstrahlung, die zugehérigen Me8punkte 
sind durch Kreise markiert. Als gerissene Linie ist 
die theoretische Kurve fiir koharente Erregung ein- 
getragen. Der untere Teil der Abbildung (Kurven ec 
und d) zeigt ahnliche Kurven fiir Zonen zwischen 
#=60° bis 90° und #=75° bis 90°. Die zur ko- 
harenten Erregung gehérenden Kurven aus Abb. 10 
sind zum besseren Vergleich hier nochmals eingetra- 
gen. In allen diesen Fallen erkennt man einmal die 
gute Ubereinstimmung zwischen Theorie und Experi- 
ment, insbesondere auch bei den Messungen mit 
phasenunabhangig betriebenen Lautsprechern. Fer- 


398 


5°6"O C50 0S 0100 G'0'0'60107 


4 ee A mk 
q0055090000000 0; 
Abb. 14. 


Autokorrelogramm fiir Anregung 
der Versuchshalle mit Rauschen. 
Terzband 2,85 bis 3,6 kHz. 


ner wird die Auswirkung der Multiplikation der fur 
inkoharente Erregung erhaltenen Korrelierten mit 
der Ortsfunktion sehr deutlich wiedergegeben, und 
damit sind die beiden Erregungsarten in ihrem Kor- 
relogramm deutlich unterscheidbar. Man sieht die 
Ubereinstimmung der Nullstellen und findet die Un- 
gleichung (6) fiir die Absolutbetrage der Funktions- 
werte sowohl durch die theoretischen Kurven als 
auch durch die Experimente bestatigt. 

Auch fiir die Ergebnisse reiner Autokorrelationen 
seien zur Erganzung noch zwei Beispiele gezeigt. Nach 
Gl. (17) sollte bei Verwendung von Rauschbandern 
das Ergebnis einer Autokorrelation z. B. des Schall- 
drucks am Aufpunkt P von den Einstrahlungsrich- 
tungen unabhangig sein. Es ist ohne weiteres ein- 
leuchtend, da Gl. (17) fiir das freie Schallfeld bzw. 
im reflexionsfreien Raum gelten muh, wenn die Kin- 
strahlung auf das Mikrophon nur aus einer Richtung 
erfolgt. Da aber auch bei Messungen in normalen 
Raumen mit grofen Diffusitaten diese Funktion ent- 
steht, ist aus Abb. 14 ersichtlich. Sie stellt durch die 
Kinhillende der aufgezeichneten Figur die Auto- 
korrelierte fiir stationare Anregung der grofen Ver- 
suchshalle mit Rauschen von Terzbreite (2,85 bis 
3,6 kHz) dar, aufgezeichnet im linearen Mafstab 
als Funktion der Laufzeitverschiebung t. Man er- 
kennt deutlich die von sin(4@ 1/2)/(M4m 71/2) um- 
hillte Funktion cos w,, 7. Die Ubereinstimmung mit 
den theoretischen Werten, die der Ubersichtlichkeit 
wegen nicht eingezeichnet wurden, ist gut. In Abb. 15 
ist schlieBlich das Ergebnis der Autokorrelation fiir 
ein mit einem Echo iiberlagertes Signal dargestellt. 
Dazu waren zwei Lautsprecher und ein Mikrophon 
im reflexionsfreien Raum kollinear aufgehangt, die 
Lautsprecher mit 390 cm Abstand voneinander. Sie 
wurden mit dem gleichen Signal betrieben, einem 
Rauschband zwischen den Grenzfrequenzen 2,8 und 
4 kHz. Die Aufzeichnung der Korrelierten erfolgte 
in diesem Falle ohne Beriicksichtigung der Vorzei- 
chen, daher sind auch die negativen Aste nach oben 
umgeklappt. Es erscheint bei etwa 11 ms, entspre- 
chend der Schall-Laufzeit zwischen den beiden Laut- 
sprechern, ein nach Gl. (20) verkleinertes Bild der 
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Abb. 15. 


Autokorrelogramm fiir Uberlagerung eines Signals mit einem Echo. 


Autokorrelierten” des von jedem der Lautsprecher 
abgestrahlten Signals, und man findet damit Gl. (19) 
gut bestatigt. 


Herrn Professor Dr. E. Myer danke ich fiir die 
Anregung zu dieser Arbeit und ihre standige Férde- 
rung. [hm verdanke ich die Bereitstellung vieler Mit- 
tel zur Durchfthrung der Untersuchungen und wert- 
volle Hinweise. 


Der Deutschen Forschungsgemeinschaft habe ich 


fiir zur Verftigung gestellte Sachmittel zu danken. 


(Eingegangen am 20. November 1956.) 
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The problems dealt with in this note arise in the 
use of electrical delay lines to deflect the beam of 


an 


THE RANDOMIZATION OF COMBINATIONS OF NOISE 
SOURCES BY MEANS OF DELAY LINES 


by D. G. Tucker 


Electrical Engineering Department, University of Birmingham 


Summary 

A discussion is given of the correlation between the outputs of the two ends of a delay 
line, or between the corresponding outputs of two different delay lines, when a number of 
narrow-band noise sources, identical except that their outputs are all uncorrelated, are 
connected to points uniformly distributed along the delay lines. It is shown that when the 
total phase-shift in a delay line has a suitable value (e. g.ta radians if the number of 
sources is large), the outputs at the two ends are quite uncorrelated, and that if the diffe- 
rences between the total phase-shifts of two delay lines has a suitable value (+ x radians 
as before, if the number of sources is large, but not the same as before if the number is 
small), then the outputs of the two lines are uncorrelated. The application of these results 
to the calculation of the noise factor of acoustic strip arrays is also discussed. 


Sommaire 

On étudie la corrélation entre les signaux de sortie aux deux extrémités d’une ligne de 
retard, ou entre les signaux de sortie correspondants de deux lignes de retard différentes, 
quand on connecte en des points distribués uniformément le long des lignes de retard un 
certain nombre de sources de bruit, a bande étroite et identiques, sauf qu'il n’y a pas de 
corrélation entre les signaux qu’elles émettent. On montre que si le déphasage total dans 
une ligne de retard a une valeur appropriée (par ex.: +z radians si le nombre de sources 
est grand), il n’y a pas du tout de corrélation entre les signaux de sortie aux deux extré- 
mités, et que si la différence entre les déphasages totaux des deux lignes de retard a une 
valeur appropriée (+ x radians comme ci-dessus, si le nombre de sources est grand, mais 
non comme ci-dessus, si le nombre de sources est petit), il n’y a alors pas de corrélation 
entre les signaux de sortie des deux lignes. On montre aussi l’application de ces résultats 
au calcul du facteur de bruit des réseaux acoustiques. 


Zusammenfassung 

Die Korrelation zwischen den Ausgangsfunktionen an den beiden Enden einer Verzége- 
rungsleitung oder den entsprechenden Ausgangsfunktionen an den Enden zweier verschie- 
dener Verzogerungsleitungen wird diskutiert fiir den Fall, daB eine Anzahl von schmal- 
bandigen Rauschquellen in aquidistanten Punkten an die Verzégerungsleitung angeschlos- 
sen ist. Die Rauschquellen sind einander gleich bis auf die Tatsache, daf ihre Spannungen 
unkorreliert sind. Es wird gezeigt, dali die Ausgangsfunktionen an den beiden Enden 
einer Verzdgerungsleitung unkorreliert sind, wenn die gesamte Phasenverschiebung auf der 
Leitung einen geeigneten Wert hat (z.B. + x falls die Anzahl der Quellen grof ist). 
Ebenso sind die Ausgangsfunktionen zweier Verzégerungsleitungen unkorreliert, falls die 
gesamte Phasenverschiebung der beiden Leitungen einen geeigneten Wert hat (+ x falls 
viele Quellen gegeben sind, ein anderer Wert bei Vorhandensein weniger Quellen). Die 
Anwendung dieser Ergebnisse auf die Berechnung der Rauschzahl von linearen Anordnun- 
gen akustischer Wandler wird diskutiert. 


1. Introduction delay line. With a suitable choice of the amount 
of phase-shift in the line, and assuming it to be 
constant over the range of frequency concern- 
ed, these two outputs may be made quite un- 


acoustic strip array or line transducer [1], but : 
* y i] correlated with one another. 


they are of some philosophical interest and may be 
of importance in other applications. They are: 


(a) 


Sources of noise 


If n independent (i.e. uncorrelated) noise 
sources of equal mean power and impedance, 
with identical energy spectra, are connected 
to points uniformly distributed along a delay 
line, as shown diagrammatically in Fig. 1, then 
two resultant noise outputs of equal R.MLS. 
voltage are available, one at each end of the Fig. 1. Noise sources connected to delay line. 


4.00 


(b) If the same noise sources as before are con- 
nected (through suitable decoupling circuits) 
simultaneously to two delay lines of different 
total phase shift, as shown in Fig. 2, then the 
resultant noise outputs at the corresponding 
ends of the two delay lines can be made quite 
uncorrelated by a suitable choice of the dif- 
ference between the phase shifts of the two 
lines. 


Sources of noise 


Decoupling 
ar at “capacitances 
Vp sate $ Delay line_No. 1 Vg 


Vig +— ps ae Delay line No. EDs 


Fig. 2. Noise sources connected to two delay lines. 


2. Analysis of the correlation between the two 
ends of the same delay line 


If the output voltages at the two ends of the delay 
line are V 4 and Vz respectively, then the correlation 
coefficient between them is given by 


VAV 
fi Se (1) 
Vee 


where the horizontal bar signifies ‘time-average’. This 
relation may be re-arranged thus 


a Wat Vs) 204 
oN eh 


(2) 


since clearly VV. It is thus necessary to con- 
sider the mathematical form of the sum of the two 
outputs. 

Let the sources be numbered zero in the centre, 
and +1, +2,...+7,...+m one side of the centre 
and —1, —2,...—r,...—m on the other side, so 
thatn=2m+1. 

Let the phase-shift be ® radians in each section 
of the delay line over the whole of the spectrum 
concerned. (Usually the spectrum is very narrow, 
and this condition holds very closely.) 

Let the source voltage of the r‘" source be repre- 
sented as 


WO: 
>» Pq cos(qt+4,), (3) 
qo, 


where @, and @, are the lower and higher limits of 
the spectrum in rad./s. 

The sum of the outputs (V4+Vg) contains two 
series of terms from this source: 


> vq cos[g t+ 49+ (m—r) P| (4) 
and >) v,q cos[g t+ag+(m+r) DP]. (5) 
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These combine to give 
2cosr P )'v;qcos(qt+a,+m®) , (6) 
(It is assumed that the delay line is free of loss.) 
Then (6) may be written as 
2cosr ®-V,’, (7) 


where V,’ has the same mean-square value as V,. 


Therefore the sum of the outputs of the two ends — 


of the delay line, counting all sources, is 
rTM 


Vat+Ve=) 2cosr®-V, , 


i> 


ces are all the same. (and equal to Vy say), the 
mean-square sum of V 4 and Vz is 


(V44+V5)? =Vx Y 2cosr B-V,', 


r=—m 


(9) 


(8) 


and since the mean-square values of the noise sour- — 


The mean-square value of V4 (or Vg) separately is 


Vi=nV*. 
Substituting now in equation (2) 


r=+m 


Ve De (2 cosr ®)? — a 
: 2nV%, % 


r=+m 


= dics rr 


¥ Seth P) 


n sin @ 


by a well-known trigonometrical transformation. 
Thus, when n is large, so that 0, but n @ re- 
mains finite, then 

_,sin n® sin®y 

Or ee Dy 

The relationship between @ and @,y for n large is 
plotted in Fig. 3. The outputs V4 and Vg are uncor- 
related when 0=0, i.e. when @p= ta, £27, ete. 


(13) 


Fig. 3. The relationship 0 _ sin Pp 


Dy * 


When n is not large, the relationship between @ 


(10) — 


(11) | 


(12) 


and @ is the repetitive pattern usually associated 


with diffraction effects. 


Clearly o=0 when 


n®—+x, +2zx, etc. But as there are only (n—1) 


sections in the delay line, so that p= (n—1) ©, 


then the outputs V 4 and Vz are uncorrelated when 


a Py = n—1 


wy lc. 


The assumptions made in the above working should 
be carefully observed. They are (1) that the phase- 
shift of the delay line is constant over the frequency 
_band of the noise, this implying that the bandwidth 
-is very small, and (2) that the time over which 

averages are taken in obtaining mean-square values 
_is much larger than the reciprocal of the bandwidth. 


_ 3. Analysis of the correlation between the 
- corresponding ends of different 
| delay lines 


Let the output V 4 from delay line No. 1 be V4; , 
and let the output V 4 from delay line No. 2 be V 42. 
Then following the lines of the previous analysis, 
the correlation coefficient between V 4; and V 42 is 


aS VanV ae _ (Vart Vaa)? ~2Vi41 


Vin Vie 20% 


The sum of the outputs (V4;+V 42) contains two 
series of terms from the r‘” source: 


: D vq cos[q t+ q+ (m—r) Py] 


j 


(14) 


(15) 
and 

Di rrqcos[qt+aq+(m—r) Pz], (16) 

where ®, and ®, are the phase-shifts per section in 
delay lines Nos. 1 and 2 respectively, and it is as- 
sumed that 2, + 2. 

_ There two series combine to give 

Oa ee 

en r)) 2 

Fi P,+ ®, 

x > v;q cos fata (m—r) rae (17) 


Thus, following the peevions analysis, 


r=+rm = 
| t eases [2 cos(m—r) - (Ps 5 x) “oy —2nVx 
8 Re aey 2n VR mae 
4 a (18) 


cos(m—r) (®,— M,) = 


j ey) See 29) 
icy  (G,-G,) +1 (19) 
sin ——*— 
2 
When n> w , (%,— ,) +0, while 
(2n— 1) (2,-¢ ; 2) mains finite, then 
~. sinn(P, — D,) sin (Py. — | Py) (20) 


n(P, — Dy—Dy, ’ 
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so that V 4; and V 42 (as also Vz; and Vgy) are un- 
correlated if ®y:— Opy = Ea, £27, ete. 

When n is not large, the effect is slightly different 
from that at the two ends of the same delay line. 
From eq. (20), 02=0 when 


sins Iya 
a i 
jf eee ey 
i 2 
i.e. when ein Se ?) i ye ela 2. oe TOES 


i. e. when 5 (©, — ®,) =n, 2n etc. 


n—] de Shay 
or 5 (®,—@,) = 3°? 9? ete. (21) 
The latter expression can be written as Py, — Py; 
=, 3%, etc., so that zero correlation is always ob- 
tained at these values independently of n. 


For example, when n=3, 0=0 when Py. — Dyy 


Ax 8x 


Sterner 3 3° 33, etc. 


4. Discussion: Noise factor of strip arrays 


Reference [1] gives an account of how delay 
lines can be used in conjunction with an acoustic 
(or, for that matter, electromagnetic) strip array 
or line transducer, not only to deflect beams, but 
also to enable special beamshapes (or directional 
patterns) to be synthesized. The transducer is di- 
vided up into sections along its length, and each 
section is connected to its corresponding tapping 
points on the delay line, so that Fig. 1 represents the 
arrangement if the noise sources are replaced by 
the transducer sections. It is shown in the paper 
cited that the process of synthesis consists of ad- 
ding together (with amplitude adjustments) the out- 
puts of one or more delay lines; if there are more 
than one, the arrangement becomes like Fig. 2. Al- 
though the transducer sections have replaced the 
noise sources of Figs. 1 and 2, yet they are, in fact, 
uncorrelated noise sources since noise arises in the 
effective resistance of the transducer sections, and 
noise from the medium of transmission is often un- 
correlated from one section to another. As the total 
phase-shift in the delay lines, assuming n is large, 
is always x radians or an integral multiple of x, the 
outputs of the lines are all uncorrelated as far as 
transducer noise is concerned. This makes the cal- 
culation of the noise factor of the array very simple, 
since all the noise outputs from the various delay 
lines can be added on a R.MLS. basis. 

In practice the synthesized directional patterns 
are usually achieved by adjustment of the sensiti- 
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vities of the various sections into which the array 
is divided, and the delay lines are used only ‘on 
paper’ as a step in the calculations. The noise factor 
can, therefore, be calculated directly by summing 
the R.M.S. noise voltages from the various sections 
of the transducer after the sensitivities have been 
adjusted, or, if the number of sections is large, the 
noise may be integrated according to the distribu- 
tion function of transducer sensitivity; but it will 
usually prove much quicker to add the noise outputs 
of the delay lines. A simple example was given in 
Section 8 of Reference [1], where the results of 
the present paper were used without proof. A fairly 
comprehensive study of signal/noise performance 
in strip arrays, where the present work is used, is 
given in a companion paper [2]. 
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List of symbols 


0 correlation coefficient between the two outputs V4_ 
and Vg of a delay line or between the correspond- | 
ing outputs of two different delay lines 

V4 output of one end of delay line 

Vg output of the other end of delay line 

n=2m-+1 number of noise sources 

@ phase-shift per section of delay line (radians) 

@y total phase-shift in delay line (radians). 

Vy R.M.S. voltage of noise source 

Subscript 1 refers to delay line No. 1 

Subscript 2 refers to delay line No. 2 

q angular frequency of typical component of noise 

dq phase angle of-typical component of noise 

ran integer, used for numbering the noise sources 

V, yoltage, as a time-function, of the rth noise source 
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SPANNUNGSOPTISCHE UNTERSUCHUNGEN 
AN BIEGEWELLEN MIT EXPONENTIELLEN. NAHFELDERN 


von O. Wets 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Das Auftreten yon exponentiellen Nahfeldern innerhalb des Geltungsbereichs der ele- 
mentaren Biegewellentheorie wird an Hand mehrerer Beispiele diskutiert. Einige dies- 
beziigliche spannungsoptische Aufnahmen werden vorgelegt. 


Summary 


The appearance of exponential near-fields within the limits of validity of elementary 
bending-wave theory is discussed with reference to several examples. Certain relevant 


stress-optical observations are quoted. 


Sommaire 


On examine sur quelques exemples la formation de champs exponentiels 4 courte distance 
dans le domaine de validité de la théorie élémentaire des ondes de flexion et on présente 


a ce sujet quelques clichés d’élasticimétrie. 


1. Einleitung 


Wegen der Vielfalt der in einem begrenzten Fest- 
k6rper anregbaren Wellentypen und ihrer teilweisen 
Konversion an Unstetigkeiten langs ihrer Ausbrei- 
tung gestalten sich Korperschalluntersuchungen in 
der Praxis meist recht schwierig. Das gilt besonders 
dann, wenn die Anteile der einzelnen Wellen ge- 
trennt erfaft werden sollen. 

Um die prinzipiellen Vorgange zu studieren, wird 
man zweckmafig am vereinfachten Modell Messun- 
gen anstellen [1]. 


MeBtechnisch ist es im allgemeinen iiblich, die 
Verschiebung s oder deren Zeitableitungen (Teil- 
chenschnelle 7= 08/02, Beschleunigung 978/90) 
auszunutzen. Leider sind diese Vektorfelder in ein- 
facher Weise nur an der Oberflache der Modelle zu- 
ganglich. 

In der vorliegenden Arbeit wird nun der Versuch 
unternommen, die in der Statik viel benutzte Span- 
nungsoptik auch bei der Schwingungsuntersuchung 
heranzuziehen. Dieses Verfahren erlaubt einen di- 
rekten Einblick in die Spannungsverteilung im In- 


‘neren eines Festkérpers. Das besondere Interesse 
_ gilt den exponentiellen Nahfeldern, die von Biege- 
_wellenuntersuchungen her bekannt sind. Sie treten 
an Storstellen wie Querschnittsspriingen, Ecken, bei 
~ Materialwechsel und dergleichen auf. 
Die Auswertung spannungsoptischer Aufnahmen 
gestaltet sich etwas kompliziert. Der tiefere Grund 
hierfiir ist in der Tensoreigenschaft der Spannun- 
_ gen. zu suchen. 


2. Zur Methode der Spannungsoptik 


Seit Davip Brewster weil man, daf alle durch- 
-sichtigen — urspriinglich optisch isotropen — Stoffe 
‘eine kiinstliche Doppelbrechung aufweisen, wenn 
_ man sie verzerrt oder verspannt. Auf diesem Effekt 
_basiert die Spannungsoptik. Sie wird fast ausschlief- 
lich in der technischen Statik angewandt. Wegen 
Einzelheiten sei auf das umfangreiche Schrifttum 
4 hingewiesen, z. B. [2], [3]. 

Die Spannungsdoppelbrechung wird genau wie 
in der Kristalloptik mit Hilfe von Polarisations- 

interferenzen beobachtet. 

| 


Bringt man ein Modell mit ebenem Spannungs- 
zustand in einen Zirkularpolarisator, so erhalt man 
das Interferenzliniensystem der Isochromaten: eine 
Art Hohenliniensystem fiir die maximale Schub- 
spannung Tp,x. Es gilt 


0,;—< Oxz — Oy \~ E 
Site y/(02 Pe) oat 


2 2 


Tmax = 


Dabei bedeuten 06,, 0, die Hauptspannungen in der 
_ay-Ebene, die senkrecht zum Lichtstrahl liegen soll, 
Orr, Syy, Try die kartesischen Spannungskomponen- 
ten. in der xy-Ebene. 

Im Linearpolarisator treten neben den Isochro- 
~maten noch die Isoklinen auf. Diese Interferenz- 
 linien dienen zur Ermittlung der Hauptspannungs- 
-_richtungen. 

Mit Hilfe der Polarisationsinterferenzen lassen 
“sich also bestimmen: 1. tax bzw. 0;— 0 aus den 
_Isochromaten, 2. die Richtungen von o, und oy aus 
den Isoklinen. Eine vollstindige Auswertung ver- 
 langt die Hauptspannungen o, und 9, einzeln, Da- 
zu ist noch eine dritte Messung nétig. Ist jedoch 
der Spannungszustand von der Theorie her bekannt, 
so geniigt der Vergleich mit den beobachteten Iso- 
chromatenbildern (Anwendung obiger Formel). [so- 
nenaufnahmen erginzen oft vorteilhaft die Iso- 


_ Ebene Spannungszustinde in der Statik konnen 
durch Angabe der Airy-Spannungsfunktion  voll- 
_stiindig beschrieben werden. Im dynamischen Fall 
existiert keine solche Spannungsfunktion. Die Ma- 
ialeigenschaften spielen hier die entscheidende 


¢ 


Rol 


. 


0. WEIS: UNTERSUCHUNGEN AN BIEGEWELLEN 


403 


3. Dynamisch-spannungsoptische 
Untersuchungsverfahren 


Die erste spannungsoptische Untersuchung an 
schwingenden Glasplatten stammt von Kunpr [4]. 
Im Drehspiegel konnte er periodische Aufhellungen 
beobachten. Macu erganzte diese Arbeit durch quan- 
titative Messungen beziiglich Schwingungsgesetz und 
Eigenfrequenzen. 1901 konnte Kénic [5] die zeit- 
liche Anderung der Isochromaten in einem Schwin- 
gungsbauch einer zu Biegeschwingungen angeregten 
Glasplatte photographieren. Er durchleuchtete dazu 
den zu untersuchenden Querschnitt intensiv (Licht- 
schlitz) und zog eine Photoplatte mit Hilfe eines 
Pendels daran vorbei. Tuzr und Nismpa [6] ersetzten 
schlieBlich das Pendel durch eine Trommelkamera. 
Dieses Prinzip wird immer dann verwendet, wenn 
der zeitliche Verlauf der Spannungen in einem be- 
stimmten Querschnitt interessiert. 

Ortliche Spannungsverteilungen lassen sich am 
besten durch kinematographische Aufnahmen oder 
durch stroboskopische Beobachtungen erfassen. Fiir 
einmalige Vorgange — wie Impulsausbreitungen, 
Bruchuntersuchungen — kommen nur kinematogra- 
phische Aufnahmen in Frage [7], [8]. 

Rein periodische Vorgange beobachtet man zweck- 
maf ig stroboskopisch. Von den Vorteilen seien nur 
angefiihrt: die beliebig lange Beobachtungszeit, ein 
verlangsamter Bewegungsablauf und die Moglichkeit 
der Aufnahme von Isochromaten und Isoklinen bei 
gleichem Schwingungszustand. Handelt es sich um 
stehende Wellen und interessiert nur die Frage nach 
dem vorliegenden Wellentyp, so geniigt die Beob- 
achtung im Gleichlicht. Beremann [9] hat eine Reihe 
solcher Aufnahmen von zu héheren Schwingungs- 
moden angeregten Glaszylindern verdffentlicht. 


4. Versuchsaufbau 


Die Moglichkeiten einer spannungsoptischen Un- 
tersuchung an elastischen Wellen sind in erster Linie 
durch das vorhandene spannungsoptische Material 
gegeben. Dieses soll bei den zur Verftigung stehen- 
den kleinen akustischen Leistungen eine mdglichst 
groBe Spannungsdoppelbrechung zeigen, anderer- 
seits aber durch die statische Beanspruchung (z. B. 
Eigengewicht) keine St6rungen im Isochromatenbild 
aufweisen. Weiterhin sind die auftretenden Material- 
verluste fiir die Wahl entscheidend. Ferner die Al- 
terung: Es treten Randspannungen auf, die auf 
Polymerisation oder Austrocknung zuriickgefiihrt 
werden kénnen. 

Benutzt wurde der auf ungesattigtem Polyester- 
harz basierende Kunststoff VP 1527 (Dynamit A.G., 
Troisdorf). Die Kunststoffmodelle wurden aus einer 
1,5 cm starken Platte geschnitten und auf genaues 
MaB gefrast und solange bei 90° C getempert, bis 


4.04. 


sie vollstandig spannungsfrei waren. Kigenspannun- 
gen machten sich bei Biegewellen vor allem an den 
Momentennullpunkten stérend bemerkbar. Sie tra- 


ten schon wenige Stunden nach dem Tempern wie-, 


der auf. Eine Isochromatenordnung (Gelblicht) ent- 
sprach bei VP 1527 rund 25 kp/cm? bei einer Mo- 
delldicke von lcm (statisch gemessen). Mehrere 
Ordnungen lassen sich am einfachsten tiber stehende 
Wellen erreichen. Das Isochromatenbild ist fiir fort- 
schreitende und stehende Wellen gleich, solange man 
geniigend weit von der Reflexionsstelle entfernt ist. 
Bei Untersuchungen an Querschnittsspriingen und 
an Ecken muBte allerdings im Stabteil hinter der 
Unstetigkeitsstelle eine rein fortschreitende Welle 
vorliegen, damit an der Sprungstelle definierte Ver- 
haltnisse bestanden. Da wegen der optischen Beob- 
achtungsmethode jegliche Reflexion sichtbar wurde, 
war eine moglichst grofe Reflexionsfreiheit in jenem 
Stabteil anzustreben. Versuche mit Sand als Damp- 
fungsmittel scheiterten an den zu grofen Biege- 
wellenamplituden. Auch ein Entdrohnungsmittel, auf 
eine Stabverlangerung von 4m aufgebracht, erwies 
sich bei einer Biegewellenlange von 50 cm als nicht 
ausreichend, Dagegen bewahrte sich folgendes Ver- 
fahren: Auf einer Strecke von etwa 2m wurde die 
Stabverlangerung aus Plexiglas mittels einer Heiz- 
wicklung auf iiber 100° C erwarmt, bis das Plexi- 
glas weich wurde. Die dann auftretenden wesentlich 
groBeren Verluste sorgten fir die notige Reflexions- 
freiheit. Der Stab war an Gummibandern aufge- 
hangt. Der Anfang der Verlangerungsleitung ist in 
Abb. 1 noch zu sehen. 

Die Polarisatoren waren Farbungsfilter von 25 em 
Durchmesser. Eine durch ein Lichtblitzstroboskop 


Abb. 1. Versuchsanordnung bei der Untersuchung an 
Stabecken. Zwischen den Polarisationsfolien 


befindet sich die Stabecke, dariiber das elektro- . 


dynamische Anregungssystem, nach links aus- 
laufend die Verlaingerungsleitung aus Plexi- 
glas. 
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beleuchtete Milchglasscheibe gab ein gleichmabig 
diffuses Licht. Als 2/4-Blattchen erwies sich 0,5 mm 
starkes Zelluloid, das durch den Herstellungsprozef 
meist eine gleichmaBige Doppelbrechung erhilt, als 
vollig ausreichend. Polarisationsfolie und 2/4-Blatt- 
chen waren fest zu einer Einheit zusammengefaBt. 
Je nachdem, welche Seite dieser Kombination zum 
Modell wies, wurde es von linear oder zirkular po- 
larisiertem Licht durchsetzt. 


Sinus- Kraft- 
generator verstarker system 
| Phasen= Impuls- Lichtblitz- 
schieber erzeugung stroboskop 


Abb. 2. Blockschaltbild der elektrischen Apparatur. 


Anregungs- 


Der elektrische Teil der Apparatur soll an Hand 
der Abb. 2 erlautert werden. Die von einem Sinus- 
generator erzeugte Spannung wurde einmal zur An- 
regung der Proben benutzt und zum anderen nach 
geeigneter Phasenschiebung zur Synchronisation des 
Stroboskops herangezogen. 

Die Isoklinenaufnahmen der Abb. 5 sind bei 
7,5 kHz aufgenommen. Hier konnte die phasenver- 
schobene Spannung einem Phasensummer entnom- 
men werden. Als Stroboskop diente eine helligkeits- 
gesteuerte Fernseh-Dia-Abtastrodhre (kurze Nach- 
leuchtzeit) , deren Leuchtfleck stark defokussiert war. 
Die Anregung geschah durch einen magnetostrik- 
tiven Schwinger. Alle anderen [sochromatenaufnah- 
men entstanden bei Anregungsfrequenzen unter 
350 Hz. 

Da sich Belichtungszeiten bis zu 10 Minuten er- 
gaben und teilweise mehrere Phasenlagen beriick- 
sichtigt wurden, muBte dafiir gesorgt werden, dah 
wahrend einer Zeitdauer von etwa zwei Stunden 
keine Anderung im Schwingungszustand auftrat. Das 
waren recht hohe Anforderungen, wenn man bertick- 
sichtigt, daB mit den groBten erzielbaren Anregun- 
gen gearbeitet wurde. Bei tiefen Frequenzen kam 
wegen der notigen Frequenzkonstanz nur ein RC- 
Generator in Betracht. Die Anregung besorgte ein 
kraftiges elektrodynamisches System. Das Lichtblitz- 
stroboskop besa eine Quecksilberhochdrucklampe, 
so dal} beinahe monochromatisches Licht vorlag. 


5. Spannungsoptische Beobachtungen 


5.1. Exponentielles Nahfeld bei der Reflexion einer 
Biegewelle am freien Stabende 


Der iibersichtlichste Fall fiir das Auftreten eines 
solchen Nahfeldes bietet sich bei der Reflexion einer 
Biegewelle am freien Stabende (Abb. 3). Man 


mu8 der einfallenden und reflektierten Biegewelle 


- 
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(~e*i** el’) ein vom freien Ende auslaufendes 


exponentielles Nahfeld (~e**e'’) iiberlagern, um 
dort die Randbedingungen der Spannungsfreiheit zu 
erfiillen: Biegemoment M,=0, Querkraft Q,—=0 
bzw. Normalspannung o,,=0, 
Tay 0 = 


Schubspannung 


exponentielles 
Nahfeld 
———_ —_———— —— 
einfallende und ae, 0 
reflektierte Biegewelle 


Abb. 3. Auftreten eines exponentiellen Nahfeldes bei 
der Biegewellenreflexion am freien Stabende. 
Die Amplituden sind normiert. Das Vorzeichen 
yon Orz, Vz bezieht sich auf y >0. 


Nun sind die gewohnlich benutzten vier Feldgré- 
Ben M., Q,, w. (Winkelschnelle) und v, = 97/9¢ 
(Transversalschnelle) gerade so definiert, daf in 
ihnen keine y-Abhangigkeit mehr auftritt. Fiir die 
spannungsoptische Untersuchung miissen wir jedoch 
auf die differentiellen GroBen zuriickgehen: 


M,=—fomydf, Qy=[tydf, we=-~”, 


0g7 = — 7! , 
I, Flachentragheitsmoment, ! Stabdicke. 

Abb. 4 zeigt die fiir //A=1/10 berechneten Iso- 
chromaten, Hauptspannungslinien und _Isoklinen. 
Die Abb.5 gibt die stroboskopischen Aufnahmen 
wieder. 

An Hand der Abb. 3, 4 und 5 soll nun erlautert 
werden, wie einige Eigenschaften der Biegewelle aus 
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werden koénnen. Die Spannungsfreiheit an dem 
oberen und unteren Stabrand (y= +//2) verlangt, 
daB dort t,, verschwindet. Die Hauptspannungs- 
linien miinden daher dort senkrecht ein oder laufen 
parallel. Am Rande kann nur o,, einen Beitrag zu 
Tmax liefern; die Isochromaten geben also ein direk- 
tes Abbild von o,,.. Man kann deutlich den sinus- 
formigen Verlauf samt Nahfeldanteil verfolgen. In 


Isochromaten 
Va=1V10 


Verschiebungen 


A/2 2/6 


Abb. 5. Spannungsoptische Aufnahmen von stehenden 
Biegewellen, die durch Reflexion am freien 
Stabende erzeugt wurden. 


den Stabquerschnitten, in denen t,,, Null ist, hat o,, 
seinen groBten Wert. Die Isochromatenverteilung 
iiber den Querschnitt ist dort aquidistant. Das ist 
eine Folge des Bernoulli-Geradliniengesetzes 0,,~ y. 
Im Querschnitt der Momentennull,,punkte“ ist nur 
die Schubspannung t,,, vorhanden. Lagen mehr Ord- 
nungen vor, so gabe die Isochromatenverteilung 
iiber den Querschnitt den parabolischen Verlauf von 
Tr, wieder. Im Isochromatenbild der Biegewelle fiir 
l/A=1/10 treten deutlich die singularen — _ hier 
spannungsfreien — Punkte hervor. Das ist einmal 
die zu einem Punkt zusammengeschrumpfte neutrale 
Faser, zum anderen sind es die Randpunkte im 
Querschnitt der Momentennullpunkte. Bei kleinem 
Verhaltnis //A erscheinen die Momentennullpunkte 
als gréBere schwarze Stellen, denn dann sind die 
Schubspannungen t,, im Verhaltnis zu o,, recht 


den spannungsoptischen Aufnahmen entnommen _ klein und bringen keine Aufhellung. 
Isochromaten i a/2 1234 5 AB 
_———————— 
Se ES seeeere ees | Reena eee 1 0 
Hauptspannungslinien Isoklinen 


Abb. 4. Berechnete Isochromaten, Hauptspannungslinien und Isoklinen fiir //Z=1/10 
(l Stabhohe, 4 Biegewellenlinge). Einige Isochromatenordnungen sind eingetragen. 
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exponentielle yy 
Nahfelder | 
x 
einfallende und durchgehende 
reflektierte Biegewelle Biegewelle 
Transveralschnelle Vy | Winkelschnelle wz is Querauslenkung 7 Biegemoment Mz Querkraft Qy 
t=to Wee ee Sy + maa NE we 
ner ea | 


w~N 


tote a aT TT ™® 


Abb. 6. Zusammenwirken von Biegewellen und exponentiellen Nahfeldern zur Erzwingung der Stetigkeit der 
FeldgroBen an einem Querschnittssprung 2 : 1. Das Verhaltnis der Feldgrofen vor und hinter dem Quer- 
schnittssprung ist von der Frequenz unabhiangig. 


5.2. Exponentielle Nahfelder am Querschnittssprung 


Die Biegewelleniibertragung an Querschnittsspriin- 
gen ist ausfthrlich theoretisch von Cremer [10] 
behandelt worden. Experimentelle Untersuchungen 
liegen ebenfalls vor [11]. 

Die Stetigkeitsforderung fiir alle vier FeldgroBen 
vy, Wz, VQ, und M, fihrt auf zwei exponentielle 
Nahfelder, die von der Stofstelle gegensinnig aus- 
laufen (Abb. 6). Das in negative x-Richtung 
fende Nahfeld und die reflektierte Biegewelle be- 
sitzen beztglich der einfallenden Welle unabhangig 


lau- 


vom Querschnittsverhaltnis immer die gleiche Pha- 
senlage; sie entspricht derjenigen bei der Reflexion 
am freien Ende, dem Grenzfall eines Querschnitts- 


Abb. 7. Isochromatenaufnahmen der Biegewelleniibertragung am Querschnittssprung 2 : 1. 
folgten im Abstande einer Achtel Periodendauer. 


sprunges. Ihre Amplituden hangen natirlich vom 
Querschnittsverhaltnis ab. Die durchgehende Biege- 
welle und das zugeh6rige Nahfeld sind ebenfalls in 
den Amplituden, aber auch in den Phasen von die- 
sem Verhaltnis abhangig. 

Beim Vergleich der spannungsoptischen Aufnah- 
men der Abb. 7 mit dem theoretisch zu erwartenden 
Feld (Abb. 6) ist zu beriicksichtigen, dai zwar Q,,. 
M, stetig sind, die auftretenden Spannungen jedoch 
wegen der verdnderten Querabmessungen in den 
beiden Stabteilen verschiedene Werte besitzen. Beim 
Vergleich der statischen Momenteniibertragung der 
Abb. 8 mit Abb. 7 zeigt sich, da am Querschnitts- 
sprung das gleiche Isochromatenmuster vorliegt; der 
Kinflu8 von t,, ist zu klein, um eine Abweichung 


Die Aufnahmen er- 
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Abb. 8. Statische Momenteniibertragung an verschieden 
gestalteten Querschnittsspriingen. 


zu bewirken. In der unmittelbaren Nahe des Quer- 
schnittssprunges ist der lineare bzw. parabolische 
Spannungsverlauf gestort. In einiger Entfernung 
sind jedoch nur die Resultierenden M,, Q,, wirksam, 
und damit liegt der normale Spannungszustand vor 
(Prinzip von Sr. Venant). Daher ist es gleichgiiltig, 
auf welche Weise der Querschnittssprung ausgebildet 
ist (Abb. 8). Es mu nur durch eine starre Verbin- 
dung fiir die Stetigkeit — mikroskopisch gesehen — 
aller FeldgroBen gesorgt werden. AuBerdem mu 
der Ubergangsbereich natiirlich sehr klein gegen- 
iiber den Biegewellenlangen sein. 


Abb. 9. Uberlagerung von Biege- 
wellen und exponentiellen 
Nahfeldern bei der reinen 
Momententibertragung an 
der symmetrischen Ecke. 
Tiefe Frequenzen: keine 
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- gangswiderstande fiir Dehnwellen — die durch die 


Querkrafte der Biegewellen vor und hinter der 
Ecke in Teil 1 und 2 angeregt werden — so grof, 
da die Transversalschnellen an der Ecke, die ja 
identisch mit den Schnellen der dort angeregten 
Dehnwellen sind, vernachlaissigt werden kénnen. 
Die Energie wird dann allein durch M, und w, 
ubertragen. Nun setzt sich der Energieflu8 


einer Biegewelle aus zwei Leistungsanteilen zusam- 


aber 


men: 
N,=Mzw2+Qyvy=— [ Ore v2 Af — f tay vy Af. 


Diese addieren sich zu einem ortlich und damit auch 
zeitlich konstanten Energieflufs V,. Da in der Ecke 
selbst die Transversalschnellen Null sein sollen, sor- 
gen zwei exponentielle Nahfelder daftir, da} der ort- 
lich und zeitlich konstante Energieflu8 in Teil 1 in 
der Ecke sin?-formig (Ubertragung durch M, w, 
allein) und in Teil 2 wieder konstant wird. 

In den Isochromatenbildern der Abb. 10 ist eine 
ortliche Storung des Spannungsverlaufs in der Ecke 
zu beobachten, die wiederum mit der statischen Mo- 


exponentielle 
Nahfelder 


—_ 


durchgehende Biegewelle 


{) einfallende und reflektierte 
Biegewelle 


oy 


Dehnwellenanregung. (Bei 
den mit * gekennzeichne- 
ten Phasen verschwinden 


die Nahfelder.) 


Querauslenkungen sessile] Winkelschnelle w, | Biegemoment He Querkraft Qy 


5.3. Exponentielle Nahfelder an der Ecke 


Eine eingehende theoretische Untersuchung der 
K6rperschalliibertragung an der Ecke ist ebenfalls 
yon Cremer [10] vorgenommen worden, Messun- 
gen an der Ecke wurden bereits durchgefiihrt und 
veroffentlicht [12]. 

In der Abb. 9 sind die Verhaltnisse an der sym- 
metrischen Ecke fiir tiefe Frequenzen (kleine ///- 
Werte) dargestellt. In’ diesem Fall sind die Ein- 


| 
> 
D 


mentenbeanspruchung tbereinstimmt. Die Ecke wird 
dabei 
macht sich die Dehnwellenanregung dadurch be- 
Momenten- 


etwas deformiert. Bei héheren Frequenzen 


merkbar, daS die Symmetrie an den 


nullpunkten stark gestort wird. 


Herrn Prof. Dr. E. 
Stelle fiir die 
die Moglichkeit ihrer Ausfiihrung herzlich danken. 


Meyer mochte ich an dieser 
Anregung zu dieser Arbeit und fiir 


(Eingegangen am 2. Marz 1957.) 
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Abb. 10. Isochromatenaufnahmen des Biegewelleneinfalls an der symmetrischen Ecke. Die 
Anordnung entspricht der Abb. 9. 
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